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The  analysis  was  based  on  reduction  to  a  type-0  follower  servo  equivalent  circuit  and  includes 
transient  behavior,  bandwidth  effects,  and  the  Applebaum  hard-limiter  modification. 

The  analysis  then  proceeded  to  a  ^-element  linear  array  with  K  adaptive  control  loops.  The 
system  analysis  consisted  of  a  Q-matrix  transformation  into  orthonormal  eigenvector  space,  and 
interpreting  the  transformation  in  terms  of  an  orthogonal  beam-forming  network,  similar  in  prin¬ 
ciple  to  a  Butler  matrix  network.  The  system  was  thus  converted  into  an  equivalent  “orthonormal,” 
adaptive  control-loop  network  to  which  the  type-0  follower  servo  analysis  can  again  be  applied. 

The  Q-matrix  transformation  network  produced  a  set  of  K  orthogonal,  normalized,  eigenvector 
beams  whose  output  power's  are  proportional  to  the  eigenvalues  of  the  covariance  matrix.  These 
beams  were  used  as  the  basis  of  a  convenient  expression  for  calculating  the  time-dependent  output 
pattern  function  for  the  array.  Performance  was  calculated  for  eight  different  distributions  of 
interference  sources  to  demonstrate  the  effects  of  power  level,  source  location  with  respect  to  the 
quiescent  beam  pattern,  source  spacing  in  terms  of  airay  resolution,  source  bandwidth,  and  con¬ 
tinuous  source  distributions.  The  Applebaum  hard-limiter  modification  was  also  introduced  into 
the  control  loops  and  the  performance  recalculated. 
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AN  INTRODUCTION  TO  ADAPTIVE  ARRAYS 


BACKGROUND 

Airborne  early  warning  (AEW)  radar  systems  operating  over  land  in  an  electronic 
countermeasures  (ECM)  environment  are  severely  troubled  by  both  jamming  and  clutter 

[1] .  Ground  clutter  can  be  very  large,  typically  some  30  dB  larger  than  return. from 
targets  of  interest*  and  because  of  aircraft  motion  (plus  antenna  rotation,  if  present) -its 
power  spectrum  generally  spreads  throughout  target  doppler  regions  of  interest,  partic¬ 
ularly  in  low-PRF  (pulse  repetition  frequency)  systems.  Clutter-cancellation  techniques 

[2]  used  in  present  AEW  radar  cannot  overcome  the  ground -clutter  problem,  and  a  per¬ 
formance  improvement  of  more  than  one  order  of  magnitude  is  essential. 

Adaptive  array  antenna  systems  arc  being  considered  as  a  solution  to  the  above  seri¬ 
ous  problems  because  they  can  respond  to  their  environment,  in  real  time,  in  both  the 
spatial  domain  and  the  time  (doppler)  domain  to  optimize  the  signal-to-noise  ratio  for 
the  system  [3].  Jammer  interference  is  sensed  in  the  spatial  domain,  and  array  pattern 
nulls  are  formed  in  the  directions  of  the  interference  sources  to  attenuate  their  contribu¬ 
tions  to  output  noise.  Platform  motion  is  sensed  in  the  doppler  domain  from  the  fre¬ 
quency  spreading  of  the  clutter  power,  and  the  array  shifts  into  appropriate  separated 
“phase  centers”  in  the  time  domain  to  collapse  the  clutter  spectral  spread  and  permit 
better  cancellation.  The  optimization  achieved  by  an  adaptive  array  may  also  result  in 
secondary  benefits,  including  compensation  for  antenna  pattern  distortions  caused  by 
aircraft  structure  blockage  and  scattering  effects;  compensation  for  element  excitation 
errors,  damaged  elements,  and  radome  effects;  the  possibility  of  more  arbitrary  element 
spacing  and  arrangement,  as  in  conformal  array-mounting  designs;  and  very  low  sidelobes 
for  antennas  on  aircraft  over  large  sectors  of  interest. 


1.  INTRODUCTION 

The  term  adaptive  array  has  been  applied  to  so  many  different  types  of  processing 
systems  that  it  is  well  to  define  the  class  of  systems  to  be  discussed  in  this  report.  For 
the  purposes  of  this  report,  an  adaptive  array  is  system  consisting  of  an  array  antenna  and 
a  real-time  adaptive  receiver-processor  which,  given  a  beam-steering  command,  samples  its 
current  environment  and  automatically  adjusts  itself  as  a  matched  filter  to  optimize  the 
ratio  of  signal  to  noise  (jamming,  clutter,  or  interference)  for  the  indicated  direction, 
frequency,  and  time.  Adjustment  control  is  accomplished  by  sensing  the  correlation  be¬ 
tween  element  signals,  i.e.,  on  the  basis  of  the  covariance  matrix  of  the  set  of  system 
inputs. 

Adaptive  arrays  as  applied  to  radar  systems  are  a  relatively  new  concept,  but  they 
have  roots  in  a  number  of  different  fields,  including  retrodirective  and  self-focusing  RF 
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antenna  arrays  [4],  sidelobe  cancellers  [5],  adaptive  filters  [6],  acoustic  or  sopar  arrays 
[7,8],  and  seismic  arrays  [9,10]. 

The  first  real  contribution  in  the  RF-antenna  field  was  the  retrodirective  array  in¬ 
vented  by  L.C.  Van  Atta  in  the  1950’s.  Another  major  step  was  the  development  of 
phase-lock  loop  theory  and  practice,  which  made  possible  self-steering  arrays.  Phase-lock 
loops  can  be  used  to  phase  coherently  add  signals  from  different  antenna  elements  in  an 
array.  A  further  advance  allowed  the  phase-lock  loop  scheme  to  produce  retrodirectivity; 
the  conjugate  phase  front  required  for  retrodirectivity  is  in  the  lower  coherent  sidebands. 
In  the  early  1960’s,  a  key  development  came  in  the  form  of  an  IF  sidelobe  canceller  cir¬ 
cuit  invented  by  Howells  [11].  This  type  of  circuit  has  been  widely  used  and  developed 
into  practical  phase-conjugate  adaptive  filtering  devices.  Howells,  Applebaum,  and  their 
coworkers  at  Syracuse  University  Research  Corporation  have  so  greatly  refined  and  gen¬ 
eralized  sidelobe  canceller  analysis  and  design  [12]  that  it  now  constitutes  one  oi  the 
most  important  contributions  to  adaptive  array  concepts. 

On  the  basis  of  the  sidelobe  canceller  experience,  Applebaum  [13,14]  developed  a 
control-law  theory  (algorithm)  for  adaptive  arrays  which  maximizes  a  generalized  signal- 
to-noise  (S/N)  ratio.  He  applied  the  theory  in  the  form  of  analog  adaptive  element- 
control  loops.  Because  of  its  practical  basis,  the  Applebaum  control  loop  was  selected 
as  the  model  for  discussion  in  this  report. 

Several  other  interesting  adaptive-array  algorithms  and  techniques  are  described  in 
the  literature.  Time  did  not  permit  their  inclusion  herein,  but  an  extensive  reference  list 
has  been  prepared  for  the  interested  reader. 

This  report  is  intended  to  be  tutorial,  and  starts  out  in  Sec.  2  by  reviewing  some 
basic  ideas,  namely  phase  conjugacy,  cross-correlation  interferometers,  and  the- IF  phase- 
cancellation  circuit.  Section  3  is  an  analysis  of  a  simple  two-element  array  that  has  a 
single  adaptive  loop  of  the  Applebaum  type.  One  can  get  an  excellent  perspective  on 
adaptive  system  performance  by  studying  this  single-loop  behavior,  since  it  is  easy  to 
keep  track  of  the  various  parameters. 

The  main  section  of  the  report  is  Sec.  4,  which  is  an  analysis  of  a  IC-element  linear 
array  with  K  adaptive  loops  of  the  Applebaum  type.  An  effort  is  made  to  interpret  the 
mathematics  in  terms  of  beam-forming  networks.  In  particular,  a  retrodirective-eigenvector- 
beam  concept  is  employed  to  aid  in  visualizing  the  transient  behavior  of  the  adaptive- 
array  output  pattern. 


2.  REVIEW  OF  SOME  BASIC  CONCEPTS 
2.1  Phase  Conjugacy 

The  adaptive  array  performs  spatial  filtering  by  sensing  automatically  the  direction 
of  a  source  of  interference  and  forming  a  retrodirective  receive  beam  in  that  direction  to 
subtract  from  its  normal  (unadapted)  pattern.  The  principle  is  illustrated  in  Fig.  1,  where 
“retrodirective  beam”  denotes  the  receive  beam  automatically  formed  in  the  direction  of 
a  single  source  of  interference.  To  achieve  retrodirectivity  [4,15],  the  phase  of  each  ele¬ 
ment  of  the  array  must  be  delayed  (with  respect  lo  a  given  reference  element)  by  exactly 
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Fig.  1— Retroditfictive  beam  principle  of  op¬ 
eration  for  adaptive  array,  with  one  interference 
source 


the  same  amount  that  the  incoming  wave  was  advanced.  At  any  given  frequency,  time 
delay  may  be  represented  by  a  phase  shift,  which  may  be  ambiguous  to  within  ±2n7r. 
Thusi  at  any  element  the  retrodirective  phase  must  bear  a  conjugate  relationship  to  the 
phase  of  the  received  signal  at  that  element  when  compared  to  a  common  reference 
element. 


It  turns  out  that  phase  conjugacy  can  be  obtained  very  easily  by  feeding  a  mixer 
with  a  reference  signal  either  equal  to  or  higher  in  frequency  than  the  received  signal, 
and  then  choosing  the  difference  frequency  as  the  output.  To  see  this,  let  us  briefly  re¬ 
view  mixer  operation  for  a  typical  hot-carrier  mixer  diode  with  two  input  signals,  as 
illustrated  in  Fig.  Z  Define  input  signals  E\  and  £2  as  the  real  parts  of  corresponding 
complex  numbers  £1  and  £2.  where 


Ei  =  «/(« !*+♦!) 


jr,  _  bei(u2t+<t>2) 


Ei  =  Re  (Ei)  =  a  cos  (o>if  +  1^1) 


£ 2  =  Re  (£2)  =  b  cos  (co2^  +  02)- 


The  instantaneous  current-voltage  relationship  for  a  hot-carrier  diode  is 


I  =  Is(eau  -  1) 


where 
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L277"  (*1  +  ^2)2  =  -£27P*  (£12  +  2*1*2  +  £22.)  • 


If  we  neglect  the  multiplying  factor,  the  product  of  £1  and  E2  in  Eq.  (2.9)  may  be 
written  as 


2(£j£2)  =  2fl5  cos  (wit  +  ^1)  cos  (co2t +  ^2)  (2.10) 

2(EiE2)  =.  ab  cos  [(o>i  ■>  C02K  -+  +  $2]  +  cos  Kwl  ~  w2)f  +  4>i  ~  02]  (2.11) 

2(£i£2)  =  Re(fiE2)  +  Re(Ei£|)  (2.12) 


where  is  the  complex  conjugate  of  E2  withu>2  <  <^i- 


Thus,  the  sum  frequency  term  corresponds  to  the  simple  vector  product  of  £1  and 
£2,  whereas  the  difference  frequency  term  results  in  a  product  with. the  complex  conjugate 
of  the  input  of  lowest  frequency.  Note  that  the  complex  conjugate  must  apply  to  the 
signal  of  lowect  frequency  to  obtain  the  correct  direction  of  vector  rotation,  i.e.,  a  “posi¬ 
tive”  difference  frequency.  The  amplitude  of  the  product  is  proportional  to  the  product 
of  the  amplitudes  of  the  two  inputs  but  is  also  dependent  on  the  mixer-diode_parameters 
arid' the  mixer  circuit,  as  indicated  in  Eq.  (2.9).  In  the  special  case  in  which  £j  and  £2 
are  of  the  same  carrier  frequency,  00 1  =  <*>2>  the  difference  frequency  is  zero  (dc  output) 
arid  the  complex  conjugate  may  be  .applied  to  either  signal.  ' 


2.2.  Cross-Correlation  Interferometer 


The  adaptive  array  derives  the  phase-conjugate  element  “weights”  for  forming  a 
retrodirective  receive  bearri  by  cross-correlating  the  received  element  signals  with  a  re¬ 
ceived  reference  signal.  The  reference  signal  may  consist  of  the  output  of  a  separate 
antenna  or  the  output  of  the  array  in  which  the  particular  element  is  located.  In  either 
case,  the  basic  principle  is.that  of  the  cross- correlation  interferometer.  This  is  a  basic 
type  of  antenna,  first  used  in  the  field  of  radio  astronomy.  It  is  shown  in  simple,  sche¬ 
matic  form  in  Fig.  3  [16,17] . 


Signals  from  a  single  point  source  at  angle  0  off  boresight  arrive  at  the  two  antennas 
A  and  B  with  a  path-length  phase  difference  of  u,  where 


2nD  .  a 
u  =  — r—  sin  6 
A 


(2.13) 


with 


D  =  distance,  between  antenna  phase  centers 
X  =  wavelength. 


These  RF  signals  are  translated  into  a  convenient  IF  band  by  two  mixers  fed  from  a  com¬ 
mon  local  oscillator,  so  that  the  RF  phases  and  amplitudes  are  preserved  in  the  IF  signals. 
If  v/e  denote  the  IF  signal  from  antenna  B  as  ^2,  we  can  write  the  real  part  as 


A'S'. . 


_ J Ife,  a  .■igafcS 


m 


WILLIAM  F.  GABRIEL 


1?2  =  b  cos  (w t  +  u) , 


(2.14) 


where  b  is  the  amplitude  of  the  signal  received  at  B,  w  is  the  IF  angular  frequency,  and 
u  is  the  phase  advance  at  B  in  reference  to  A.  The  IF  signal  from  the  mixer  of  antenna  A 
is  amplified  and  then  offset  to  a  higher  frequency  by  mixing  it  with  a  constant  reference 
offset  frequency  ojo  in  a  second  mixer.  (The  advantages  of  inserting  an  offset  frequency 
are  discussed  later.)  Thus,  we  may  express  the  real  part  of  signal  E\  as 


E\  =  a  cos  [(cj  +  woK  +  M  * 


(2.15) 


where  a  is  the  amplitude  of  the  signal  received  at  A,  (w  +  cjq)  is  its  shifted  IF  angular 
frequency,  and  <Pq  is  a  phase  constant  of  the  reference  offset. 


'  The  mixer  technique  shown  in  Fig.  3  is  only  one  of  several  methods  that  have  been 
used  for  achieving  a  frequency  offset.  Other  methods  include  the  following: 


1.  A  Fox-type,  0°-360°  RF  phase  shifter  may  be  continuously  rotated  [16] . 


2.  If  antenna  A  is  circularly  polarized  (as,  for  example,  a  cavity-backed  spiral 
antenna),  the  antenna  itself  can  be.  rotated,  as  in  the  AN/ASQ-96  DF  system. 


i, 
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3.  Various  electronic  serrodyne  techniques  may  be  used,  such  as  sawtooth-swept 
travelling  wave  tubes,  “staircase-driven”  digital  phase  shifters,  and  upper-sideband  ampli¬ 
tude  modulators. 

x 

The  continuous  phase-shift  rate  u>Qt  introduced  by  the  offset  frequency  in  any  of  the 
above  techniques  results  in  continuous  scanning  of  the  multiiobe  receive  pattern  of  the 
interferometer.  Lobes  scan  through  a  given  point-source  position  at  the  rate  of  (ojoM 
lobes  per  second. 

If  signals  £*  and  £2  are  combined  in  the  final  mixer  shown  in  Fig.  3  and  the 
difference-frequency  output  is  chosen,  then  from  the  mixer  description  given  in  the 
previous  section  the  output  will  be  the  real  part  of  the  product  of  £1  and  the  complex 
conjugate  of  £2: 

Re  [EXE%]  =  “Y  cos  (“0 *  +  $0  “  “)  (2.16) 

=  —  cos  u  cos  (cuo  t  +  <t> 0)  *  Y  sin  u  sin  (wo  t  +  <t> 0) . 

This  final  mixer  operation  is  a  true  cross-correlation  of  the  signals  received  at  an¬ 
tennas  A  and  B,  with  the  output  having  a  carrier  frequency  precisely  equal  to  the  refer¬ 
ence  offset  frequency.  The  amplitude  of  the  output  is  proportional  to  the  product  of 
the  signal  amplitudes  received  at  the  two  antennas.  The  phase  of  the  output  is  equal  to 
phase  difference  u  between  the  signals  received  at  the  two  antennas,  except  for  the  arbi¬ 
trary  phase  constant  4>q  associated  with  the  reference  offset  frequency.  The  constant  <Pq 
can  be  calibrated  out  in  a  standard  phase  detector,  permitting  measurement  of  phase  dif¬ 
ference  angle  u,  as  is  done  in  some  applications. 

The  cross-correlation  mixer  signal  is  then  amplified  in  a  high-gain  amplifier  of  gain 
G,  with  passband  centered  at  offset  frequency  wo,  and  integrated  in  a  narrowband, 
high-Q  filter  to  improve  the  output  signal-to-noise  ratio.  In  addition  to  reducing  the 
noise  in  the  output  by  narrowing  its  bandwidth,  the  filter  integrates  or  averages  the  cor¬ 
relation  mixer  output  signal.  This  averaging  is  denoted  symbolically  by  an  upper  overall 
bar: 


+  T  _  _  _ 

Re  [EXE%]  dt  =  Re  (EXE%) .  (2.17) 


The  averaging  is  very  important  in  adaptive  array  processing  and  is  discussed  in  greater 
detail  in  later  sections.  For  the  moment,  it  does  not  change  the  output  as  represented  by 
Eq.  (2.16)  because  we  are  dealing  with  simple  sinusoidal  signals  in  steady  state. 

Therefore,  the  cross-correlation  interferometer  receiver  system  shown  in  Fig.  3  results 
in  an  output  at  some  convenient  reference  offset  frequency  which  is  a  constant  complex 
number,  the  amplitude  of  which  is  proportional  to  the  product  of  the  amplitudes  of  the 
signals  received  at  the  two  antennas,  and  the  phase  of  which  is  equal  to  path-length  phase 
difference  u.  In  addition,  the  phase  of  this  complex  number  happens  to  be  the  conjugate 
phase  of  antenna  B  with  respect  to  antenna  A,  so  that  it  constitutes  the  exact  phase  angle 
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required  for  this  two-element  array  to  form  a  retrodirective  lobe  pointed  toward  the  signal 
source. 


It  is  not  necessary  to  inject  an  offset  frequency,  since  we  could  get  the  same  output 
result  from  a  “zero-IF”  correlation  detector  with  dc  output.  However,  the  dc  detector 
output  must  contend  with  all  of  the  other  dc  voltage  terms  represented'  in  Eqs.  (2.8)  and 
(2.9)  and  requires  balanced  mixer  networks  and  differential  dc  amplifier  balancing  in 
order  to  get  rid  of  the  unwanted  voltages.  This  is  generally  very  troublesome  in  low-level 
dc  detection  and  requires  continual  dc  drift  corrections.  Another  disadvantage  associated 
with  low-level  dc  detection  is  the  presence  of  1  If  noise  (flicker  noise). 


2.3.  IF  Phase-Cancellation  Mixer 


The  next  step  is  to  add  to  the  cross-correlation  interferometer  a  phase-cancellation 
mixer  as  shown  in  Fig.  4  (based  on  Ref.  11).  The  output  of  the  averaging  filter  is  de¬ 
noted  as  complex  weight  V/  and  its  steady-state  value  is  gain  G  times  the  averaged  cross- 
correlatov  output  from  Eq.  (2.16): 


W  =  kGiE^l) 

W  =  Re  (IF)  =  kab  G  cos  (co0t  +  Go  ~  u) , 


(2.18) 

(2.19) 


where  ft  is  a  constant  representing  the  correlation— mixer  conversion  factor.  The  constant 
ft  can  be  further  defined,  from  constants  given  in  Eq.  (2.9),  as 


(2.20) 


where  Rm  represents_the  mixer  circuit  load  resistance.  The  constant  ft  has  units  of 
volts-1,  such  that  k{ExE%)  will  have  the  expected  units  of  volts. 

Weight  W  and  signal  £2  are  fed  into  the  additional  mixer,  and  the  sum  frequency 
output,  denoted  by 


£3  =  k(WE2)  =  kZG&^&Ez, 


(2.21) 


is  chosen  so  that  the  product  of  E 2  and  the  averaged  cross-correlator  output  is  obtained. 
In  the  case  of  simple jinusoids,  as  represented  in  Eq.  (2.16),  there  is  no  envelope  modula¬ 
tion  to  average,  and  £2  can  be  directly  multiplied  into  its  own  complex  conjugate,  thus 
changing  Eq.  (2.21)  into 


E  3  =  ft2G|E2|2£. 


(2.22) 


£3  =  Re  (£3)  =  k2b2Ga  cos  [(w  *■  coq)*  +  M 


(2.23) 


The  net  result  is  that  the  phase  of  the  £2  sjgnal  is  automatically  cancelled,  and  the  out¬ 
put  is  precisely  in  phase  with  offset  signal  E\.  Also,  it  should  be  noted  that  the  ampli¬ 
tude  of  £3  is  proportional  to  the  received  element  power. 
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FI#.  4— Cross-correlation  interferometer  with 
IF  phase-cancellation  mixer  added 


In  the  general  case  where  E\  and  E%  are  broadband  envelope-modulated  signals,  phase 
:*  is  slightly  different  for  each  sr^tral  line  and  weight  IV  represents  an  average  conjugate 
phase  of  E<i  relative  to  £\.  Thus,  the  product  of  E%  with  IV  would  not  result  in  complete 
phase  cancellation  across  the  band,  and  the  output  would  be  slightly  different  from  TL\ . 


2.4  Integrating  Filter  Considerations 

Section  2.2  briefly  discussed  the  steady-state  output  of  the  integrating  filter,  bat  there 
are  other  asoects  to  examine.  For  convenience,  consider  the  single-tuned  RLC  filter  cir¬ 
cuit-  shown  in  Fig.  5.  This  can  be  analyzed  via  the  usual  Laplace-transform  approach  [18], 
s'-  ting  with  the  integrodifferential  equation  of  the  circuit, 


W  ,  r  dW  .1  f  W/1.  _  v 
S  *  c  It  *L  {Wdt -T 


(2.24) 


where  v  is  an  input  step-function  sinusoid,  and 
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V  =  0 


abG  . 


for 

t  < 

0 

for 

t  > 

0. 

(2.25) 

The  necessary  manipulations  are  contained  in  Appendix  A,  where  W  is  found  to  be 


IF  =  abG  [ e~atA  cos  (0f  +  \p)  +  B  sin  (coot  +  0)]  . 


(2.26) 


The  various  quantities  are  defined  in  Appendix  A.  Equation  (2.26)  consists  of  the  usual 
transient  term  plus  a  steady-state  term.  Examination  of  the  expressions  for  B  and  U  shows 
that  the  steady-state  term  will  not  be  identical  to  the  input  sinusoid  unless  the  circuit 
resonance  is  tuned  exactly  equal  to  coq,  whereupon  0  =  (0o  ~  u)  and  B  =  1/2.  Otherwise, 
the  filter  circuit  introduces  a  constant  phase  shift. 


The  transient  term  is  a  decaying  sinusoid  of  frequency  fi  nearly  equal  to  coo,  since 
we  assume  a  high-Q  filter,  and  the  exponential  decay  is  governed  by 


a  = 


2 RC  ’ 


(2.27) 


where  r  is  the  circuit  time  constant  or  equivalent  integration  time. 


In  adaptive  arrays  it  is  desirable  to  be  able  to  control  the  time  constant,  yet  not 
introduce  extra  phase  shifts,  which  would  interfere  with  proper  operation  of  the  phase- 
cancellation  mixer  shown  in  Fig.  4.  The  RLC  circuit  relationships  indicate  that  this  may 
be  difficult,  and  in  practice  this  has  been  found  to  be  the  case.  Thus,  tuned-carrier  inte¬ 
grating  filters  are  not  very  desirable  here.  A  compromise  solution  is  to  convert  down 
from  offset  frequency  wo  to  dc  baseband  I  and  Q  channels  (in  phase  and  quadrature)  at 
the  output  of  the  coo  amplifier.  The  integrating  filter  can  then  be  a  simple  RC  type. 
Such  a  compromise  is  shown  in  Fig.  6.  It  attempts  to  retain  the  “best  of  both  worlds”; 
i.e.,  it  retains  the  offset  frequency  output  from  the  cross-correlation  mixer  in  order  to 
keep  the  offset  detection  advantages  noted  in  Sec.  2.2,  and  after  amplification  it  converts 
to  high-level  dc  in  order  to  take  advantage  of  the  simpler  RC  integrating  filter  circuits, 
which  avoid  unwanted  phase-shift  problems  and  permit  easy  control  of  time  constants. 
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Fig.  6— Technique  for  using  simple  RC  integrator  filters 

The  arrangement  shown  in  Fig.  6  uses  the  I  and  Q  bipolar  video  from  the  RC  filters  to 
remodulate  the  offset  reference  signal  and  reestablish  weight  W  at  the  offset  frequency 
a>o  prior  to  mixing  with  E 2. 

For  the  rest  of  this  report,  it  is  assumed  that  the  integrating  filter  is  of  the  RC  type, 
as  in  Fig.  6,  and  our  filter  transient  analysis  therefore  can  be  based  on  the  simple  RC  cir¬ 
cuit  shown  in  Fig.  7.  The  differential  equation  will  be 


W  ,  n  dW  _  u_ 
R  r  L  dt  R 


T0  +  W  =  i>. 


(2.28) 


(2.29) 


where  r0  =  RC  is  the  circuit  time  constant  and  v  is  an  input  step-function  dc  voltage; 

v  -  0  for  t  <  0 

v  =  v0  for  t  >  0.  (2.30) 
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Fig.  7— Simple  RC  integrator  filter 


Taking  the  Laplace  transform  of  Eq.  (2.29),  we  get 


r0SW(S)  -  r0W(0)  +  W(S)  = 

(2.31) 

or 

_  V(Q)  <x0v o 

W{S)~{S  +  a0)  S(S  +  a  o)’ 

(2.32) 

where  Oq  =  l/^o.  and  W(0)  is  the  initial  value  of  voltage  W  at  time  t  ~  0. 
verse  Laplace  transform  results  in 

Taking  the  in- 

W  =  W(0)e"“°f  +  u0(l-e"“0') 

(2.33) 

or 

W  =  [W(0)  -  i>o]e'a°*  +  'Jo- 

(2.34) 

Since  we  could  do  this  separately  for  the  /  and  Q  filters,  it  is  obvious  that  Eq.  (2.34)  can 
be  written  in  terms  of  the  entire  complex  numbers  W  and  Vq,  so  that 

W  =  [W(0)  -  Uo]«fao'  +  v0  .  (2.35) 

This  gives  us  a  simple  transient  equation  for  complex  weight  W  representing  the  /  and  Q 
components  of  the  bipolar  video  coming  out  of  the  integrating  RC  filters. 


3.  TWO-ELEMENT  ARRAY  WITH  SINGLE  ADAPTIVE  LOOP 

The  first  adaptive  array  configuration  to  be  discussed  consists  of  a  simple  two-element 
array  with  a  single  Applebaum  loop  (Fig.  8).  It  is  schematically  similar  to  the  single-loop 
sidelobe  canceller  described  in  Refs.  11  and  19,  except  that  a  beam-steering  signal  is 
added,  as  described  in  Ref.  13.  The  most  significant  difference  between  Fig.  8  and  Fig.  4 
is  the  addition  of  a  summing  junction  in  which  is  added  to  with  the  sum 

then  becoming  the  output  and  being  fed  back  to  the  correlation  mixer.  The  arrangement 
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Fig.  8— Two-element  array  with  a  single  adaptive  loop 

is  intended  to  result  in  a  negative-feedback  servo  loop.  For  the  sake  of  simplicity,  Fig.  8 
does  not  show  the  local  oscillator,  IF  buffer  amplifiers,  and  bandpass  filters  of  an  actual 
system.  These  n  not  essential  to  the  analysis  at  this  point. 

Beam-steering  signals  B*  and  B\  are  intended  to  steer  the  receive  beam  in  some  de¬ 
sired  azimuth  direction  6q.  For  quiescent  conditions  wherein  only  receiver  noise  is  pres¬ 
ent,  adaptive  weight  W2  will  settle  to  a  steady-state  value  denoted  by  quiescent  weight 
Wq.  We  want  Wq  and  W\  to  be  precisely  equal  to  the  weight  values  needed  to  point  the 
beam  in  direction  6q.  Thus,  define  W*  and  Wq  as  those  desired  values  with  unit  ampli¬ 
tudes,  such  that 


W1  =  eju°  and  Wq  =  Wf  =  e~ju° 


where 


«o 


—  sin  Bq  . 
A 


(3.1) 


(3.2) 
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The  quantity  uo  is  defined  on  the  basis  of  a  reference  phase  point  midway  between  the 
ph»«e  centers  of  the  two  elements,  which  are  spaced  apart  by  distance  d.  These  weights 
result  in  quiescent  beam  pattern  Gq(0),  such  that 


Gn(6) 


-  ~\ei{-u 
2  . 


-«o)  +  e-J(u-u0) 


=  cos(u-u0), 


(3.3) 


where 


u  =  ~  sin  6  ; 

A 


(3.4) 


6  is  the  far-field  angle  variable. 

Weight  Wi  is  injected  directly  by  the  beam-steering  signal  B*,  so 

B*  =  Wi  '=  e;“° . 

Beam-steering  signal  8f,  however,  is  related  to  Wq  through  constant  62 »  as 


(3.5) 


B%  =  b2Wq  =  b2e~ju° 


(3.5) 


This  constant  will  be  evaluated  later.  B*  and  B^  are  injectrd  at  some  reference  offset 
<00,  but  the  eWot  carrier  term  is  not  included  since  it  would  be  only  an  extra  nonessential 
quantity  to  be  carried  along.  In  fact,  we  assume  in  this  discussion  that  all  signals  have 
bandpass  frequency  spectra  which  are  represented  by  their  complex  envelopes,  and  RF  or 
IF  carrier  modulated  by  that  envelope  will  be  assumed;  i.e.,  it  will  not  appear  explicitly 
in  any  of  the  equations  to  follow. 


Define  element  signals  Ei  and  E2  as  consisting  of  quiescent  receiver  channel  noise 
voltages  «i  and  n2  plus  a  statistically  independent  noisy  voltage  Jj  arising  from  a  single 
external  interference  source  located  at  angle  0,-.  It  is  assumed  that  the  source  of  inter¬ 
ference  is  narrowband  unless  otherwise  stated  and  that  the  source  is  suddenly  switched  on 
at  time  t  =  0  in  a  step-function  manner.  Thus, 


for  t  <  0  (quiescent) 


(3.6) 


Ei  =  ti\  +  Jfi  *Ul 
E2  =  n2  +  Jie^u‘  J 


for  t  >  0 


(3.7) 


where 


Uf  =  sin  8j , 


(3.8) 
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3.1.  Equations  for  Weights  W2  and  Wq2 

Weight  W2  is  equal  to  the  beam-steering  signal  3$  for  element  No.  2,  minus  the  out¬ 
put  V2  of  the  correlator: 

W2  =  Bt  -  V2 .  (3.9) 

The  correlator  consists  of  the  mixer,  amplifier,  and  filter.  Its  transient  behavior  may  be 
solved  for  on  the  basis  of  the  simple  RC  filter  approach  used  in  Sec.  2.4,  provided  that 
the  input  voltage  can  be  represented  as  a  suitable  step  function  over  time.  Now  the  out¬ 
put  from  the  correlation  mixer, 

*2  =  kHW^  +  W2E2)E% ,  (3.10) 

has  very  rapid  rms  fluctuations  in  accordance  with  the  receiver  channel  passband.  The 
filter,  however,'  integrates  or  averages  these  rapid  rms  fluctuations  in  accordance  with  its 
closed-loop  characteristics.  The  important  criterion  for  our  purposes  is  to  have  this  aver 
aged  value  x2  of  the  mixer  voltage  represent  an  input  step  function.  Thus,  assume  that 
x2  remains  constant  with  respect  to  time  except  for  the  step-function  change  occurring  at 
time  t  =  0,  which  implies  that  the  averaged  rms  values  of  the  element  signals  must  obey 
this  step-function  property  in  the  analysis  to  follow: 


x2  =  WWiEi  +  W2E2)E$ 

=  fe*Wi(f[£f)  +  k*W2\E2\2 

=  +  *2|f2|2(i,*_V2)>  (3.H) 

Voltage  V2  obeys  the  RC  filter  differential  equation,  Eq.  (2.27),  discussed  in  Sec.  2.4. 
Under  the  above  assumption  regarding  averaged  values,  this  may  be  written  as 

dV2 

T°  ~dT  +  Vz  =  G*2  *  (3-12) 

Substituting  Eq.  (3.11)  into  Eq.  (3.12)  and  rearranging  yields 


r0  +  (l+*2Gi£2l2)v2  =  k*G\E2\2 

For  a  step-function  change  in  (EiE%)  and  \E2\,  Eq.  (3.13)  yields  for  V2  the  solution 

V2  =  [Va(0  )-V2M)e~«t  +  V2(~),  (3.14) 

where 


ji’.giiiun'—iM'  MBarni 


B\ 


W^EZ) 


\e2\2 


(3.13) 
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and 


V2(°°) 


k2G\E2\2 

°2  +  -  9 

L  \e2\\  J 

1  +  k2G\E2\2 


(3.15) 


1  +  k2G\E2\2 


(3.16) 


The  quantity  1^(0)  is  the  initial  value  of  voltage  V2  at  time  t  =  0,  and  V2{°°)  is  the 
steady-state  value  after  the  transient  has  died  out.  This  solution  for  V2  is  substituted 
into  Eq.  (3.9)  to  get  weight  W2. 

The  quantity  Wq2  is  defined  as  the  optimum  value  of  weight  W2 ;  it  is  that  value 
which  minimizes  the  output  noise  power  of  the  array.  The  output  noise  power  is  the 
sum  of  the  quiescent  receiver  noise  plus  the  external  interference  noise,  weighted  by 
array  weights  and  W2.  If  Yn  is  the  array  output  noise  voltage,  the  expression  to  be 
minimized  is  the  mean  square  of  Yn,  or 

IVnl2  =  l(B W  +  (W2E2)  I2.  (3.17) 


We  see  by  inspection  that  the  optimum  value  of  W2  needed  to  minimize  Eq.  (3.17)  is 


^02  =  - 


(W^i  )Et 

\E2\ 2 


WiiEtEZ) 

\e2\2 


(3.18) 


Incorporating  the  assumptions  regarding  averaged  values  for  the  element  signals  results  in 


^02  =  " 


|i’2|2 


(3.19) 


Recalling  the  discussion  in  Sec.  2.2,  we  note  that  Wfo  is  the  normalized  retrodirective 
weight  which  will  place  a  perfect  spatial  pattern  null  in  the  direction  of  an  external 
source  of  interference. 

Returning  to  Eq.  (3.15),  we  note  that  the  relationship  for  W02  is  contained  therein, 
so  that  the  steady-state  output  of  the  correlator  may  be  rewritten  as 


V2M  = 


_  k2G\E2\Z{B2  -  Wp2) 


1  +  k2G\E2l2 


(3.20) 
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3.2.  Servo-loop  Considerations 

The  adaptive  loop  in  Fig.  8  exhibits  behavior  similar  to  that  of  a  Type  0  follower 
servo  [20],  a  simple  example  of  which  is  shown  in  Fig.  9a.  This  uni.,  feedback  circuit 
is  inherently  stable  and  could  represent,  for  example,  a  cathode  follower  or  an  emitter 
follower.  If  the  same  approach  is  used  as  in  Sec.  2.4,  the  equation  for  the  circuit  may 
be  written 

dvn 

T°  ~df+  v°  =  Me  =  “  Uo) 


ro  +  (l  +  M)uo  =  m- 


(3.21) 


For  a  Vi  input  step  function,  the  solution  for  vq  will  be  found  to  be 


■  [-<0.  -  (r^)] 


where 


1  +  p 


(3.22) 


(3.23) 


If  the  above  equations  are  compared  with  Eqs.  (3,13),  (3.14),  (3.16),  and  (3.20),  it  is 
obvious  that  there  are  direct  correspondences  between 

vo  and  V2  , 

H  and  k*G\E2\2 , 


Vi i  and  ( B2  ~  Wfo)  • 


Based  on  this  correspondence,  the  adaptive  loop  of  Fig.  8  can  be  modeled  by  the 
equivalent  servo  loop  shown  in  Fig.  9b,  whereupon  the  equations  associated  with  weight 
W2  simplify  to  the'  following  set: 


W2  =  B2  -  V2 


(3.24) 


V2  =  [V2(0)  -  V2(»)]e-«t  +  V2(oo) 


(3.25) 


viM  =  rrj<B2-  Ww) 


(3.26) 


1  +  M 
To 


(3.27) 
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v»°-*-+\  2  r—s —  ^ 


x' 


(a)  Type  0  follower  servo 


(b)  Equivalent  circuit  of  single  adaptive  loop 
Fig.  9— Servo  loop  schematic  diagrams 


W02  =  -  -ii-i-ii 

W 

y  =  fe2G|£2|2. 


(3.28) 


(3.29) 


This  equivalent  servo  concept  is  very  interesting  because  it  uses  the  optimum  weight  as  its 
input  signal,  and  the  degree  to  which  the  output  approaches  optimum  then  depends  on 
servo  gain  factor  y.  Note  that  for  y  »  1,  the  steady-state  value  of  W2.  is  equal  to  Wo2. 

It  is  also  important  to  note  that  the  equivalent  servo  gain  factor,  as  defined  in  Eq.  (3.29), 
is  proportional  to  the  average  power  level  from  element  No.  2.  Two  consequences  of  this 
peculiar  power-sensitive  gain  behavior  are  that  it  is  necessary  to  establish  a  minimum  value 
for  y,  based  on  quiescent  receiver  noise  power,  and  a  maximum  value,  based  on  the  maxi¬ 
mum  interference  power  to  be  received  at  the  elements. 

Addressing  the  minimum  condition  first,  we  would  have  E2  equal  to  quiescent  re¬ 
ceiver  noise  voltage  n2  at  the  mixers,  or 


minimum  y 


yo  =  k2G\n2\2 


(3.30) 
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By  choosing  an  appropriate  value  for  amplification  gain  G,  we  can  set  the  value  of  hq  as 
desired.  For  example,  a'  convenient  choice  might  be  unity.  Therefore,  the  gain  will  be 
defined,  on  the  basis  of  Eq.  (3.30),  as 


(3.31) 


g;  g 


m  i  I 


The  quantity  fe2!^!2  represents  the  voltage  coming  out  of  the  correlation  mixer  due 
to  n2  alone,  and  since  it  is  proportional  to  power  rather  than  voltage,  it  is  important  to 
use  enough  pveamplification  to  ensure  its  dominance  over  the  thermal  noise  voltage  gen¬ 
erated  at.that  mixer.  This  is  not  an  insignificant  problem,  as  can  be  see:,  by  examining 
the  two  magnitudes.  The  available  thermal  noise  power  at  the  mixer  will  be  KT0B0, 
where 

K  =  1.38  X10-23  J IK  (Boltzmann’s  constant), 

Tq  =  equivalent  noise  temperature  (K), 


B o  =  bandwidth  of  the  JRC  integrating  filter. 

If  we  choose  typical  values  of  725 K  for.  To  and  25  cycles  for  Bo,  then  the  available  ther¬ 
mal  noise  voltage  at  the  mixer  is 


/KTqBq  =  5X10” 10  volts. 


(3.32) 


6  $ 


h 

m- 


m  m 


By  comparison,  the  quiescent  receiver  noise  power  referred  to  the  elements  would  be 
KToBCt  where  Bc  is  the  element  signal  channel  bandwidth.  Choosing  5  Me  as  a  typical 
value  for  Bc,  we  have 

yjKToBc  =  \J% X  10"7  volts,  (3;33) 

and  n2  represents  this  voltage  after  preamplification  by  some  gain  A,  or 

n2  =  A  s/KtoBc  =  A  v/5  X  10"7  volts.  (3.34) 

Thus,  the  correlation  mixer  voltage  due  to  n<i  alone  will  be 

fc2 (WqE2E$)>  =  k*Wq\n2\2  =  k*WqA*  5  X  lO'*4  volts.  (3.35) 

Since  the  amplitude  of  Wq  is  unity,  the  magnitude  of  Eq.  (3.35)  is  essentially  (fe2|r»2l") . 
Furthermore,  the  mixer  conversion  factor  k,  given  in  Eq.  (2.20),  is  unlikely  to  exceed 
unity  by  an  appreciable  amount.  Therefore,  a  considerable  preamplification  gain  A  is 
necessary  to  ensure  the  dominance  of  fe2|n2l2  over  the  mixer  thermal  noise  level  as  repre¬ 
sented  in  Eq.  (3.32).  A  preamplification  gain  on  the  order  of  60  dB  is  indicated  here. 
Figure  8  shows  this  necessary  preamplification  in  the  schematic  form  of  a  single  preampli¬ 
fier  for  each  element.  The  preamplifiers  set  the  level  of  n2,  which  in  turn  determines  the 
servoamplifier  gain  G  from  Eq.  (3.31). 
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•  The  level  of  quiescent  noise  n2  also  determines  the  magnitude  of  beam-steeringsignal 
B$,  as  may*be'seeh  by  taking  the  steady-stateexpression  for  weight  W2,  given  by  Eqs; 
(3.24)  and  (3.26): 


W2{°°)  =  B$  ~  fT 7  ~  ^02)  •  (3.36) 

In  the  quiescent  condition,  Ei  =  ti\  and  £2  =  n2-  These  independent  noise  voltages 
have  zero  average  cross-correlation;  i.e., 

r  ' 

(nmZ)  =  0,.  (3,37) 

or  Hfo  is  zero.  Also,  the  quiescent  value  of  p  is  po,  so  that  the  quiescent  expression  for 
1^2 (°°)  becomes 


quiescent  (°°)  = 


*2 

1  +  Ho 


b2 

1  +  MO 


(3.38) 


The  quiescent  value  of  H^C00)  is  equal  by  definition  to  quiescent  weight  W0  defined  in 
Eq.  (3.1),  or 


Therefore 


f>2  =  1  +  Mo  • 


(3.40) 


The  amplitude  of  62  is  greater  than  the  magnitude  of  the  amplified  quiescent  voltage 
coming  out  of  the  correlation  mixer.  Note  that  the  phase  of  the  quiescent  correlation 
voltage  output  is  precisely  the  same  as  that  of  the  beam-steering  signal,  since  we  have, 
from  Eqs.  (3.37)  and  (3.38), 


Xq2  =  k^(WlE1+W2E2)E^  =  *S(Wri»i  +  W2n2)n$  =  fc2^l«2l2 


or 


quiescent  Xq2  =  ' k2\n2\2e~iu° .  (3.41) 

Then,  using  the  value -for  gain  G  from  Eq.  (3.31),  we  have,  for- quiescent  V2t 

-  Vq2  =  GXq 2  =  Moe~iu°  ,  (3.42) 

and  we  can  .doublecheck  the  quiescent  weight 

Wq  *  B%-Vq 2  =  [(1  +  Mo)e~iu°  -  Moe~iu°]  =  .e"y“°.  (3.43) 
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Next,  consider  the  effect  on  ju  of  a  steadily  increasing  interference  power  level  at  the 
elements,  in  order  to  establish  a  maximum  gain  condition.  From  Eqs.  (3.27)  and  (3.29) 
we  see  that  as  the  power  increases,  the  response/integration  time  decreases  and  produces 
less  and  less  averaging  effect  until  the  conditions  assumed  for  the  simple  transient  solu¬ 
tion  of  Eq.  (3.13)  are  no  longer  valid. 

Both  p  and  W<i  then  tend  to  follow  the  fast  fluctuations  in  the  envelopes  of  the  ele¬ 
ment  signals,  causing  weight  to  become  “noisy.”  Reference  21  presents  the  theory 
for  control-loop  noise  and  derives  expressions  for  the  variance  of  the  array  element  weights 
and  for  the  additional  noise  in  the  array  output  due  to  this  element  weight  fluctuation. 

The  upper-bound  condition  for  avoiding  “noisy  loops”  is  most  conveniently  stated 
in  terms  of  bandwidth;  i.e.,  the  closed-loop,  two-sided  bandwidth  should  not  exceed  ap¬ 
proximately  one-tenth  the  bandwidth  of  the  element  signal  channels  [19] .  This  ensures 
enough  integration  time  to  average  out  rapid  fluctuations  in  p  and  to  permit  IV2  to  be 
reasonably  independent  in  a  statistical  sense  from  the  instantaneous  fluctuations  of  the 
signal  envelopes.  For  the  simple  circuit  of  Figs.  9a  and  9b,  the  closed-loop  bandwidth  is 
a,  as  is  easily  seen  by  assuming  that  u/.and  uq  are  sinusoids,  whereupon  Eq.  (3.21)  can  be 
written  as 


joJTQVQ  +  Vq  =  H(V{-V  0) 


or 


.  m 

V°  1  +  p  *  /cotq 


(3.44) 


The  half-power  or  3:dB  bandwidth  point  occurs  for  the  condition 

CO3T0  =  (1+p)  or  CO3  =  a.  (3.45) 

Thus,  if  the  element  signal  channel  bandwidth  is  Bc,  then  the  upper-bound  condition  may 
be  expressed  as 


2  co  3  =  2 am  = 


27r.Bc 

10 


or 

_  irBcro 
f-'m  -  1Q 


(3.46) 


(3.47) 


where  pm  is  the  upper  bound  or  maximum  value -permitted  for  the  gain  factor. 

It  is  convenient  to  express  p  in  terms  of  po  and  a  power  ratio,  and  this  can  be  done 
by  noting  that  when  an  interference  source  is  present,  |E2I2  will  be  the  sum  of  the  squares 
of  the  magnitudes  of  n2  and  the -interference  source  voltage  Jt  at  the  mixers,  so  that 
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P  PGiEjl2  _  1»2|2  +  |J,12  _  f  Vi I2 
^0  "  fe2G|n2l2  ln2l2  l«2l2  ’ 

The  ratio  of  the  squares  of  the  voltage  magnitudes  is  equal  to  the  ratio  of  interfer¬ 
ence  power  to  receiver  noise  power.  If  this  power  ratio  is  denoted  by  Pi,  then,  from 
Eq.  (3.48),  m  can  be  expressed  in  the  form 


■M  =  Mo  (1+P/),  (3.49) 

where  the  gain  factor  is  essentially  normalized  to  the  quiescent  receiver  noise-power  level. 
One  can  proceed  to  calibrate  servo-loop  response  according  to  the  practical  ratio  of  inter¬ 
ference  power  to  receiver  noise  power.  Substituting  Eq.  (3.49)  into  Eq.  (3.47)  results  in 

TtBcTQ 

Pm  =  PoQ+Pim)  ~  m  -  1 


or 


T°  =  ^  +  +  PO^im)  •  (3.50) 

This  simple  expression  relates  basic  filter  time  constant  To  to  the  maximum  interference 
power  to  be  handled,  since  channel  bandwidth  Bc  is  generally  fixed  by  the  radar-pulse 
waveform  characteristics  and  cannot  readily  be  changed.  For  example,  choose  a  maximum 
interference  power  condition  40  dB  above  receiver  noise  level,  and  a  channel  bandwidth  of 
5  Me.  Then  from  Eq.  (3.50)  we  get,  for  Mo  =  1> 

To  =  mel  (2+1o4, 

or 

r0  =  0.00637  s. 


The  corresponding  filter  bandwidth  is  (1/2ttto),  or  about  25  cycles. 

The  relationship  between  servo-loop  gain  and  bandwidth  is  best  illustrated  by  making 
a  Bode  plot  [20]  as  shown  in  Fig.  10,  where  the  servo-loop  gain  is  plotted  as  a  function 
of  frequency.  Loop  gain  is  obtained  by  cutting  the  loop  at  a  convenient  point,  in  this 
case  the  feedback  path,  whereupon  we  get 


or 


juTovo  +  vo  =  pe 


loop  gain  = 


P 

1  +  iWT0  * 


(3.51) 


The  breakpoint  for  the  20-dB-pet-decade  slope  line  occurs  at  w  =  1/t0,  the  basic  RC  filter 
bandwidth  point,  and  the  intersection  of  the  slope  line  with  the  unity  gain  axis  occurs  ap¬ 
proximately  at  to  =  m/to-  The  slope-line  intersection  is  a  very  useful  indication  of  the 
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n=^Jcap*sBSBKas:aBtt^ 


Fig.  10— Bode  plot  illustrating  adaptive  loop  bandwidth 
variation  with  power  level 


overall  3-dB  bandwidth  of  the  servo,  which  can  be  verified  by  referring  to  Eq.  (3.45). 
Figure  10  shows  the  Bode  plot  for  interference  power  ratios  of  10  dB,  30  dB,  and  40  dB 
above  receiver  noise  level  and  is  obtained  by  substituting  Eq.  (3.49)  into  Eq.  (3.51)  to 
obtain 


loop  gain  = 


Mo(H-Pj) 
1  +  /wr0 


(3.52) 


3.3.  Revised  Equations  for  and  Hfo 

The  servo-loop  considerations  discussed  in  Sec.  3.2  have  introduced  so  many  modifi¬ 
cations  to  the  equations  developed  earlier  that  it  is  desirable  to  consolidate  all  of  these 
changes  into  a  revised  set  of  equations  for  W2  and  W02-  Starting  with  the  latter,  from 
Eq.  (3.19),  we  first  evaluate  the  averaged  cross-correlation  products  from  signals  Ei  and 
Ezt  as  defined  in  Eq.  (3.7): 


(E1E2)  =  \Ji\2e~i2u‘ 


(3.53) 
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(E2E%)  =  \E2\2  =  |n2|2  +  |./;|2 
Substituting  into  Eq.  (3.19)  results  in  the  convenient  expression, 

r 


r  Wi\Ji\2e-j2un 

L  i«2t2  +  Wi\2] 

=  -Wx  Pie~i2u‘. 


(3.54) 


(3.55) 


Next,  we  convert  the  transient  equations  in  V 2,  Eqs.  (3.25)  and  (3.26),  over  to  W2 
by  using  Eq.  (3.24);  i.e., 

W”)  - 1«5  -  v„~)i  =  [(i4t)  w* +  (rrir)  *«>]  <M6> 


W2(0)  =  flf  -  V2(0)  =  If,  «  Wf , 


(3.57) 


where  V<2(0)  is  the  quiescent  Vq2  evaluated  in  Eq.  (3.42).  Then  the  transient  equation  in 
W2  becomes 


where 


W2  =  [W2(0)  -  ^2(<»)]e-«f  +  1*2(00), 


1  +  M 

a  =  . —  -r  and  M  =  Mo(l  +?t)  ■ 


(3.58) 


(3.59) 


Note  that  steady-state  part  W72(°°)  has  two  distinct  components: 

Bo  / l  + Po\  * 
beam-steering  component  -  y ^  -j  W* 

and 

retrodirective  component  =  TJfj  ^02  • 

When  power  ratio  P,-  is  close  to  zero,  Wq2  is  also  close  to  zcro  and  the  beam-steering  com¬ 
ponent  is  dominant.  As  P;  increases,  \x  increases,  and  the  beam-steering  component  is 
attenuated  while  the  source  retrodirective  component  increases  in  magnitude.  For  large 
Pi  ratios  where  n  »  1,  the  beam-steering  component  becomes  negligible  and  the  source 
component  dominates,  essentially  “capturing”  adaptive  weight  W2. 
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3,4.  Adaptive  Loop  Performance 


initial  Conditions  and  Assumptions— Recall  from  Eqs.  (3.6)  and  (3.7)  that  quiescent 
receiver  noise  prevails  up  to  time  t  =  0,  when  a  single  narrowband  interference  source  is 
switched  on  in  a  step-function  manner.  Ratio  P,-  of  interference  power  to  receiver  noise 
power  is  assumed  to  have  a  maximum  limit  of  40  dB.  When  coupled  with  a  receiver 
channel  bandwidth  Bc  of  5  Me  and  a  minimum  servo  gain  factor  Mo  of  unity,  Eq.  (3.60) 
results  in  an  RC  time  constant  7o  of  6,370  /is.  Summarizing  these  conditions,  together 
with  G,  m.  JPi,  J?2»  and  1^(0),  we  have  • 


maximum  P/  =  P;m  =  10,000 


Bc  ■  5  Me 


Tq  =  6,370  ps 


Mo  =  1 


fe2!n2l2  fe2|n2l2 


M  =  Mo(l  +?i)  =  (1  +Pi) 


W1  =  eiu° 


B%  =  (1  +Mo)H^  =  2e"'“° 


W2(0)  =  Wf  -  e~ju°. 


Let  us  assume  that  the  elements  are  spaced  apart  by  a  half  wavelength  d  =  A/2,  to 


simplify  the 


expressions  for  no  in  Eq.  (3.2),  u  ir»  Eq.  (3.4),  and  u;  in  Eq.  (3.8). 


Transient  Behavior  of  Adaptive  Weight  The  first  performance  characteristic  of 
interest  is  the  transient  oehavior  of  the  adaptive  weight  l V2.  This  can  be  calculated  from 
Eqs.  (3.65)  through  (3.59),  using  the  mitral  conditions  and  assumptions  listed  above. 


The  transient  behavior  of  W2  depends  on  two  factors:  the  obvious  power  ratio  P,- 
contained  in  m  and  Hfo,  and  the  less  obvious  phase-angle  rotation  which  W2  undergoes  in 
reaching  its  adapted  value.  To  illustrate  the  latter  factor,  Figs,  11a  through  lid  con¬ 
tain  four  plots  of  the  amplitude  and  phase  of  vs  time,  with  the  same  power  ratio 
P,  =  100,  but  with  source  direction  0/  varied.  Starting  with  Fig.  lia,  where  the  source  is 
directly  on  boresight  (0,-  =  0  degrees),  W2  must  undergo  a  complete  180°  phase  reversal  to 
place  a  pattern  null  where  it  previously  had  a  maximum,  and  it  will  be  noted  that  it  does 
so  by  rapidly  decreasing  in  amplitude  to  zero  ip  about  45  ms,  abruptly  reversing  phase 
180°,  and  then  increasing  in  amplitude  toward  its  steady-state  value  of  0.96.  In  Fig.  lib 
the  source  is  10°  off  boresight,  which  requires  R'2  to  rotate  through  a  phase  angle  of  148°. 
For  this  rotation,  the  amplitude  drops  sharply  for  about  45  ms  but  does  not  decrease  to 
zero,  and  the  phase  changes  smoothly.  In  Fig.  11c  the  source  is  .30°  off  boresight,  which 
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requires  to  rotate  through  a  phase  angle  of  89°.  Note  that  the  amplitude  dip  has  now 
become  rather  shallow.  In  Pig.  lid  the  source  is  50°  off  boresight,  which  requires  W2  to 
rotate  through  a  phase  angle  of  only  42° ,  and  here  the  amplitude  varies  only  slightly  from 
unity.  If  the  source  is  moved  further  over  into  the  quiescent  pattern  null  at  90°  off  bore- 
sight,  V/2  is  already  correct  in  both  amplitude  and  phase,  so  that  no  transient  would  occur. 

The  effect  of  changing  the  power  ratio  amounts  to  changing  the  time  scale  of  the 
transient,  as  illustrated  in  Fig.  12,  where  we  have  the  same  conditions  as  in  Fig.  lib  ex¬ 
cept  that,  the  power  ratio  is  now  Pt-  =  1,000.  Note  that  the  transient  curves  are  almost 
identical  to  Fig.  lib,  but  the  time  scale  has  been  reduced  by  a  factor  of  10;  i.e.,  the  in¬ 
crease  in  source  power  level  by  a  factor  of  10  has  caused  the  adaptive  loop  to  respond 
ten  times  faster  in  changing  W2. 
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Fig.  12— Transient  behavior ’of  W2  for  Pj  *  1000,  6q  =  0°,  and  9t  =  10° 

Transient  Behavior  of  Spatial  Pattern— We  easily  obtain  the  behavior  of  the  spatial 
pattern  of  our  two-element  array  by  “freezing”  weight  W2  at  various  instants  of  time  dur¬ 
ing  its  transient  change  and  calculating  array  output  voltage  Y  produced  by  a  far-field 
test  source  of  variable  spatial  angle  0;  i.e., 


Y  =  S(Wie-Ju +  W2eiu),  (3.60) 

where  u  is  the  phase  factor  defined  in  Eq.  (3.4),  and  S  is  the  voltage  amplitude  at  the 
element  mixers  produced  by  the  test-  source.  As  a  practical  note,  it  would  be  necessary 
to  design  the  test-source  signal  and  receiver  system  so  that  the  test-source  signal  would 
not  be  sensed  by  the  adaptive  loop;  otherwise  it  might  influence  weight  W 2 .  The  spatial- 
array  pattern  is  proportional  to  the  absolute  value  of  Y,  and  for  convenience  in  plotting 
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we  generally  normalize  the  function  to  the  maximum  value  of  |  Yj.  Since  the  maximum 
value  for  two  elements  would  be  2 S,  the  normalized  spatial-array  pattern  function  G(8,  t) 
becomes 


G(6,  t) 


Y_ 
2  S 


=  -h-{Wie~iu  +  W2eiu): 
£  \ 


*  (3.61) 


Figure  13  contains  five  plots  of  |G(0),  t)|  under  the  condition  P,-  =  100  for  c  ^  0,  Qq  -  0°, 
and  0,-  =  15°.  The  plots  are  for  time  t  =  20,  50,  100,  and  200  ps,  and  for  a  steady  state. 
The  quiescent,  beam-steered  pattern  at  time  t  =  0  is  repeated  in  each  plot  to  serve  as  the 
beginning  reference.  Recall  that  this  quiescent  pattern,  denoted  as  Gq(8),  was  defined  in 
Eq.  (3.3).  It  may  be  verified  by  substitution  of  Eq.  (3.57)  into  Eq.  (3.61).  Thus,  start¬ 
ing  from  this  cosine  function  at  t  =  0  with  6q  =  0°,  the  pattern  plots  demonstrate  the 
progressive  development  of  a  pattern  null  in  the  direction  of  a  source  of  interference  at 
0;  =  15°. 


It  may  be  of  interest  to  note  that  the  final  Steady-state  pattern  can  also  be  manipu¬ 
lated  into  a  trigonometric  expression  if  the  steady-state  part  of  W2  from  Eq.  (3.56)  is 
substituted  into  Eq.  (3.61).  This  yields 


G(0,oo) 


=  cos  (u  -  Uo)  - 


Pi  cos  (u0  “  ) 

2  +  Pi 


=  G9(0)  - 


PiGq(8i) 

2  +  P{ 


gj(u-Ui) 


(8.62) 


where  Gq(0,)  is  the  value  of  the  quiescent,  beam-steered  pattern  function  in  direction  0,\ 
Note  that  for  the  particular  direction  0  =  0*,  G(0,°°)  reduces  to 


G(0f,  °°)  = 


2  Gq(0j)  2  cos  (u, --Mo) 

2  +  Pi  ~  2  +  Pi 


(3.63) 


which  gives  us  the  depth  of  the  pattern  null  in  the  direction  of  the  source  of  interference. 
For  the  conditions  in  Fig.  13,  the  depth  of  the  null  would  be  about  -35  dB,  which  seems 
remarkable  since  Pf  is  only  20  dB.  This  excellent  null  performance  occurs  because  the 
voltage  null  in  Eq.  (3.63)  is  inversely  proportional  to  power  ratio  P,-. 


•  Transient  Behavior  of  Output  Noise  Power— The  performance  factor  of  ultimate 

i  interest  in  an  adaptive  array  is  the  improvement  in  output  signal-to-nois  •  ratio  as  com¬ 

pared  to  a  conventional  array  subjected  to  the  same  interference  conditions.  In  this 
I  ratio,  the  signal  portion  can  be  readily  calculated  from  the  change  produced  in  G(0,  f), 

f  as  discussed  in  the  previous  section.  The  denominator  output  noise  portion,  however,  is 

(  more  fundamental  to  the  improvement  that  can  be  obtained  and  will  therefore  be  treated 

in  some  detail  in  this  section  and  the  one  to  follow. 

!  The  output  noise  power  is  the  sum  of  the  quiescent  receiver  noise  and  the  interfer- 

S  ence  noise,  weighted  by  array  weights  Wi  and  From  the  expressions  for  signals  Pi 

I  and  E2  as  given  in  Eq.  (3.7),  output  noise  voltage  Yn  may  be  written  as 

i 

I 

! 
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Fig.  13— Transient  behavior  of  spatial  patterns  for  Pj  •  100,  Bq  ■*  0°,  and  0;  *  16 
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Fig.  13  (Continued)— Transient  behavior  of  spatial  patterns  for  Pt  =  100, 0o  *■  0°,  and  0j  =■  16° 


Yn  =  WXEX  +  W2E2 


=  7l^l«ll2  +  \W2n2\2  +  Ji  [Wxe~jui  +  W2eiu)  , 


(3.64) 


where  the  square-root  term  expresses  the  fact  that  «i  and  n2  derive  from  independent  re¬ 
ceiver  noise  sources.  This  term  can  be  simplified  because  the  rms  amplitudes  are  assumed 
to  be  equal;  i.e., 


^nil2  +  |WV*2l2  =  \n2\Vm2  +  W2. 


(3.65) 


The  value  of  this  term  prior  to  time  t  =  0  represents  the  quiescent  output  receiver  noise 
voltage,  and  since  both  Wx  and  W2(0)  have  an  amplitude  of  unity, 


\n2\y/\Wx\2  +  irfLo  =  y/z\n2\. 


(3.66) 


Thus,  the  increase  in  output  noise  power  caused  when  an  interference  source  is 
turned  on  can  be  expressed  as  a  ratio  of  the. square  of  the  amplitude  of  Yn  to  the  square 
of  the  quiescent  output  receiver  noise  voltage;  i.e., 
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I 

I 


m 


m\ 


;  |2  • 


l. 


Mil 

l«2l? 


=  f  p  +  I W2!2 '  +  ^—5  \W1e-}ui>*  W2e'u<\ 


2\hi\ 

where  the  ratio  of  |«7/|2  to  |n2l2  can  be  replaced  by  power  ratio  P,-. 


(3.67) 


Figure  14  illustrates  the  plotting  ofEq.  (3.67),  when  06  =  0°  and  0,-  =  15° ,  for 
power  ratios  of  10  dB,  20  dB,  30  dB,  and  40  dB.  These  curves  show  clearly  the  direct 
dependence  on  power  ratio  P,-  of. the  speed  of  transient  response.  If  one  examines  the 
amount  of  time  elapsed  for  output  noise  power  to  drop  to  the  3-dB  increase  point,  the 
following  tabulated  are  obtained:  * 


Pi  Recovery  Time  to  3-dB  Increase 
(dB)  {{is) 


10- 

20 

30 

40 


1,000  ; 
175 
24 
3 


There 'are  two  additional  points  . of  interest  concerning  the  curves  in  Fig.  14.  First, 
the  starting  value  of  the.  output  noise  power  at  time  t  =  0  will,  of  course,  depend  on  the 
vaiue  of  the  array  space  factor  in  direction  0/;  i.e.,  the  location  of  the  interference  source 
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Fig.  14— Transient  behavior  of  output  noise  powe*  when  0O  =  0°,  0t  ■  15°, 
for  Pf  values  of  1.0  dB,  20  dB,  30  dB,  and  40  dB 
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in  the  quiescent  spatialpower  pattern  of-  the:  array.  For  our  tw;o-element  array,. this  is 
readily  shown  if  the  value  of  ^(0)  from  Eq.  (3.57,)is  inserted  into  Eq..  (3.67);  simpli¬ 
fying  shows  that  ~  •  *1 


JSy  =  1  +  2P,  cos*  (U0  -  Ui)  (3.68) 

!W‘i,  ... 

Note.that  the  squared  cosine  term  is  quiescent  power  pattern  value  |G9(0,-)|2. 

Second,  the  end  or  steady-state  value  of  the  increase  in  output  noise  power  after  the 
transient  has  died  out.  We  can  obtain  this  by  taking  the  steady-state  part  of  W2  as  given 
in  Eq.  (3.56)  and  substituting.it  into  Eq.  (3.67)  to  get 


{ 


1  +  |W2|2 


+ 


[4  y/Pj  COS  (Uq  “  «i) 


2  +  Pi 


f) 


(3.69) 


This  steady-state  equation  results  in.  unity  when  P,-  =  0  or  when  P,-  »  1,  but  it  is  greater 
than  unity  when  P/  is  near  unity.  This  steady-state  increase  in  output  noise  power  is 
shown  in  Fig.  15  plotted  vs  ratio  P,-  for  the  worst  condition,  where  0,-  =  0O.  Note  that  it 
peaks  at  an  increase  of  1.75  dB  when  P;  =  2,  and  that  steady-state  ^2(0°)  is  zero  for  that 
particular  ratio. 


Fig.  15— Steady-state  increase  in  output  noise  power  vsP( 
for  worst  condition,  0/  m  0q 
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Bandwidth  Effects— The  bandwidth  of  the  power  received  from  the  source  of  inter*, 
ference  will  affect  the  performance  of  our  adaptive  array  because  the  array  is  frequency 
sensitive.  Recall  that  adaptive  weight  W2  can  assume  only  one  amplitude  value  and  one 
phase  value  at  any  given  instant  of  time.  The  delay  (or  advance)  of  an  incoming  signal 
at  one  element  with  respect  to  another  involves  a  true  time-delay  distance,  as  shown  in 
Fig.  16,  where  time-delay  distance  ut  is  referenced  to  the  geometric  phase  center  of  the 
array.  For  any  given  wavelength  X  we  defined  ut  in  Eq.  (3.8)  as 


ird  .  0 
Uj  -  —  sm  Of. 


We  have  assumed  an  element  spacing  of  half  a  wavelength,  which  can  now  be  more  rigor¬ 
ously  defined  as  wavelength  A0  corresponding  to  center  frequency  f0  of  the  RF  bandwidth. 


Therefore  specify 


(3.70) 


JT  .  a 

ut  =T2sm9' 


=  —  Jt  • 


U*  ~  fo  2  ®in  * 


(3.71) 


ELEMENT- 


-  ELEMENT 


Fig.  16— Diagram  for  true  time-delay  distance  tq 


1  ' r' - ** 


a'A3±AXi,.Jts*r?Jki.<crr' 
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Fig.  17— Interference-source  poorer  spectrum,  uniform  amplitude 
lines  spaced  e  apart 


Frequency  f  can  be  further  defined  as  equal  to  f0  plus  some  offset  Af,  so  that  ut  finally 
becomes 


(3.72) 


With  this  frequency-dependent  relationship  for  u,-,  we  can  now  handle  an  interference 
source  with  bandwidth  by  dividing  its  power  into  a  number  of  discrete  spectral  lines,  the 
mth  line  of  which  has  associated  with  it  offset  frequency  Afm ,  voltage  at  the  mixers  Jm , 
and  power  ratio  Pm .  For  convenience  in  calculation,  assume  a  uniform  amplitude  spec¬ 
trum  of  lines  spaced  apart  by  a  constant  frequency  increment,  as  shown  in  Fig.  17.  Fur¬ 
ther  assume  that  the  lines  are  not  coherent  with  one  another;  i.e.,  they  will  r.ot  cross- 
conrelate.  Under  these  assumptions,  element  signals  E l  and  E2  >n  Eq-  (3.7)  may  be 
rewritten  as  the  summations 


' 


En 


nl  +  Z ]  Jm  e 


M 

Z 

m  » 1 
M 


[  for  t>  0, 


tin  + 


y  j- 


Jur 


(3.73) 
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where 


u"!  =  i1  +  ?r)i sin  6‘ 


(3.74) 


iff 

W2  -  £  IJ,,'" 


(3.75) 


M  =  total  number  of  spectrum  lines 
|«7/|  =  total  voltag.-.  agnitude  at  mixers. 


From  the  same  arguments  on  cross-correlation  as  used  previously,  the  cross-correlation 
products  from  the  new  E1  and  E2  would  be 


(3.76) 


in 

(¥2£|)  *  ln2|2  +  £  4 12  “  l«2l2  +  K-l2 


(3.77) 


The  new  expression  for  (EXE2)  k-.ust  be  incorporated  into  IF02,  changing  Eq.  (3.55)  to 


Wn  - - 


^1(5^2*) 


u  M 
I1 0  V1 


(3.78) 


The  expression  for  p  remains  the  same  as  before,  since  we  interpret  Pt  as  the  sum  of  all 
the  spectral-line  power  ratios;  i.e., 


M  =  M0(l  +  P() 


jn 

where  Pt  ~  P» 


It  is  evident  that  1F02,  although  optimum,  no  longer  represents  a  perfect  solution  for 
iV2  because  we  now  have  a  different  phase  angle  2 um  associated  with  each  spectral-line 
contribution,  and  we  will  get  a  single  resultant  vector  representing  the  sum  of  all  of  these 
small  vector  contributions.  Thus  W0  2  represents  an  adaptation  to  the  power  centroid  of 
the  interference  spectrum. 


The  increase  in  output  noise  power  can  be  expressed  by  a  ratio  similar  to  Eq.  (3.67): 
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-  +  \W2\2  +  £  Pm \(W1elu'»  +  W2efUm) |2j. 


(3.80) 


The  sum  term  consists  of  a  summation  of  individual  power  residues  which  correspond  to 
the  various  spectral  lines  of  the  source  of  interference.  Figure  18  shows  a  plot  of  Eq. 

(3.80)  vs  interference  source  bandwidth  when  P(  =  2,000,  6(  =  45°,  and  0O  =  0°.  Also 
plotted  are  the  phase  of  steady-state  W2  and  the  output  noise  that  would  result  from  the 
quiescent,  beam-steered  pattern  alone.  The  abscissa  scale  is  percentage  bandwidth,  and 
each  successive  point  represents  an  increase  of  one  more  spectral  line,  which  increases 
bandwidth  by  0.2%.  For  example,  starting  at  the  origin,  we  have  the  entire  power  ratio 
of  2,000  concentrated  into  a  single  spectral  line  of  “zero  bandwidth,”  and  steady-state 
adapted  weight  W2(°°)  from  Eq.  (3.52)  nulls  it  out  to  a  negligible  residue.  The  next 
abscissa  point  then  adds  a  spectral  line  at  a  frequency  0.2%  higher  than  the  first  line, 
giving  us  a  “bandwidth”  of  0.2%  for  the  two  lines  and  dividing  the  power  equally  into 
1,000  for  each  line.  The  next  point  adds  a  third  spectral  line  0.2%  higher  in  frequency 
than  the  second,  resulting  in  a  “bandwidth”  of  0.4%  for  the  three  lines  and  an  equal  power 
of  666  per  line.  This  process  of  adding  spectral  lines  is  continued  up  to  51  lines,  resulting 
in  a  bandwidth  of  10%  and  an  equal  power  of  39.2  per  line.  Figure  18  shows  clearly  the 
deterioration  (increase)  in  output  noise  power  as  the  constant  power  jf  the  interference 
source  is  spread  over  an  increasingly  wider  bandwidth. 

To  prevent  confusion,  recall  that  two  different  bandwidths  are  involved  here:  first,  the 
fixed  channel  bandwidth  Bc  of  5  Me,  which  determines  a  fixed  receiver  noise-power  level; 


444444444444444*4»44444«4*44 

444444444444' 


-UNADAPTED,  33 -dB  JAMMER 


.... 


W2  PHASE 


STEADY-.  'ATE  ADAPTED 


3  4 


8  9  10 


BANDWIDTH  IN  PERCENT 
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Fig.  19— Output  noise  power  contributed  by  each  of  the  61  spectral  lines  for  the 

10%  bandwidth  case 


second,  the  bandwidth  of  the  interference  source,  expressed  in  percentage,  which  must 
always  be  contained  within  Bc  in  order  to  be  meaningful. 

The  7-dB  increase  in  output  noise  power  at  10%  bandwidth  for  this  case  is  a  summed 
result.  It  is  of  interest  to  see  what  happens  to  each  spectral  line  in  tne  output  for  this 
10%  bandwidth  case,  and  Fig.  19  shows  the  increase  in  output  noise  power  contributed 
by  each  of  the  51  spectral  lines  within  its  own  0.2%  subbandwidth.  Here  we  see  the  ef¬ 
fect  of  adapting  W2  to  the  “power  centroid”  of  the  interference  spectrum;  i.e.,  the  center 
frequencies  of  the  power  spectrum  are  nulled  out  nicely  by  the  single  vector  weight,  and, 
even  though  the  nulling  must  deteriorate  for  frequencies  away  from  band  center,  the  ex¬ 
tremes  ate  well  balanced  to  achieve  the  best  overall  noise-output  compromise. 


3.5.  Hard-limiter  Modification 

The  adaptive  loop  configuration  of  Fig.  8.  which  we  have  been  discussing,  has  two 
drawbacks  because  voltage  output  X2  from  the  correlation  mixer  is  proportional  to  the 
power  received  at  the  elements: 

1.  The  dynamic  range  of  voltage  X2  is  the  square  of  the  element  signal  dynamic 
range.  Thus,  to  handle  the  40-dB  range  of  interference  power  discussed  in  previous  sec¬ 
tions,  the  correlator  branch  components  must  have  a  linear  dynamic  range  of  80  dB, 
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Fig.  20— Hard-limiter  modification  of  single  adaptive  control  loop 


which  is  uncomfortably  close  to  practical  component  limits.  Conversely,  a  severe  restric¬ 
tion  is  placed  on  the  system  input  signal  dynamic  range. 

2.  The  speed  of  response  is  proportional  to  power,  as  shown  in  Fig.  14,  and  this  re¬ 
sults  in  very  sluggish  response  for  weaker  sources  of  interference. 

The  best  solution  so  far  discovered  for  alleviating  these  drawbacks  has  been  to  modify  the 
adaptive  loop  by  incorporating  a  hard  limiter  in  the  conjugate  signal  branch  of  the  corre¬ 
lation  mixer  [19] .  This  hard-limiter  modification  is  shown  in  Fig.  20,  where  it  will  be 
noted  that  an  amplifier  of  gain  A2  precedes  the  limiter.  The  purpose  of  the  amplifier 
is  to  boost  the  amplitude  of  the  conjugate  signal  so  that,  after  limiting,  it  will  be  of 
normal  local  oscillator  level  when  fed  to  the  correlation  mixer.  This  amplitude  level  is 
denoted  by  constant;  h. 

The  modification  changes  the  equations  developed  in  the  previous  sections  because 
instead  of  the  correlator  product  as  given  in  Eq.  (3.10),  we  now  have 
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*i  -  ^OV,  +  >Mj>iit- 


(3.81) 


The  amplitude  variation  in  the  conjugate  signal  has  been  removed  and  only  the  phase 
variation  is  retained.  If  we  start  with  loop  gain  ju  in  Eq.  (3.29),  we  find  that  this  becomes 


UUn  _ 

M'  =  k*G'E21^  =  hk*G'\E2\, 


(3.82) 


so  that  loop  gain  is  now  proportional  to  voltage  level  rather  than  to  power  level.  If  we 
again  define  amplifier  gain  G'  as  in  Eq.  (3.31),  based  on  the  quiescent  loop  gain  of  h'q, 
the  new  gain  setting  will  be 


hk*\n2\ 


(3.83) 


There  is  no  longer  a  problem  in  getting  the  magnitude  of  the  voltage  coming  out  of  the 
correlation  mixer,  due  to  n2  alone,  to  dominate  over  the  thermal  noise  generated  there, 
because  enough  preamplification  has  been  inserted  to  raise  the  output  of  the  limiter  to 
local  oscillator  level  (on  the  order  of  1  mW). 

The  loop  gain  can  be  expressed  in  terms  of  a  power  ratio  as  in  Eq.  (3.49)  if  the 
above  two  equations  for  and  G'  are  used  together  with  Eq.  (3.54),  so  that 


u'  hk*G’\E2\  /|n8|»  +  |J,|» 

Mq  hk2G'\n2\  Y  ln2l2  ' 

Thus,  mVMq  is  simply  the  square  root  of  the  previous  expression. 


(3.84) 


Turning  next  to  the  relationship  between  maximum  interference  power  and  basic 
filter  time  constant  r0,  as  in  Eq.  (3.50),  we  have 


«;  -«iVi  +  ?,» - 


(3.85) 


If  we  again  take  the  maximum  interference  power  ratio  of  40  dB  and  a  channel  band¬ 
width  of  5  Me  for  Br,  then  the  value  for  Tq  is 


7q  »  0.637  X  10'6(1  +  100  m'o). 


(3.86) 


For  Mo  =  1,  Tq  =  64.3  ns,  and  the  corresponding  filter  bandwidth  is  2,475  cycles.  Note 
that  because  of  the  square  root  of  *lm>  the  basic  time  constant  is  faster  now  by  the 
factor  of  100. 

Optimum  weight  WQ2,  as  given  in  Eq.  (3.19),  remains  the  same,  but  when  expressed 
in  terms  of  new  servo  gain  factor  m\  it  must  change  to  the  new  expression 
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w2  w 


(3.87) 


The  equations  for  new  weight  W2  will  be  the  same  as  Eqs.  (3.56),  (3.58),  and  (3.59)  if 
ju'  is  substituted  for  y  hi  those  expressions;  i.e., 


W'2  =  [B^(0)  -  W'2(oo)]e-at  +  H^(°°) 


(3.88) 


(3.89) 


H^2v0)  =  B2*-  V2(0)  =  W* 


(3.90) 


where 


«'  =  and  y!  =  y^y/THT. 

To 

Note  that,  because  of  the  square  root  in  y',  a  considerable  change  has  occurred  in  both 
a'  and  W^f00)-  This  is  best  illustrated  if  Eq.  (3.67)  is  used  to  plot  the  transient  behavior 
of  the  increase  in  output  noise  power  for  the  same  conditions  as  in  Fig.  14.  The  new  set 
of  curves  for  the  hard-limiter  modification  is  shown  in  Fig.  21,  where  0O  =  0°  and  0,  = 

15°,  for  Pj  of  10  dB,  20  dB,  30  dB,  and  40  dB.  The  transient  decay  for  a  40-dB  ratio 
is  identical  to  the  response  in  Fig.  14  because  40  dB  was  the  maximum  power  condition 
for  both  cases  and  results  in  the  same  value  of  ct'm  =  am  =  1.58  X  106  s-1.  However, 
the  decay  for  ratios  less  than  40  dB  is  obviously  much  faster  with  the  limiter  modification. 
Another  large  difference  evident  in  Fig.  21  is  the  rather  high  steady-state  residue  (  approxi¬ 
mately  7  dB)  remaining  after  the  transient  term  has  died  cut.  This  is  principally  because 
the  new  loop  gain  y!  involves  the  square  root  of  P,-  and  is  therefore  much  smaller  than  the 
previous  y  for  the  same  power-ratio  values. 

To  reduce  the  unacceptably  large  output  noise  residue,  it  is  necessary  to  increase 
quiescent  loop  gain  y'Q  above  the  value  of  unity  chosen  previously.  Although  this  in¬ 
creases  the  value  of  Tq  almost  in  direct  proportion  to  y'o,  as  seen  in  Eq.  (3.85),  it  has 
only  a  minor  effect  on  a',  and  in  turn  on  the  speed  of  response.  Thus,  we  can  improve 
the  retrodirective  amplitude  portion  of  W^00),  by  increasing  y'0>  without  incurring  a 
penalty  in  response  time.  Figure  22  shows  the  increase  in  steady-state  output  noise 
power  residue  vs  Pt,  where  0,-  =  0O,  for  four  values  of  quiescent  loop  gain  y^  =  1,  3,  10, 
and  100.  This  shows  that  it  is  necessary  to  have  y'Q  >  10  to  hold  the  residue  to  a  reason¬ 
ably  small  value.  To  illustrate  that  the  response  time  is  unaffected,  Fig.  23  shows  the 
transient  behavior  of  output  noise  power  under  the  same  conditions  as  Fig.  21,  except 
that  now  y'Q  =  100.  If  h'q  can  be  thus  increased,  the  hard-limiter  modification  results  in 
a  satisfactory  output  noise  residue  and  retains  its  advantages  of  much  faster  response 
time,  double  the  dynamic  range  in  decibels,  and  better  loop  stability.  lTie  disadvantages 
include  the  inherent  limiter  problems  of  small  signal  suppression  when  more  than  one 
signal  is  received  and  the  possible  generation  of  spurious  response. 
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3.6.  Correlation  Coefficient 

In  the  study  and  analysis  of  adaptive  arrays,  it  is  desirable  to  become  familiar  with 
the  concept  of  a  correlation  coefficient.  This  can  be  done  by  starting  from  the  expression 
for  the  mean  square  of  output  noise  voltage  Yn,  as  given  by  Eq.  (3.17): 

\Yn\2  =  +  (W.2E2)  I2.  (3.91) 

If  the  expression  for  W02  in  Eq.  (3.19)  is  used,  \Yn  |2  can  be  manipulated  into  the  inter¬ 
esting  form 

I Y„\2  =  l«W2  -  (HV02i2l^2i2)  +  (1*2  "  ^02l2l^2l2)-  (3-92) 

Then,  for  W2  equal  to  P/02,  we  would  have  the  minimum  output  noise  residue  expressed 
as 
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or 


v  |2 
min 


I*. 


-  1 


I*. 


02  • 


\W,E 


v  IWi^l2 

Next,  define  a  coefficient  p  for  the  ratio  contained  in  the  brackets;  i.e.,  let 


Ipl2  = 


WQ2\2\E2\ 


\E\E2 


¥|2 


^il2|£2l 


or  define 


_  (E i£2*) 

p  ifiii®2r 


Thus,  we  have  a  normalized  cross-correlation  coefficient  for  element  signals  island  E2  or, 
simply,  a  correlation  coefficient  in  the  statistical  sense,  because  the  time-averaged  product 
(E[£f)  is  the  covariance  of  E\  and  E2,  and  it  is  being  divided  by  the  time-averaged  mag¬ 
nitudes  of  Ei  and  E2  which  are  the  square  roots  of  their  variances.  This  is  a  useful  con¬ 
cept,  because  the  real-life  signals  either  consist  of  noise  or  are  corrupted  by  noise,  so  that 
statistical  theory  is  needed  for  co.rect  mathematical  interpretation. 


Using  the  correlation  coefficient,  then,  makes  the  minimum  output  noise  residue 
simply 


\Y„\ 


2  _ 


(1  -  Ipi^lHW' 


(3.97) 


Note  that  (1  -  |p|2)  is  equal  to  the  ratio  of  minimum  output  noise  residue  to  the  power 
in  channel  1  and  therefore  defines  a  maximum  cancellation  ratio, 


\r. 


mm 


1*1*1 


=  (1  -  ipl2)- 


(3.98) 


Obviously,  a  good  cancellation  ratio  demands  a  high  degree  of  correlation  between  the 
two  channels. 


If  we  now  substitute  the  complete  expression  for  B^00)’  Eq.  (3.56),  into  Eq.  (3.92) 
to  incorporate  the  steady-state  servo  error  contribution,  we  find  that  the  output  noise 
residue  increases  to  the  value 


ir„l2  -  [i  -  IPI2  +  <!♦%> 


IV,* 


1*2 1 


w  |E,| 


11 


i  WM2. 


(3.99) 
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Fig.  24— An  analog  linear  adaptive  array  configuration  of  six  ele¬ 
ments  and  six  loops  with  beam  steering 


For  the  hard-limiter  modification,  n  must  be  replaced  by  y! . 

These  expressions,  based  on  the  correlation  coefficient,  permit  evaluation  of  decorre¬ 
lation  effects  caused  by  mismatches  between  the  transfer  functions  of  the  two  channels, 
differences  in  arrival  times  at  the  elements  (array  bandwidth  effects),  and  multipath 
effects.  Reference  19  is  recommended  to  the  reader  for  further  discussion  of  such  de- 
correlation  effects. 


4.  K-ELEMENT  ARRAY  WITH  K  ADAPTIVE  LOOPS 

Now  that  the  basic  principles  of  operation  and  the  performance  characteristics  of  a 
single  adaptive  servo  loop  have  been  discussed,  we  are  ready  to  consider  the  multiple-loop 
case,  wherein  each  element  of  a  IC-element  linear  array  has  an  associated  adaptive  servo 
loop.  A  possible  configuration  for  such  an  adaptive  array  is  shown  in  Fig.  24  for  six  ele¬ 
ments.  Note  that  each  servo  loop  is  arranged  in  the  same  manner  as  the  single  loop  of 
Fig.  8. 
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Define  element  signal  column  vector  E  in  which  feth  element  component  Ek  is  similar 
to  Eq.  (3.7 )  and  consists  of  quiescent  receiver  channel  noise  voltage  nk  plus  a  summation 
of  voltages  associated  with  /  external,  narrowband  interference  sources: 


where 


E'  =  [Ei.E2.E3,  ...,Ek) 


i=l 

(nd\  .  ,, 

■  (T;s,n  *<■ 


(4.1) 

(4.2) 

(4.3) 


Note  that  E'  is  the  transpose  of  matrix  E.  The  sources  are  assumed  to  be  statistically  in¬ 
dependent;  Jj  is  the  element  channel  voltage  amplitude  associated  with  the  ith  source,  and 
Oj  is  its  azimuth  angle  direction  from  array  boresight.  It  is  assumed  that  a  given  source 
will  induce  equal  voltage  amplitudes  at  all  of  the  array  elements.  Element  chase  is  refer¬ 
enced  to  the  geometric  center  of  the  array. 

The  beam-steering  signals  are  intended  to  set  up  a  shaped  receive  beam  which  is 
steered  in  some  desired  azimuth  direction  00.  For  quiescent  conditions  wherein  only  re¬ 
ceiver  noise  is  present,  the  adaptive  weights  will  settle  to  steady-state  values  denoted  by 
the  quiescent  weight  column  vector  Vfq .  We  want  the  components  of  this  vector  to  be 
precisely  equal  to  the  array  weights  one  would  choose  to  generate  the  desired  quiescent, 
shaped-beam  pattern  Gq(6).  Thus,  define  W?  as  this  desired  vector,  so  that 

y*q'  -  [WqVWq2,  wq3 . WqKy  (4.4) 


i 

I 

'i 

s 


where 


w  _  _  -ju0(2k-K-l) 

wqk  ake 

(4.5) 

(itd\ 

u°  \T/sin00, 

(4.6) 

The  values  of  element  coefficients  ak  are  chosen  to  achieve  the  desired  beamshape  and 
sitMob)  levelo.  The  quiescent  beam  pattern  can  then  be  expressed  as 

o,m .  ,s'w,)  - 1 

k  =  l 

(4.7) 

where  S  is  a  column  vector  representing  element  signals  of  unit  amplitude. 
u  is  associated  with  far-field  angle  variable  0;  i.e., 

Phase  factor 

S'  =  [Sx.S2.S3 . sk) 

(4.8) 
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where 
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Sk  =  ei“W-K-V 

ird  .  Q 
u  =  -jj-  sm  6. 


(4.10) 


The  components  of  input  beam-steering  column  vector  B*  are  related  to  the  components 
of  Vfq  by  constants  bk ,  so  that 


r>*/  -  TD*  D*  »* 

a  >  ^3  ,  . . .  9  u§ 


where 


K  -  bkWqk. 


(4.11) 


(4.12) 


The  6a  will  be  evaluated  in  the  following  section,  where  it  is  shown  that  they  are  deter¬ 
mined  by  the  quiescent  servo  gain  factor. 

Voltage  outputs  from  the  correlator  filters  are  represented  by  column  vector  V,  so 
that  we  can  define  an  adaptive  weight  column  vector  VV  for  the  configuration  of  Fig.  24: 


W  =  W3  =  [B*  -  V]  =  B3*  -  V3 


(4.13) 


where  Vk  is  the  voltage  output  from  the  kth  correlator  filter. 

Analysis  of  multiloop  adaptive  arrays  requires  some  familiarity  with  the  theory  of 
linear  differential  equations,  matrix  algebra,  and  the  solution  of  eigenvalue  problems. 
References  22  and  23  are  recommended  for  readers  who  desire  to  review  the  mathematics. 


4.1.  Adaptive  Weight  Equations 

The  adaptive  servo  loops  in  Fig.  24  are  the  same  as  the  single  loop  in  Fig.  8,  which 
was  described  in  Sec.  3,  so  that  formulation  of  the  adaptive  weight  equations  can  proceed 
in' much  the  same  manner.  Thus,  weight  Wk  associated  with  the  feth  clement  is  equal  to 
beam-steering  signal  B*  minus  output  Vk  of  the  associated  correlator  filter,  so  that 


Wt  =  B*  -  VL. 


(4.14) 


Hie  averaged  r  orrelation  mixer  voltage  Xk  is  agmn  the  averaged  product  of  element  signal 
Ek  with  the  summed  output  of  the  array, 
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%  =  *2fel>f*/i. 


(4.15) 


and  voitage  Vk  will  obey  the  same  RC  filter  differential  equation,  Eq.  (3.12)  and  will 
result  in  the  similar  averaged  value  expression, 


dVk  /  *  A 

r»  df+  v*=ykLw& 


where 


7  »  fe2G. 


(4.16) 


(4.17) 


The  quantity  y  is  defined  as  a  conversion-factor  gain  constant,  assumed  to  be  the  same 
for  each  of  the  servo  loops.  However,  whereas  Eq.  (3.12)  contains  only  one  unknown 
and  can  be  solved  immediately,  the  above  equation  contains  K  unknowns  and  must  be 
handled  as  a  member  of  a  set  of  K  simultaneous  linear  differential  equations,  one  for  each 
adaptive  loop.  Foi  convenience,  let  us  use  Eq.  (4.14)  to  convert  from  Vk  to  Wk,  noting 
that 


dt  dt 


(4.18) 


whereupon  Eq.  (4.16)  may  be  rewritten  in  terms  of  V/k  as 


'o  nr'**  -  K  -  7(3 


't  "wY 


(4.19) 


The  complete  set  of  these  weight  equations  for  subscript  ft  values  of  1  through  K  may 
then  be  expressed,  in  terms  of  the  more  convenient  matrix  notation,  as 


r0  lit  +  W  =  B*  "  7[E*W'E] 

where  W'  is  the  transpose  of  matrix  W.  Recall  from  matrix  multiplication  that 


(4.20) 


A 

(W'E)  =  (EfW)  =  £  WjEi , 


(4.21) 


so  that  the  product  of  the  three  matrixes  in  Eq.  (4.20)  may  also  be  written  as 

[E*W'E]  =  [E*E'W]  =  [E*E'JW. 


(4.22) 
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lhe,  averaged  product  of  il»e  conjugate  of  signal,  column  vector  E*  and  its  transpose  results 
in  a  matrix  whose  components  represent  the  correlations  between  the  various  element 
channel  signals  or,  in  other  words,  the  covariance  matrix  of  the  set  of  system  inputs.  De¬ 
fine  this  covariance  matrix  as  M; 


M  =  [E*E'J  = 


F*  P  P*  F  P*P 


!»♦  P 

av 2  />2  "2  ^3 

p#  p  p*  p  p^p 

"3"1  "3  "2  A3A3 


^  ^  ”*J 


(4.23) 


Note  that  the  components  of  matrix  M  may  be  obtained  from  the  channel  signals  defined 
in  Eq.  (4.2),  whereupon  these  averaged  correlations  are  found  to  be 


<??)  -  £  Wl2e'2“‘(,'‘> 

i-1 


for  /  =/=  k 


(4.24) 


or 


(^)  *  |£*  |2  =  (In* l2  +  J  | //  for  /  =  ft, 

/-l 


(4.25) 


where  Eq.  (4.25)  represents  the  diagonal  elements  of  M.  It  is  evident  that  the  covariance 
matrix  may  be  written  as  the  sum  of  quiescent  receiver  noise  matrix  plus  individual 
interference  source  matrixes  M,-,  such  that  H 


where 


and 


/ 

M  =  Mg  +  £> 
**1 


Mo  - 


I«i  I2  0  0 

0  (n2 12  0 

0  0  |n3|2 
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(4.27) 
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M;  =  \T,\2 


1 

el2u‘ 

efAu‘ 

-j2ut 

1 

el2Ui 

-/Au, 

e-i2Ul 

1 

(4.28) 


Note  th?v  M  is  a  positive,  definite,  Hermitian  matrix. 

Substituting  M  into  Eq.  (4.20)  and  rearranging,  we  obtain  a  concise  final  form  of  the 
adaptive-weight  matrix  equation,  , 


dW 


T0  ~  +  [I  +  7M1W  =  B  , 
where  I  is  the  identity  matrix,  such  that 


I  =  [By]  = 


10  0 
0  10 
0  0  1 


(4.29) 


(4.30) 


The  solution  of  Eq.  (4.29)  can  be  accomplished  via  a  special  transformation  which  con¬ 
sists  of  the  eigenvectors  of  the  covariance  matrix.  From  the  theory  of  matrixes,  we  know 
that  a  positive,  definite,  Hermitian  matrix  such  as  M  can  be  diagonalized  by  a  nonsingular, 
or+honormal,  modal  matrix  transformation  which  shall  be  defined  as  matrix  Q.  Further¬ 
more,  we  know  that  the  resulting  diagonal  components  are  the  eigenvalues  of  matrix  M. 

In  accoiUance  with  the  usual  eigenvalue  problem  statements, 

|M  -  &2I|  =  0  and  Me,  =  ft2e,.  (4.31) 

the  /3,2  are  the  eigenvalues  (real,  positive  numbers)  of  M,  and  e,  are  the  associated 
eigenvector!,, 


'ell 

ei2 
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These  eigenvectors,  which  are  normalized  to  unit  Hermitian  length  and  are  orthogonal  to 
one  another,  make  up  the  rows  of  the  Q  matrix: 


ell  el2  el3  '  ‘  ‘ 
e21  e22  e23 

Q  =  e31  232  e33  ‘  ‘  • 


(4.33) 


.  ek  1  ek  2  ek3 


Diagonalization  of  M  by  the  Q-matrix  transformation  is  expressed  in  the  form 


[Q*MQ‘]  =  [(?,%•]  = 


0l2  0  0 

0  022  0 


0  0  03 


2  ... 


(4.34) 


Since  M  is  a  product  of  the  conjugate  of  signal  vector  E  and  its  transpose,  Eq.  (4.23),  it 
is  evident  :nat  Eq.  (4.34)  may  be  written  as 

[Q*MQf]  =  [Q*E*E'Qf]  =  [E*Ef] 


where 


[E*Ef]  =  [ft2$y] 


E  =  QE. 


(4.35) 


(4.36) 


Thus,  the  Q  matrix  transforms  real  signal  vector  E  into  new  orthonormal  signal  vector  E. 
The  components  of  E  are  determined  by  the  eigenvectors  of  M;  i.e.,  from  Eq.  (4.33), 


Sk  -  (efc'E). 


(4.37) 


Note  that  these  components  have  two  special  characteristics.  They  are  decorrelated, 
so  that 


(EkE{)  =  0  for  l*k, 

and  their  amplitudes  are  the  square  root  of  the  eigenvalues,  so  that 


(4.38) 
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{E*E,)  =  ^2  for  /  =  k 


i4l  *  h- 


(4.39) 


The  transformation  is  equivalent  to  the  operation  of  an  orthogonal  beam-forming  network, 
and  this  analogy  is  pursued  further  in  Sec.  4.3  in  terms  of  physical  antenna  feed  networks. 

The  transformation  of  the  E  vector  in  Eq.  (4.36)  suggests  that  a  similar  transforma¬ 
tion  may  be  performed  on  input  beam-steering  vector  B*  defined  in  Eq.  (4.12),  since 
this  vector  can  be  viewed  as  the  conjugate  of  an  equivalent  signal  vector  B.  Therefore, 
define  transformed  beam-steering  vector  B*  as 


B  =  QB 


B*  «  Q*B>( 


(4.40) 


where  the  ftth  component  is  determined  by  the  kVn  eigenvector, 

B*  =  (e*'B*). 


(4.41) 


The  Q-transformation  operations  on  both  E  and  B*  suggest  an  equivalent  circuit 
representation  for  the  system  illustrated  in  Fig.  24,  in  which  Q-transformation  networks 
would  be  used  to  achieve  an  “orthonormal  adaptive  array”  system.  This  new  equivalent 
circuit  representation  is  shown  in  Fig.  25,  alongside  a  simplified  schematic  diagram  of  the 
real  system.  There  will  be  a  new  set  of  weights  W  in  the  orthonormal  system,  and  if  one 
develops  the  adaptive  weight  matrix  equation  in  a  manner  similar  to  Eq.  (4.29),  it  will 
be  found  that 


where 


T0  [I  +  tM]W  =  B* 


M  =  [E*E']  =  \fit%f]. 


(4.42) 


(4.43) 


Thus,  in  the  orthonormal  system  we  obtain  a  set  of  independent  linear  differential  equa¬ 
tions,  each  of  which  has  a  solution  if  the  eigenvalues  can  ^e  determined.  Each  of  the 
orthonormal  servo  loops  will  behave  as  if  the  other  loops  did  not  exist,  because  the  Ek 
signals  are  orthogonalized  and  have  zero  correlation.  The  feth  servo  loop,  then,  can  be 
viewed  in  much  the  same  manner  as  the  single  servo  loop  discussed  in  Sec.  3,  and  from 
Eq.  (4.42)  its  weight  equation  will  be 


+  a +  yPkWk  =  b;. 


(4.44) 


The  first  important  similarity  to  note  is  the  equivalent  servo  gain  factor,  which  may  be 
defined  from  the  above  equation  as 


-  Tftfc2 


(4.45) 
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Fig.  25— Equivalent  circuit  representations  for 
six-element  adaptive  array  configuration 
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so  that  our  equivalent  servo  gain  factors  are  now  determined  by  the  eigenvalues  of  the 
covariance  matrix.  When  compared  against  Eq.  (3.29),  note  that  positive,  real  eigenvalues 
/3j.2  correspond  to  the  square  of  a  signal  voltage  amplitude.  This  analogy  is  carried 
further  in  Sec.  4.3,  where  it  is  shown  that  any  given  eigenvalue  is  proportional  to  the 
power  appearing  at  its  orthonormal  network  output  port. 

There  is  also  an  optimum  weight  Wok,  which  will  be  derived  in  Sec.  4.2  and  found 
to  be  equal  to 

*.*  “  (%)*»•  (4  46) 


If  we  substitute  nk  and  W0 k ,  into  Eq.  (4.44),  it  can  be  rewritten  in  the  form 


+  d  +  =  nkWok. 


(4.47) 


This  equation  is  now  of  exactly  the  same  form  as  Eq.  (3.21),  and  for  a  step-function 
change  in  the  input  signal  a  similar  solution  may  be  written: 


where 


Wk  -  [^*(0)  -  #*(°°)]e‘“*'  +  #*(«») 


(4.48) 


(4.49) 


(4.50) 


Also,  W^(<»)  is  the  steady-state  weight,  Wk(0)  is  the  initial  weight  value  at  f=0,  and  ak 
is  the  transient  decay  factor.  The  transient  responses  will  now  be  determined  by  the 
eigenvalues.  This  feth  orthonormal  servo  loop  may  be  modeled  as  the  simple  Type-0 
follower  servo  shown  in  Fig.  26. 


M: 


Fig  26— Type-zero  follower  servo  model  for  the  rth  ortho¬ 
normal  adaptive  control  loop 
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The  next  step  is  to  relate  these  Wk  solutions  to  real  weights  Wk.  Note  that  the  out¬ 
puts  of  the  two  systems  shown  in  Fig.  25  must  be  identical;  i.e., 

L  -  £  «A 

*"i  *  =  1  (4.51) 

W'E  =  WfE  =  WfQE. 

Thus, 

W'  =  W'Q  or  W  =  Q'W,  (4.52) 

and  the  solution  for  the  ftth  real  weight  becomes 

Wk  ~  («l^l  +  e2fc^2  +  e3k  ^3  +  *  ’  ’  +  eK  k  ^ K  )  •  (4.53) 

Each  real  weight  is  therefore  a  summation  of  all  the  orthonormal  weights. 

Using  Eq.  (4.52),  we  can  go  back  to  Eq.  (4.29)  to  verify  the  orthonormal  Eq.  (4.42). 
Substituting  for  W  in  Eq.  (4.29)  results  in 

TojjjW'Vt]  +  IQ'  +  7MQ']W  =  B*.  (4.54) 

Multiplying  through  on  the  left  by  Q*,  we  get 

r0^[Q*Q'W]  +  [Q*Qf  +  7Q*MQ']W  =  q*b*.  (4.55) 


I  Because  of  the  orthonormal  properties  of  the  eigenvectors  in  the  Q  matrix, 

?  [Q*Q'J  =  l,  (4.56) 


and  by  Eq.  (4.34)  the  covariance  matrix  will  be  diagonalized  by  the  Q  transformations, 
which  results  in 


v0^  +  [(1  +  7ft2)5//]W  =  B*. 

Next  evaluate  the  bk  coefficients  which  relate  B*  to  Vfq  (Eq.  (4.12)).  To  do  this  in 
a  simple  manner,  assume  quiescent  conditions  in  which  only  receiver  noise  is  present,  so 
that  the  element  channel  signals  are  decorrelated  and  there  is  no  difference  between  trans¬ 
formed  weights  and  real  weights;  i.e.,  the  quiescent  Q  matrix  would  be  an  identity  matrix, 

quiescent  Q  =  Qq  -  I.  (4.58) 

The  quiescent  covariance  matrix  from  Eq.  (4.27)  is  already  diagonalized,  and  if  we 
further  assume  that  the  receiver  noise  power  in  all  element  channels  is  equal  and  denoted 
by  |n0|2,  then  from  Eq.  (4.34) 
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QfMq%‘  -  Wo2*//] 


0O2  =  |n0|2,  (4.59) 

so  that  we  have  a  smallest  eigenvalue  |302  equal  to  the  receiver  noise  channel  power.  This 
smallest  eigenvalue  then  defines  a  minimum  servo  gain  factor  p0,  from  Eq.  (4.45): 

m0  =  yfio2-  (4.60) 

Since  quiescent  steady-state  weight  W(°°)  must  be  equal  to  W„  by  definition,  we  can  apply 
Eqs.  (4.49),  (4.46),  and  (4.12)  to  get 

w'*  ’  (rr^K* '  (rr^K 


or 


h  =  (1  +  M0). 


(4.61) 


Thus,  the  coefficients  of  the  input  beam-steering  signal  B*  are  always  greater  than  tho 
coefficients  of  Vfq  by  the  factor  (1  +  p0). 

Combining  Eqs.  (4.49),  (4.46),  (4.41),  (4.12),  and  (4.61),  one  can  rewrite  Wk(°°)  in 
terms  of  Wq  to  obtain 


(4.62) 


where 

Kk  =  («**%).  (4.63) 

Wq  k  is  the  fcth  component  of  quiescent  beam-steering  weight  vector  V/q  in  orthonormal 
space. 


If  we  assume  quiescent  conditions  up  to  time  t  =  0,  with  only  receiver  noise  present, 
so  that  the  external  interference  sources  are  switched  on  at  t  =  0,  then 


wk(  0)  =  Wqk,  (4.64) 


and  Eq.  (4.48)  may  be  manipulated  into  the  convenient  form 
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This  equation  is  convenient  because  (nk  -  p0)  will  be  zero  for  all  eigenvalues  that  are 
equal  to  0O2,  and  those  terms  are  thus  eliminated  from  the  computation.  The  expression 
will  be  used  in  following  sections. 


4.2/  Signal-to-Noise  Optimization 

It  is  well  known  that  a  uniformly  weighted  array  gives  the  maximum  signal-to-noise 
ratio  when  the  noise  contributions  from  the  element  channels  have  equal  power  and  are 
uncorrelated.  These  conditions  are  approximately  valid  when  receiver  noise  and  uniformly 
distributed  sky  noise  are  the  predominant  noise  contributions;  they  pertain  exactly  in 
linear  array  antennas  with  half-wave  spacing.  However,  when  there  is  directional  interfer¬ 
ence  from  other  in-band  transmitters,  from  jammers,  or  from  natural  phenomena,  the 
noise  out  of  the  element  channels  will  be  cor  related,  and  uniform  weighting  will  not 
optimize  signal-to-noise  ratio.  The  solution  to  the  general  problem  is  readily  obtained 
by  the  elegant  mathematical  approach  of  maximizing  ratios  of  quadratic  forms.  Refer¬ 
ences  24,  25,  and  26  are  recommended  for  a  discussion  of  the  optimization  procedure 
and  a  more  detailed  mathematical  treatment.  Other  optimization  procedures  and  alternate 
performance  measures  may  also  be  found  described  in  the  literature  [14,21,27-39]. 

Maximum  signal-to-noise  (S/N)  ratio  per  se  is  not  really  the  desired  object  of  our 
optimization,  because  we  are  willing  to  compromise  on  S/N  ratio  to  buy  some  control 
over  the  quiescent  steered-beam  characteristics.  For  example,  it  may  be  desirable  to  con¬ 
trol  the  main-beam  shape,  sidelobe  levels,  pattern  null  placements,  or  array  phase  center. 
Since  these  desired  constraints,  or  controls,  must  be  incorporated  in  the  input  beam¬ 
steering  vector  B*,  it  follows  that  one  should  optimize  on  equivalent  signal  vector  B. 

Thus,  let  us  assume  that  the  array  output  signal  power  desired  is  given  by  the  equivalent 
expression 


s  =  |W'B|2. 


(4.66) 


The  array  output  noise  power  is  assumed  to  derive  from  the  quiescent  receiver  channel 
noise  plus  the  noise  signals  received  from  external  sources  of  interference,  as  defined  in 
Eq.  (4.2)  for  the  element  channel  components  of  E.  Output  noise  power  is  therefore 
given  by 


n  =  |W'E|2, 


(4.67) 


and  we  can  formulate  our  signal-to-noise  performance  index  as  a  ratio  of  these  two  quad¬ 
ratic  forms: 


|W;B12 

IW'EI2 


(4.68) 


Equation  (4.68)  can  be  manipulated  readily  into  a  ratio  of  Hermitian  matrix  forms,  as, 

£  =  [W'B]*[B'W]  =  W*'[B*B']W 
”  [W^E]  *[E^W]  W*f  [E^E7  ]  W 
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s_  =  W*'AW 
n  W*'MW  ’ 


(4.69) 


where  A  is  a  positive,  semidefinite,  Hermitian  matrix  (a  one-term  dyad)  and  M  is  the 
positive,  definite,  Hermitian  covariance  matrix  described  in  Eqs.  (4.23)  through  (4.28). 

In  accordance  with  the  procedure  outlined  in  Ref.  26,  the  optimization  of  Eq.  (4.69)  re¬ 
sults  in  the  eigenvalue  expression 


AW  =  -  MW, 
n 


(4.70) 


where  s/n  now  represents  an  eigenvalue.  The  maximum  value  of  s/n  which  shall  be  de¬ 
noted  with  a  zero  subscript,  as  ( s/n)0 ,  is  the  largest  eigenvalue  of  Eq.  (4.70).  Further¬ 
more,  because  of  the  properties  of  the  matrixes  involved,  it  also  happens  to  be  the  only 
nonzero  eigenvalue.  The  unique  eigenvector  W0  associated  with  eigenvalue  (s/n)0  there¬ 
fore  represents  the  optimum  element  weights.  Thus  we  have 


'•  ■  (4 


(4.71) 


Substituting  for  (s/n)0  from  Eq.  (4.69)  and  cancelling  the  common  term  (BfW0),  one 
obtains 


'W*'B*  > 


MW0. 


(4.72) 


The  quotient  on  the  right-hand  side  is  just  a  complex  number,  which  we  shall  denote  by 
C.  The  desired  optimum  weight  vector  is  then  obtained  by  inversion  of  Eq.  (4.72),  or 


W0  =  4-(M‘1B*J. 


(4.73) 


Hence,  the  optimum  weights  may  be  'obtained  directly  by  inverting  the  known  (or  esti¬ 
mated)  covariance  matrix. 

It  is  interesting  to  examine  Eq.  (4.72)  under  the  assumption  of  quiescent  conditions, 
whereby  the  covariance  matrix  would  be  ,  as  in  Eq.  (4.27),  and  the  “optimum” 
weight  mu3t  be  W?,  from  Eq.  (4.4),  by  definition.  The  quotient  term  reduces  to  unity 
for 


B*  =  MgW?  , 


(4.74) 


and  Eq.  (4.72)  is  then  identical  to  Eq.  (4.74).  This  establishes  the  optimum  relationship 
for  B*  when  W„  has  been  chosen  to  satisfy  de.ired  quiescent  beam  characteristics.  For 
the  simple  quiescent  noise  conditions  assumed  in  this  discussion,  is  a  diagonal  matrix 
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and  results  in  B*  components  being  equal  to  the  Vf„  components  multiplied  by  constants. 
Recall  that  these  constants  were  denoted  as  bk  in  Eq.  (4.12)  and  that  wo  evaluated  the 
bk  in  Eq.  (4.61). 

If  Eq.  (4.74)  is  substituted  into  Eq.  (4.72)  for  the  general  case,  our  control  law  for 
optimum  weights  becomes 


where 


B*  =  MqV/q  =  CMWq 

'w$'m9w91 


C  = 


WofMW0 


(4.75) 


(4.76) 


This  optimum  control  law  can  be  converted  into  orthonormal  system  variables  by  multi¬ 
plying  through  from  the  left  with  Q*  and  substituting  for  W0  from  Eq.  (4.52), 
whereupon 

Q*B*  =  CQ*MQfW0 
B*  =  C[ft.2S..]W0 


or 


Wn  «  4r\ 


SU 


ft2 


B* 


(4.77) 


Equation  (4.77)  is  the  basis  for  the  Wok  defined  in  Eq.  (4.46),  since  the  feth  component 
of  the  above  W0  vector  may  be  written  as 


WO 


ok 


•m 


B* 


(4.78) 


where  nk  =  70*2.  The  constant  7/C  may  be  ignored,  since  W0  may  be  multiplied  by  any 
nonzero  constant  without  changing  the  value  of  (s/n)0.  Thus,  we  have  derived  the  neces¬ 
sary  relationships  for  calculating  the  optimum  weights  from  a  known  (or  estimated)  co- 
variance  matrix. 


4.3.  Q-Transformation  Physical  Network  Analogy 

The  Q-matrix  transformation  as  defined  in  Eq.  (4.33)  is  a  rather  complicated  mathe¬ 
matical  matrix  operator  composed  of  normalized  and  mutually  orthogonal  eigenvectors. 
References  22  and  23  discuss  the  related  eigenvalue  problem  and  the  procedures  for 
.  ring  these  eigenvectors  from  known  covariance  matrix  M.  Fortunately,  computer 
p.  grams  are  available  fo*-  carrying  out  the  laborious  calculations  involved.  Despite  the 
complexity  associated  with  their  mathematical  evaluation,, however,  Q-matrix  eigenvectors 
have  a  relatively  simple  interpretation  in  term-.  oi.  physical  feed  networks,  and  this  inter¬ 
pretation  will  be  presented  herp  ^  ^ive  a  better  insight  into  the  operation  of  adaptive 
arrays. 
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Essentially,  the  components  of  the  eigenvectors  may  be  interpreted  as  array  element 
weights,  giving  rise  to  a  set  of  orthogonal,  normalized  eigenvector  beams.  The  tth  eigen¬ 
vector  beam  may  be  expressed  in  the  form 


K 

gt{B)  -  (S'e,)  =  £  eikSk 


k= l 


(4.79) 


where  S  and  its  components  Sk  were  defined  in  Eqs.  (4.8)  and  (4.9),  respectively.  Next  I 

define  a  variable  Z ,  related  to  spatial  angle  6  as  ! 

Z  =  e;2u  (4.80)  j 

where  u  =  (jidfk)  sin  0,  as  defined  in  Eq.  (4.10).  The  locus  of  Z  is  the  unit  circle  in  the  j 

complex  Z  plane.  If  we  factor  out  the  term  Sx  from  the  summation  of  Eq.  (4.79),  g,(0 )  { 

may  be  rewritten  in  the  variable  Z  as  j 


/  1  \K'1 

gi(6)  =  -  [en  +  e/2Z  +  e,.3Z2  +  •••  +  e^Z*’1],  (4.81) 

which  expresses  the  eigenvector  beam  in  familiar  array  polynomial  form  wherein  the 
eigenvector  components  become  the  coefficients  of  the  array  polynomial. 

From  the  work  of  Schelkunoff  f40],  we  know  that  an  array  space  factor  F(Z)  has 
the  two  related  forms 

F(Z)  =  aQ  +  a^Z  +  a2Zz  +  a3Z3  +  •  •  •  +  (4.82) 

or 

F(Z)  =  <Ik.\(Z  -  Z])(Z  -  Z2)(Z  -  Z3)  •  (Z  -  ^jf_i);  (4.83) 

i.e.,  it  may  be  expressed  either  in  the  polynomial  form  or  as  the  product  of  zero  factors 
containing  the  roots  of  the  polynomial.  The  roots  Zlt  Z2,  Z3,  . . . ,  ZK.1  are  the  zeros 
or  null  points  of  the  array  space-factor  pattern.  Knowing  the  null  points,  one  can  solve 
for  the  array  polynomial  coefficients,  or  conversely,  knowing  the  coefficients  one  can 
find  the  null  points.  All  of  the  (Jf-1)  null  points  will  be  located  on  the  Z-plane  unit 
circle. 


Schelkunoff’s  null-point  concept  is  particularly  applicable  to  the  patterns  associated 
with  adaptive  arrays  for  the  following  reasons: 


1.  Adaptive  arrays  form  nulls  in  the  directions  of  interference  sources. 

2.  An  array  oi  K  elements  possesses  (K-l)  degrees  of  freedom,  which  are  repre¬ 
sented  by  the  (K- 1)  nulls.  The  behavior  of  adaptive  arrays  relates  to  how  many  degrees 
of  freedom  (or  nulls)  are  “captured”  by  the  interference  environment. 
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A,  S,  A,  2 


A,  2  A  A, 


(a)  Constraint  of  one  null 


wW 


(b)  Constraint  of  two  nulls 
Fig.  27— Davies  null-control  networks 


3.  Controlled  null  placement  may  be  incorporated  into  the  input  beam-steering 
array  weights. 

4.  Constrained,  null  positions  are  usually  associated  with  the  eigenvector  beams. 

Networks  in  which  independent  control  of  pattern  nulls  is  leadily  achieved  have 
been  described  in  the  literature  [41  ],  and  it  is  instructive  to  review  their  operation  for  a 
simple  six-element  array.  Consider  first  the  network  shown  in  Fig.  27a,  wherein  one  de¬ 
gree  of  freedom  is  consumed  to  accommodate  the  constraint  of  placing  one  null.  We 
form  an  array  of  (jK-1)  “elements”  in  which  the  new  “elements”  are  subapertures  con¬ 
sisting  of  adjacent  element  pairs  phased  by  <l>i  to  the  direction  of  an  interference  source. 
Each  subaperture  has  a  Ej  (sum)  and  Al  (difference)  output  port.  All  of  the  ports 
contain  power  from  the  source,  but  the  Aj  ports  would  receive  little  or  no  power  because 
of  the  null  being  directed  toward  the  source. 
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The  network  shown  in  Fig.  27b  consumes  two  degrees  of  freedom  to  accommodate 
the  constraint  of  placing, two  nulls.  We  now  form  an  array  of  (K-2)  “elements”  in  which 
the  new  “elements”  are  subapertures  consisting  of  three  adjacent,  elements  phased  by  0j 
and  02  to  the  directions  of  two  interference  sources.  Each  subaperture  has  three  outputs 
of  interest,  labeled  SiA2,  22A1,  and  AjA2.  The  2jA2  ports  receive  power  from  the 
source  at  6X  but  none  from  the  source  at  02,  and  vice-versa  for  the  22A1  ports.  The 
unique  AXA2  ports  receive  little  or  no.  power  from  either  source  because  of  the  two 
directed  nulls. 

This  network  implementation  of  controlled  nulls  can  be  continued  until  ali  (Jf-3) 
degrees  of  freedom  have  been  consumed,  whereupon  one  obtains  a  complete  “Davies 
tree”  matrix  network,  as  shown  in  Fig.  28.  The  network  has  six  output  ports  of  interest 
to  our  discussion,  and  we  can  write  the  associated  array  pattern  functions  by  inspection 
for  half-wavelength  element  spacing,  because  2y  and  Ay  are  simple  cosine  and  sine 
functions: 

*l(0)  (^i^2^3^6^5) 

*2(0)  k 

*3(0)  *  A223A4A5) 

*4(0)  ^  (Ai2223A3A2) 

*g(0)  88  (^iA^-jAjjAj) 

*6(0)  (^1^2^3^4^s) 

where 

2y  =  COS  (U  -  My) 

Ay  =  Sin  (U  -  My). 

We  have  K  beams  which  may  be  used  to  resolve  up  to  (ff-1)  interference  point  sources. 
Note  that  the  last  beam  g'6(d)  is  unique  in  that  it  alone  contains  all  CK~JL)  controlled  nulls. 
This  network  represents  an  intuitive  first  approximation  to  the  operation  of  the  Q-matrix 
transformation  because  it  is  capable  of  achieving  a  shaped  spatial  coverage  incorporating 
constrained  null  suppression  of  interference  sources  up  to  the  limits  of  its  degrees  of 
freedom.  However,  except  for  those  special  cases  involving  orthogonal  sets  of  uniform 
illumination  beams,  the  network  cannot  serve  as  an  exact  analogy  because  all  of  the 
beams  represented  in  Eq.  (4.84)  are  tied  together  in  a  fixed  relationship  by  the  same 
(JC-l)  null  points,  i.e.,  by  the  same  phase  shifters.  This  restrictive  relationship  prevents 
the  beam  array  vectors  from  possessing  the  optimizing  characteristics  associated  with 
true  eigenvectors. 

To  achieve  an  exact  representation  for  the  Q-matrix  transformation,  each  eigenvector 
:..ust  have  its  own  (Jf-1)  null-controlled  network,  i.e.,  a  network  similar  to  the  one  as¬ 
sociated  with  beam  g'6(d)  in  Eq.  (4.84).  In  addition,  a  multiplicative  constant  must  be 
added  to  permit  normalization.  Figure  29  illustrates  a  simplified  schematic  of  such  a 
netv,  ork  for  the  *th  eigenvector,  where  the  orthonormal  output  voltage  is  Et  and  the 
eigenvector  beam  may  be  written  as 


(4.84) 


(4-85) 
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Fig.  29— Simplified  schematic  of  (JC  -  1)  controlled-null  network  representing  one 
eigenvector  of  the  Q-matrix  transformation 


gi(0)  -  C(A1A2A3A4A5  .  (4.86) 

The  output  power  from  this  network  will  be  exactly  equal  to  the  eigenvalue  &f“,  recall 
from  Eq.  (4.39)  that 


{E*E,)  =  tf. 


(4.87) 


Each  null  in  Eq.  (4.86)  corresponds  to  a  zero  of  the  associated  array  space  factor  poly¬ 
nomial  located  on  the  E-plane  unit  circle.  If  we  denote  Zr  as  the  rth  zero,  then  from 
Eq.  (4.80), 

(Z  -  Zr)  =  e’2u  -  e'ZUr 


-  j2e'^u+u^  sin  (u  -  ur), 


(4.88) 


and  we  can  readily  convert  Eq.  (4.86)  into  the  form  of  products  of  zero  factors: 


VI 


(6)  = 


_C£i_ 
(2 sJZf 


[(Z  -  ZX)(Z  -  Z2)(Z  -  Z3hZ  -  Z4)(Z  -  Z6)] 


(4.89) 
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where 


r 


-f(u1+u2+u3+u4+u6) 


(4.90) 


To  demonstrate  an  application  of  the  eigenvector  network  of  Fig.  29,  let  us  take  the 
case  of  a  single  narrowband  source  located  at  angle  d1.  Eq.  (4.31)  shows  that  this  single 
source  would  result  in  only  one  unique  eigenvalue  /3-j2;  the  other  eigenvalues  would  be 
multiple  roots  assumed  equal  to  |302,  the  output  receiver  noise  power  (Eq.  (4.59)).  There 
would  be  only  one  unique  eigenvector  ei  associated  with  0J2;  the  other  (K- 1)  eigenvec¬ 
tors  are  not  unique.  Eigenvector  ei  is  found  to  be  equal  to 

-^[e/6\e,3\e'\e-'\e-'S\e-'6">]  (4.91) 


where 


7r  .  a 

«1  =  2  sm  01- 


(4.92) 


Note  that  the  eigenvector  phasing  is  simply  the  complex  conjugate  of  the  element  signals 
received  from  the  single  source  interference.  From  Eqs.  (4.81)  and  (4.83),  one  obtains, 
for  unique  eigenvector  beam  g'1(0). 


*l(0)  =  (S'e,)  =  — ==rr  F(Z') 

•JWf 

F(Z')  =  1  +  Z'  +  (Z') 2  +  ( Z ')3  +  ( Z ')4  +  (Z')6 


(4.93) 


(4.94) 


F(Z')  =  (Z'  -  Z0)(Z'  -  Z02)(Z'  -  Z03)(Z'  -  Z04)(Z'  -  Z05) 


(4.95) 


where 


Z'  =  (^)=  e/2(“'Ul)  and  Z0  =  (“)  =  e> (4.96) 

Equation  (4.94)  is  recognized  as  a  uniform-illumination  array  factor,  and  one  can  readily 
convert  g1  (0)  to  the  familiar  'r  gonometric  form 


fl,  sy  6  (<■  -  Hi) 

y/6  sm  (u  -  Uj) 


(4.97) 


Thus,  the  unique  eigenvector  beam  is  a  uniform-illumination  beam  centered  on  the  source 
at  flj.  Figure  30  shows  the  Z-plane  unit  circle  with  source  point  S1  plus  all  five  null 
points  equally  spaced  around  the  circle  by  2n/K  rad.  If  the  network  phase  shifters  in 
Fig.  29  are  set  to  produce  the  null  points  shown  in  Fig.  30,  the  correct  eigenvector  beam 
will  be  obtained. 

The  formation  of  the  remaining  nonunique  eigenvectors  can  be  accomplished  easily 
for  this  single-source  ease  by  using  the  null  points  of  Fig.  30  to  successively  center  uniform 
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Fig.  30—2-plane  unit  circle  null  locations  for  single¬ 
source  case 


illumination  beams  at  N1,  N2,  iV3,  N4,  and  N5.  Note  that  each  of  these  beams  incor¬ 
porates  constrained  null  point  Sx,  which  is  the  location  of  the  jammer.  This  will  result 
in  a  set  of  orthogonal,  uniform  illumination  beams  such  as  that  obtained  from  a  Butler 
matrix  feed  network.  The  associated  eigenvectors  are  similar  to  ex  and  result  in  the  fol¬ 
lowing  complete  Q  matrix: 


Q  = 


V6 


'ei5u' 

tH 

3 

CO 

’•N 

*> 

e/“i 

g-/«x 

C-/3"l 

e'/5“l' 

e/5“* 

e/3“2 

e/“2 

e-/«2 

e‘/3“2 

c‘/B“2 

e/6u3 

e/3“3 

e/«s 

e/5“4 

„/5«  5 

• 

• 

• 

• 

• 

(4.98) 

*'6“6 

If  we  use  element  signal  components  as  defined  for  Eq.  (4.2)  with  the  single  source  lo¬ 
cated  at  0j,  i.e.,  if 


r  -  n  4.  t  Ju i<2*-7> 
“k  ~  nk  +  "le  » 


(4.99) 


then  the  transformed  orthonormal  output  signals  are 
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E  =  QE  = 


(»01  +  >/®"  *^l ) 
n02 
n03 
«04 
n05 
n06 


(4.100) 


where  nQj  represents  the  summation  of  quiescent  receiver  noise  voltages  nk  at  the  ith  out¬ 
put  port.  Define  an  input  noise  vector  N  consisting  of  the  nk, 


N'  =  [nlt  n2,  n3,  n4,  n5,  n6]  . 


(4.101) 


Then 


N  =  QN  or  n0i  =  (e/N). 

Cross-correlating  any  two  of  the  output  noise  voltages,  results  in 


(4.102) 


(*0%/>  =  (S*'N*)(N  »6/) 
=  ef'  [N*N,]e/ 


(4.103) 


If  we  assume  equal  magnitudes  of  receiver  noise  power  in  each  channel,  such  that  |n^|2  = 
|n0|2  for  all  ft,  then  Eq.  (4.103)  may  be  written  as 


(»o/"o/)  =  l"ol2(e?'e/) 


(4.104 


Thus,  although  each  .output-port  noise  consists  of  a  mixture  of  all  the  input  noise  voltages, 
the  outputs  have  zero  cross-correlation  because  of  the  orthogonality  of  the  Q-matrix 
eigenvectors.  Also,  the  noise  power  at  any  output  port  equals  the  input  receiver  noise 
power,  since  eigenvectors  are  normalized  to  unity  Hermitian  length: 


l»o/i2  =  («w«o«)  =  K!2- 


(4.105) 


The  squares  of  the  absolute  ••alues  of  the  orthonormal  output  signals  in  Eq.  (4.100) 
equal  the  eigenvalues  of  the  covariance  matrix, 


0i2  =  (ln0l2  +  6U1I2) 

0/2  ~  ln0l2  =  0o2  for  *  =£  1» 


(4.106 
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where  we  note  that  the  power  from  the  interference  source  appears  at  only  one  output 
port  for  this  single-source  case. 


It  is  also  of  interest  to  see  what  is  obtained  for  transformed  input  beam-steering 
vector  B*,  from  Eqs.  (4.98)  and  (4.41).  Since  the  eigenvectors  are  all  uniform- 
illumination  array  factors,  we  have,  from  Eqs.  (4.12),  (4.61),  and  (4.5), 


B*  =  (e*'B*) 


=  (l+MoXef'W  ) 


1  +^o  £ 

La*e 


/(«r"0)(  2*-7) 


(4.107) 


From  Eq.  (4.7),  we  see  that  the  summation  equals  the  value  of  the  quiescent  beam  pattern 
at  angle  dh  or  G_  (0,  ),  so  that  B*is  simply 


B*  =  QB* 


Gq(h) 

_  /I  +  Up\  G?(02> 

=  U/6  /  G?(03) 


(4.108) 


G?M 


A  . 

So  we  find  that  the  components  of  B  are  proportinal  to  the  quiescent  beam  pattern 
sampled  at  the  orthogonal  eigenvector  beam  positions. 


We  have  evaluated  the  complete  Q  transformation  for  the  case  of  a  single  source  of 
interference  and  have  interpreted  the  eigenvalues  and  eigenvectors  in  terms  of  the  physical 
network  analogy  of  Fig.  29.  Further  application  of  this  concept  will  be  made  in  follow¬ 
ing  sections  for  more  complicated  distributions  of  interference  sources. 


4.4.  Retrodirective  Eigenvector  Beam  Concept 


A  valuable  insight  into  a  fundamental  principle  of  operation  for  adaptive  arrays  may 
be  gained  by  examining  the  formation  of  retrodirective  beams,  as  illustrated  in  Fig.  1. 

We  saw  in  the  previous  section  that  for  a  single  narrowband  source  of  interference,  we 
obtain  one  unique  eigenvalue  and  one  unique  eigenvector  which  produces  a  uniform- 
illumination,  retrodirective  eigenvector  beam  centered  on  the  source  at  9X,  as  given  in 
Eq.  (4.97).  Note  that  even  though  a  complete  set  of  K  orthogonal  uniform-illumination 
beams  was  set  up  by  the  Q-matrix  transformation  network,  only  one  of  those  beams  was 
retrodirective  toward  the  jammer,  and  it  was  produced  by  the  one  unique  eigenvector.  It 
will  be  shown  in  this  section  that  adaptive  array  pattern  performance  can  be  character¬ 
ized  by  considering  only  the  unique,  retrodirective  eigenvector  beams;  the  arbitrary, 
nonunique  eigenvector  beams  are  not  essential  and  need  not  be  evaluated. 
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The  necessary  array  output  pattern  function  G(d,t)  can  be  derived  most  readily  by 
considering  the  output  of  the  orthonormal  system  of  Fig.  25b  for  array  input  signal 
vector  S,  defined  in  Eq.  (4.8)  and  recalling  that  the  output  must  be  identical  for  both 
the  real  system  and  the  orthonormal  system;  i.e., 


K 

g(0,o  -  £  w&  -  £  m 

/=!  i=l 


(4.109) 


or 


G(9,  t)  =  (W'S). 


(4.110) 


The  vector  S  of  course  results  from  the  Q-matrix  transformation  operating  on  the  input 
signal  vector  S: 


S  =  QS. 

From  Eq.  (4.33)  we  see  that  ith  component  S{  would  be  given  by 


(4.111) 


K 


S,  -  (e/S)  -  £  «»«*• 
Jt=l 


(4.112) 


But  this  summation  defines  the  ith  eigenvector  beam  as  in  Eq.  (4.79),  so  that 

^  =  (e/S)  =  gi{6) 
or 


(4.113) 


K 

G(9,t)  -  £  Wtgf[9). 

i=l 


Thus,  our  output  pattern  function  is  a  summation  of  the  K  eigenvector  beams  weighted 
by  the  orthonormal  system  adaptive  weights. 


A  convenient  equation  for  Wj  has  already  been  developed  in  Sec.  4.1  under  the  as¬ 
sumption  that  quiescent  noise  conditions  hold  up  to  time  t  =  0,  when  the  external  sources 
of  interference  are  suddenly  switched  on.  Repeating  Eqs.  (4.65)  and  (4.63),  it  was  found 
that  Wt  may  be  written  as 


w,  =  Wqi  -  (1  -  e'01**) 


(4.114) 
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*  =  1 

"H  is  the  Hermitian  scalar  product  of  the  quiescent  beam-steering  vector  Vtq  and  the  ith 
eigenvector.  It  represents  the  ith  component  of  W?  in  transformed  orthonormal  space. 
Note  that  at  time  t  =  0,  Eq.  (4.114)  reduces  to  Ws  =  Wqj,  and  Eq.  (4.113)  results  in 

K 

G(6, 0)  =  £  Wqtgt(6)  =  (fy'S)  =  (W/QS).  (4.116) 

i=l 

But  from  Eq.  (4.52)  we  see  that 

W,'  =  Wq'Q  (4.117) 

or 

0(6,0)  =  (W/S)  =  Gq(6),  (4.118) 

where  quiescent  beam  pattern  Gq(9)  was  defined  in  Eq.  (4.7).  This  result  could  be  an¬ 
ticipated,  of  course,  from  the  sampling  properties  of  the  orthonormal  eigenvectors  or 
merely  from  the  fact  that  the  quiescent  outputs  from  the  two  systems  must  be  identical. 

The  final  step,  then,  is  to  substitute  Eqs.  (4.114)  and  (4.118)  into  Eq.  (4.113), 
whereupon  we  obtain  the  desired  relationship, 

0(9, t)  =  Gq(9)  -  g  (1  -  «■“'')(— 7y)^|(0).  (4.119) 

Recall  that 

1  +  Hi 


=  yPo2- 

Thus,  the  output  pattern  function  of  our  adaptively  controlled  linear  array  consists  of 
two  parts.  The  first  part  is  quiescent  beam  pattern  Gq(6),  and  the  second  part,  which 


I 


and 


Mi 
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is  subtracted  from  Gq(6),  is  a  summation  of  weighted,  orthogonal,  eigenvector  beams. 

An  important  point  in  the  weighting  is  that  the  numerator  (p,-  -  p0)  will  be  zero  tor  all 
eigenvalues  that  are  equal  to  quiescent  eigenvalue  0O2. 

Therefore,  one  may  disregard  all  0q2  associated  eigenvector  beams  in  the  summation, 
retaining  only  the  unique  eigenvector  beams,  which  also  happen  to  be  retrodirective. 
Another  important  point  in  the  weighting  is  that  the  transient  response  time  of  the 
unique  eigenvector  beams  is  controlled  by  ctj,  which  is  proportional  to  the  eigenvalue. 

A  large  eigenvalue  implies  a  fast  transient  response  for  its  associated  eigenvector  beam, 
whereas  a  small  eigenvalue  results  in  slow  response. 

Let  us  first  apply  the  above  pattern  function  to  the  case  of  the  single  narrowband 
source  located  at  angle  0lf ; which  was  discussed  in  the  previous  section.  Since  there  is 
only  one  unique  eigenvalue  0i2,  we  have  only  one  nonzero  term  in  the  summation,  so 
that  Eq.  (4.119)  reduces  to 

g(M)  -  oq(d)  -  (i  -e'v)(~rirr)  ^i*i(0)>  <4-120> 

where  =  (1  +  px/t0)  and  -  70x2.  From  Eq.  (4.97)  in  the  previous  section, 

(4.121) 


1  [sin  K(u  -  u-i )' 

*<*>  = 


and  from  Eqs.  (4.107)  and  (4.7), 


*,l  = 


-Mil 


(4.122) 


From  Eq.  (4.106),  0X2  was  evaluated;  it  can  be  used  to  convert  the  servo  gain  factor  term 
to  a  more  meaningful  form: 


01  =  (fh0r  +  K\h\l) 

/v  |2  j.  tr\~r  |2 

Mo  =  “ 


Ml  _  2HwqI  +  f  !£ii )  -  n  +  KP, 

7IS0I2  ^ 


(4.123) 


where  is  now  the  ratio  of  jammer  power  to  receiver  noise  power  at  the  preamplifiers. 
If  Hq  is  set  equal  to  unity,  which  would  be  normal  practice  for  the  circuit  of  Fig.  24, 
the  servo  gain  term  becomes 


/Ml  -  M0\  , 

f  \ 

Ui+irl 

(2  +  JCPj 

(4.124) 


and  Eq.  (3.120)  can  be  rewritten  as 
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Fig.  31— Steady-state  adapted  pattern  for  single  source.  Case  A;  power  ratio  Py  =  1,250, 

located  at  6y  =  21° 


„  //lx  -attJPiGeViAfdnKiu-uJ 

0(e,t)  =  0,(6)  -  (1  - .  )(-2,.ip1  »iT. 


(4.125) 


This  gives  us  the  performance  for  the  single-jammer  case  in  easily  understood  variables, 
without  any  need  for  either  eigenvalues  or  eigenvectors. 


To  illustrate  the  use  of  Eq.  (4.125),  Fig.  31  is  a  plot  of  the  quiescent  Gq(9)  pattern 
(dotted  line)  and  the  steady-state  adapted  pattern  (solid  line)  computed  for  an  eight- 
element  linear  array  (if =8)  with  a  jammer  of  power  ratio  Py  =  1250  in  the  first  sidelobe, 
at  dy  =  21°.  Gq(6)  was  chosen  as  a  uniform-illumination  pattern  steered  to  broadside. 
Figure  32  illustrates  the  two  parts  of  the  adapted  pattern  on  an  expanded  angular  scale: 
the  quiescent  Gq($)  pattern  (dotted  line)  and  the  retrodirective  eigenvector  beam  pattern 
g[(0)  (solid  line)  which  has  been  multiplied  by  its  weighting  factors,  i.e.. 


/  Pi  \  /sin 


/sin  K(u  -  ut) 
\  sin  (u  -  iiy) 


(4.126) 


Note  that  the  weighting  factors  cause  the.  ^euk  of  the  retrodirective  eigenvector  beam  to 
become  alined  exactly  with  the  sidelobe  magnitude  corresponding  to  the  angular  position 
of  the  jammer.  This  produces  the  deep  null  in  the  adapted  pattern  when  the  two  parts 
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are  subtracted,  of  course,  and  one  can  readily  compute  the  depth  of  that  null  at  the 
jammer  position  from  Eq.  (4.125): 


'sin  K  (u  -  u1)\  _ 

^(u-Uj)  / 


for  9  =0i 


or 


001,0 


1  -  (1  -  e 


0,0!). 


(4.127) 


In  the  steady  state  t  =  Eq.  (4.127)  will  reduce  even  further  to  the  simple  form 


(4.128) 


Thus,  for  strong  jammers  the  gain  is  reduced  by  the  considerable  factor  of  the  square  of 
KP±  in  the  direction  of  the  jammer,  effectively  eliminating  the  jammer  power  from  the 
output  of  the  adapted  array. 
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Fig.  33— Transient  response  for  single  source,  Case  A,  for  source  power  ratios  of  21  dB  and 
31  dB;  eigenvalues  1,001  and  10,001,  respectively 
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Transient  response  can  be  shown  conveniently  by  plotting  tne  increase  in  output 
noise  power,  proportional  to  ^i[G(0i,f)]2  for  this  case,  vs  time.  The  values  of  G(01,  t) 
are  calculated  from  Eq.  (4.127).  Figure  33  illustrates  such  a  transient  response  (solid 
line)  where  P j  =  1,250  and  T0  =  12,750  ps,  which  result  in  0^  =  0.784.  The  dashed- 
line  plot  was  ccjrpule'l  1c.-  10  dB  less  jammer  power  Pj  =  125,  whereby  a*  =  0.078 
and  the  convergence  time  is  correspondingly  ten  times  longer. 

Next  apply  Eq.  (4.119)  to  the  case  of  two  narrowband  interference  sources  of 
nearly  equal  power  ratios  Pj  =  1,250  and  P2  =  1,200,  close  together  in  the  first  sidelobe 
of  the  quiescent  pattern,  at  01  =  18°  and  02  =  22°,  respectively.  If  one  forms  the  co- 
variance  matrix  (Eq.  (4.26) j  for  this  case  and  solves  for  the  unique  eigenvalues  and 
eigenvectors,  the  solutions  listed  in  Table  1  are  obtained.  The  exact  procedure  for  ar¬ 
riving  at  these  solutions  will  be  discussed  in  Sec.  4.5.  Note  that  there  are  two  unique 
eigenvalues  and  two  associated  unique  eigenvectors,  but  that  the  ratio  of  the  two  eigen¬ 
values  bears  no  relationship  to  the  jammer  powers,  which  are  essentially  equal  in  this 
case.  Such  widely  different  eigenvalue  solutions  will  be  found  to  be  typical  of  situations 
in  which  sources  are  close  together  in  terms  of  array  beamwidth.  Using  the  solutions  for 
@l2<  $%'  ei>  anc*  e2  fr°m  Table  1,  we  can  evaluate  the  associated  pf,  a,-,  Wqi,  and  g((0) 
for  substitution  into  Eq.  (4.119).  Figure  34a  illustrates  the  resulting  steady-state  adapted 
pattern  (solid  line)  for  our  eight-element  array,  and  Fig.  34b  shows  an  expanded  plot  of 
the  sidelobes  in  the  immediate  vicinity  of  the  jammers. 
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Table  1 

Eigenvalues  and  Eigenvectors  for  a  Two-Source  Case* 


Eigenvector 

Coefficient 

Amplitude 

Phase  Angle 
(deg) 

1,1 

0.340120 

69.7477 

1,2 

0.351501 

8.3626 

1,3 

0.359156 

-53.03 

1,4 

0.363004 

245.573 

1,5 

0.363004 

184.175 

1,6 

0.359156 

122.778 

1,7 

0.351501 

61.3852 

1,8 

0.340120 

0 

2,1 

0.537157 

250.308 

2,2 

0.387808 

188.701 

2,3 

0.234347 

126.989 

2,4 

0.078413 

64.5464 

2,5 

0.078413 

185.762 

2,6 

0.234347 

123.319 

2,7 

0.387808 

61.607 

2,8 

0.537157 

0 

♦Unique  eigenvalues:  18,544.4  and  1,057.58. 


Note  that  the  pattern  nulls  are  very  closely  alined  with  the  positions  of  the  two  strong 
jammers  in  this  case.  In  later  examples,  we  will  see  cases  where  the  nulls  are  not  so  well 
alined. 

Figure  35  shows  the  components  of  the  adapted  pattern  on  an  expanded  angular 
scale:  the  quiescent  Gq(Q)  pattern  (dotted  line)  and  the  two  retrodirective  eigenvector 
beam  patterns  g'x{d)  (solid  line)  and  (dashed  line): 


g[m  =  wqlgl{6) 
g'2(0)  =  vqi8z(e). 


(4.129) 

(4.130) 


Beam  gi{6)  covers  both  sources  in  the  manner  of  a  centered  sum  beam,  and  its  power 
gain  of  approximately  7.57  at  the  source  locations  leads  to  a  total  output  power  equal 
to  the  first  eigenvalue: 


=  1+Pi*i20i)+P2*12(®2) 


=  1  +  (1250  +  1200)7.57 
=  18,544. 


(4.131) 
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Fig.  35— Retrodirective  eigenvector  beams g[(0)  and  g'2(0)  for  two  sources.  Case  B 


Note  that  the  weighting  factor  Wqi  does  not  aline  the  peak  of  ^(0)  exactly  with  the 
sidelobe  magnitude  of  either  source,  although  it  comes  close  enough  that  its  subtraction 
would  cause  a  null  structure  close  to  the  final  nulls  in  that  sidelobe  region.  Eigenvector 
beam  g2(0)  spl»*s  the  sources  in  the  manner  of  a  difference  beam,  and  its  power  gain  of 
approximately  0.43  at  the  source  locations  leads  to  a  total  output  power  equal  to  the 
second  eigenvalue: 


2,  =  1+Pi*2*<#1>+*W<*2) 


=  1  +  (1,250  +  1,200)0.43 


=  1,057. 


(4.132) 


An  interesting  point  here  is  that  both  eigenvector  beams  contain  power  from  both  of  the 
sources,  so  that  one  may  be  curious  as  to  how  decorrelation  is  accomplished.  To  explain 
this,  we  start  with  the  ith  eigenvector  beam  output  from  Eq  (4.37)  and  substitute  the  E 
vector  components  from  Eq.  (4.2): 
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E{  =  (e/E) 


2  AT 


-(VN)  +  i;jrrZ:«*.^“‘Jr-l) 


(4.133) 


/■“l  *■! 


where  N  is  the  quiescent  receiver  noise  vector  defined  in  Eq.  (4.101).  The  noise  part 
is  nQi,  as  defined  in  Eq.  (4.102),  and  from  Eqs.  (4.79)  and  (4.9)  we  see  that  the  summa¬ 
tion  in  k  is  simply  the  value  of  the  eigenvector  beam  pattern  at  angle  0r,  or 


\  -  V  *  Jur(2k-K-\) 

L  eike 


(4.134) 


where 


nd  .  a 

ur  =  Y  sin  er 


Thus, 


*  n0l  +  £  Jr8t(8r)- 
r*  1 


(4.135) 


From  Eqs.  (4.104,  (4.39),  and  (4.59),  the  averaged  correlations  from  (4.135)  become 


(E*Et)  =  ft2  =  |n0!2  +  £  |7,I2*M) 


=  1  +  £  Prefer) 


(4.136) 


(E*E2)  =  0  =  \Ji\2giW  1X2(61)  + 


(4.137) 


Equation  (4.13/)  shows  that  the  zero  cross-correlation  can  be  related  to  the  products  of 
she  voltage  patterns;  i.e.,  the  product  of  [gi(0)g2(0)]  is  positive  when  9  =  6L  but  nega¬ 
tive  when  6  -  b2.  This  equation  may  be  further  generalized  to  a  summation  for  R 
jammer  sources: 


(E*E.)  =  |,.'0|2  £  Pr^MO,)  =  0.  i 


(4.138) 
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Fig.  36— Transient  response  for  two  sources,  Case  B;  eigenvalues  18,644  and  1,068 


Equation  (4.136)  is  the  basis  for  Eqs.  (4.131)  and  (4.132).  Since  all  of  the  terms  in  its 
summation  must  be  positive  or  zero,  the  nonunique  eigenvector  beams  that  would  result 
in  fy2  =  0O2  must  incorporate  nulLs  located  precisely  on  the  jammer  positions. 

Figure  36  shows  the  transient  response  for  this  two-jammer  case,  and  one  can  readily 
see  the  two  distinct  slopes  associated  with  the  two  different  eigenvalues.  The  increase  in 
output  noise  power  for  this  case  is  proportional  to 


PxG2{dr,t)  +  P2G2(02,f). 


(4.139) 


Eigenvector  beam  g\{0),  with  =  1.45,  attenuates  both  terms  rapidly  to  the  point 
where  second  eigenvector  beam  g2(0),  with  a2  =  0.083,  takes  over  and  completes  the 
attenuation  at  its  slower  rate.  An  exact  expression  for  the  output  noise  power  will  be 
developed  in  the  next  section. 


4.5,  Performance  Characteristics 

Initial  Conditions  and  Assumptions— The  performance  characteristics  are  calculated 
from  the  equations  developed  in  Secs.  4.1  and  4.4.  One  initial  condition  already  mentioned 
is  that  quiescent  receiver  noise  is  assumed  to  bf  the  only  system  noise  present  up  to  time 
t  =  0,  when  the  entire  selected  distribution  of  external  interference  noise  sources  is 
switched  on  in  a  single  step  function.  Another  assumption  already  discussed  is  that  the 
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receiver  noise  rms  voltage  in  all  element  channels  is  equal  and  denoted  by  n0,  such  that 
quiescent  eigenvalue  (302  is  the  square  of  its  rms  value,  as  in  Eq.  (4.59),  and  thereby 
defines  minimum,  quiescent  servo  gain  factor  /i0.  as  in  Eq.  (4.60).  Repeating  these  for 
convenience,  we  have 


Po2  ~  {no\2 
M0  =  yPo2  =  k2G\n0\2. 


(4.140) 

(4.141) 


For  the  circuit  of  Fig.  24,  it  is  convenient  to  choose  unity  for  the  value  of  p0,  and  we 
assume  that  amplifier  gains  G  will  be  set  accordingly. 


Once  Hq  has  been  defined,  it  is  convenient  to  express  M;  from  Eq.  (4.45)  as  a  ra.U; 
of  eigenvalues,  such  that 


M/ 

Mo 


(4.142) 


or 


This  expression  can  be  converted  to  jammer  power  ratios  by  Eq.  (4.136),  so  that 


M/  Mo 


R 


i  +  £  prg?(dr) 


r=  1 


(4.143 


This  gives  an  expression  similar  to  Eq.  (3.49). 


The  next  initial  condition  is  to  specify  quiescent  steered-beam  pattern  Gq  (6)  and 
its  associated  set  of  quiescent  weights  W?,  as  defined  in  Eqs.  (4.4)  through  (4.7).  For 
purposes  of  this  report,  Gq(6)  is  chosen  to  be  a  simple  uniform-illumination  beam 
formed  by  an  eight-element  linear  array  with  elements  spaced  X/2  apart; 


K  =  8  elements 
d  -  X/2  element  spacing 
ak  =  1  for  all  elements 


r  m  a  /sin  8  (u  -  u0)\ 
9  )  \sin  (w  -  u0)  ) 


(4.144) 
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where 


n  .  a 

u  -  g-  sm  0 


»,6  =  e-'“°<2‘-9). 


Then  the  coefficients  of  input  beam-steering  vector  B*  may  be  evaluated  from  Eqs. 
(4.61)  and  (4.12)  as 


bk  ■  (1  +  M0)  =  2 

A*  -  -  2C‘'"»(2#'9). 


(4.145) 


The  maximum  power  condition  must  be  considered  in  a  manner  similar  to  the  single¬ 
loop  case  discussed  in  Sec.  3.2,  and  a  relationship  similar  to  Eq.  (3.47)  can  be  derived 
for  the  orthonormal  servo  loops  of  Fig.  25b,  wherein  the  maximum  servo  gain  factor 
is  given  by 


Mm  =  Mol 


Kl\  _  frVo\ 
A?r  \  io  r  x* 


(4.146) 


where  0m2  represents  the  maximum  eigenvalue  to  be  handled,  or  the  maximum  power  to 
be  delivered,  at  any  of  the  orthonormal  output  ports.  Channel  bandwidth  Bc  and  basic 
filter  time  constant  r0  are  assumed  to  be  the  same  for  all  element  channel  servo  loops. 
Thus  solving  for  r0  yields 


r0 


/10  \ 

VBc) 


1  +  Mo 


po  J  } 


(4.147) 


or 


T0 


-(«.) 


1  +  Mq  +  Mo  £  Pr8m2(er) 


r=  1 


(4.148) 


Note  that  the  maximum  power,  or  maximum  eigenvalue,  will  be  much  larger  than 
jammer-to-receiver-noise  power  ratios  Pr  because  the  Pr  are  multiplied  by  t:  -•  power  gain 
of  the  retrodirective  eigenvector  beams.  For  example,  in  the  single-source  case  discussed 
in  connection  with  Fig.  32,  Pi  =  1,250,  and  from  Eq.  (4.121), 


SiHed  =  K  =  8 


or 


=  [l+Pj*!2^)]  =  10,001. 
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Similarly,  in  the  two-source  case  discussed  in  th?  previous  section,  we  saw  in  Eq.  (4.131) 
that  the  largest  eigenvalue  was  18,544,  whereas  P1  was  1,250  and  P2  was  1,200. 


For  the  cases  to  be  considered,  a  value  of  20,000  for  nm  will  cover  the  largest 
eigenvalues,  and  if  channel  bandwidth  Bc  is  kept  at  5  Me  as  in  the  single-loop  case,  then 
the  value  for  r0  from  Eq.  (4.146)  should  be  approximately  12,750  ps.  Summarizing 
these  selected  constants,  we  have 


= 


20,000 


Bc  =  5  Me 


(4.149) 


Tq  =  12,750  ms. 


Output  Noise  Power  and  S/N  Degradation— The  performance  factor  of  ultimate  in- 
terest  in  an  adaptive  array  is  the  improvement  in  output  signal-to-noise  ratio  as  com¬ 
pared  t< .  a  conventional  array  subject  to  the  same  interference  conditions.  In  this  ratio, 
the  output  noise  power  is  fundamental  to  the  improvement  obtained  and  usually  is  suffi¬ 
cient  by  itself  for  illustrating  the  transient  behavior  of  the  system.  To  calculate  it,  we 
take  advantage  of  the  fact  that  the  receiver  noise  is  statistically  independent  of  the  ex¬ 
ternal  interference  noise  sources,  so  that  we  can  add  their  separate  output  powers 
linearly. 


Starting  with  receiver  noise,  the  output  contribution  can  be  expressed  in  terms  of 
either  the  real  system  or  the  orthonormal  system  in  Fig.  25: 


K  K  A 

\Y0n{t)\2  -  £  l"W2  -  £  IV/»0/l2 

fc=l  1=1 


(4.150) 


where  n0f  was  defined  in  Eq.  (4.102).  Substituting  for  W(  from  Eq.  (4.65)  and  for  n0i 
from  Eq.  (4.105)  results  in 


\Y0n(t)\2  =  ln0!2£ 

i=  1 


iv2 


(4.151) 


or 


K 

ir0*(t)l2  “  l«ol2£[i-A(012IV2 

i=i 


(4.152) 


where 


At(t)  =  (1 


_e-a,t  (tLZlo 

V  +  1 


(4.153) 
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Aj(t)  will  be  zero  for  1  =  0  and  for  M/  =  Mo 5  i-e.»  for  the  nonunique  eigenvalues.  For 
quiescent  conditions  at  t  -  0,  note  that 


ir0„(0)l2  =  l«0l2i;  l</l2  =  l«ol2  E 

i=i  *=i 


K 


(4.154) 


since  the  noise  power  output  must  be  the  same  for  either  system.  Thus,  Eq.  (4.152) 
may  be  rewritten  as 


K  K 

ir0„(f)!2  -  l«ol2  E  l^l2-E  [2-Af(01A(f)i^il2  •  (4.155) 

*=1  1=1 


This  form  is  convenient  because  the  Wqi  associated  with  nonunique  eigenvectors  need  not 
be  evaluated,  since  A((t)  -  0  lot  them. 


The  noise  power  contributed  by  R  external  interference  sources  is  simply  a  summa¬ 
tion  of  their  output  power  pattern  levels, 


R 

\Y0i(t)\2  =  |n0l2  E^OW), 

r= 1 


(4.156) 


where  Pr  is  the  power  ratio  of  the  rth  source,  dr  is  its  location,  and  G{0n  t)  is  given  by 
Eq.  (4.119). 


Total  output  noise  power  is  the  summation  of  Equ.  v4.155)  and  (4.156),  and  the 
increase  in  output  noise  power  is  the  ratio  of  these  two  to  quiescent  noise  in  Eq.  (4.154): 


W)!2 


JYonm* 


"  1+^ 


R 

E  PrGHen  t)  -  E  [2  -  Ai(t)Wt) i  wqt I5 

r=l 


K 

I 
1  =  1 


l 


IX*  i! 


/t=i 


(4.157) 


This  increase  in  output  noise  power  is  the  quantity  usually  plotted  for  illustrating  the 
transient  behavior  of  the  system.  Figures  33  and  36  are  examples  of  its  application.  A 
pertinent  characteristic  of  this  performance  index  is  that  it  indicates  the  general  magni¬ 
tude  of  the  adapted  weights  upon  convergence  to  steady-state  conditions;  i.e.,  Eq.  (4.157) 
can  also  be  expressed  in  a  form  using  the  real  weights, 


IWI2 


ir0„<o)r 


R 

E^2(0„o+  e  ws 

r=l  *=t 


£  i>vs 


*= i 
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in  which  the  effect  of  the  magnitudes  of  the  weights  is  obvious. 

To  normalize  the  effect  of  adapted-weight  magnitude  level,  one  can  calculate  the 
degradation  Dsn  in  the  S/N  ratio.  This  is  simply  the  quiescent  S/N  ratio  divided  by  the 
adapted  S/N  ratio  and  leads  to  the  convenient  expression 


'a,2(8,)\./iy0(t)i2\ 


(4.159) 


where  the  ratio  in  the  second  term  is  Eq.  (4.157)  or  (4.158),  the  increase  in  output  noise 
power.  Thus,  we  simply  multiply  the  previous  performance  index  by  the  ratio  of  the 
power  pattern  values  in  the  direction  of  the  signal  6S . 

Covariance  Matrix,  Eigenvalues,  and  Eigenvectors— The  eigenvalues  and  eigenvectors 
of  the  covariance  matrix  are  evaluated  as  solutions  to  Eq.  (4.31): 

\M  -  (3,-2I|  =  0  and  Me,-  =  0,-2e,-,  (4.158) 

where  the  Hermitian  covariance  matrix  M  is  formed  as  indicated  in  Eqs.  (4.23)  through 
(4.28).  Receiver  noise  power  is  assigned  a  level  of  unity  for  convenience  in  computation, 
since  ail  noise  powers  are  expressed  as  ratios  to  receiver  noise  power.  Thus,  quiescent 
noise  matrix  Mq  becomes  an  identity  matrix,  and  the  individual  rth-jammer  covariance 
matrix  is  then  multiplied  by  its  power  ratio  Pn 

Mq  =  I  (4.159) 

R 

M  =  I  +  £>rM..  (4.160) 

r=l 

To  incorporate  bandwidth  into  the  interference  sources,  the  summation  in  Eq.  (4.160)  is 
further  refined  by  dividing  the  jammer  power  spectrum  into  a  number  of  discrete  spectral 
lines  as  described  in  Sec.  3.3.  We  assume  a  uniform  amplitude  spectrum  of  uncorrelated 
lines  spaced  apart  by  a  conotant  frequency  increment,  as  illustrated  in  Fig.  17.  If  Br 
denotes  the  percent  bandwidth  of  the  rth  jammer  and  Lr  its  total  number  of  spectrum 
lines,  the  power  ratio  and  frequency  offset  of  the  1th  line  are  given  by 


Prl  « 

A/i- 


ft)- 


power  ratio  of  spectrum  line 


»)  “  2  +  =  freqUenCy  0ffeet 


B, 

100 


(4.161) 

(4.162) 


ft 

4 


The  element  spacing  of  half  a  wavelength  in  Eq.  (4.144)  must  now  be  defined  in  terms 
of  wavelength  X0,  which  corresponds  to  the  RF  bandwidth  center  frequency  f0,  and  a 
new  phase  factor  for  the  ith  spectrum  line  is  obtained, 


d  = 


(4.163) 
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u"  *  To  f sin  e'  *  1  +  7T  I  si“ 8' 


■  (4 +  (m)  [-  i +  (£= t)]  )f sin  6'  •  (4-164) 


The  covariance  matrix  may  now  be  rewritten  as  a  summation  of  the  spectrum  line 
matrixes  Mrt, 


R  Lr 

M  =  I  +  L  LPrlMr 

r=  l  1=1 


1  e 


;2ur,  /  iurl 


(4.165) 


1  e' 


MW  =  [sSSrf]  “  g-/4«r,  e-/2 url  1 


(4.166) 


where 


Sr[k  =  exp  \iuri(2k  -K-l)], 


(4.167) 


with  url  defined  in  Eq.  (4.164)  and  Prl  defined  in  Eq.  (4.161).  If  m  denotes  the  row  and 
n  the  column,  the  mnth  component  of  the  Mr/  matiix  is  equal  to 


-  KmSr,J  - 


(4.168) 


Note  that  the  calculations  require  four  data  values  to  be  specified  for  each  jammer: 

Pr  =  ratio  of  total  jammer  power  to  receiver  noise  power, 

Qr  =  spatial  angle  location  off  boresight, 

Br  -  percent  bandwidth  of  spectrum, 

Lr  =  number  of  discrete  spectrum  lines. 

Since  the  system  cannot  respond  to  noise  outside  of  its  element  channel  receiver 
bandwidth  Bc ,  it  is  assumed  that  Br  would  not  be  chosen  to  exceed  Bc. 
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Having  formed  the  complete  covariance  matrix  for  some  chosen  distribution  of  inter¬ 
ference  sources  enables  us  to  solve  it  for  the  eigenvalues  and  eigenvectors,  using  the  For¬ 
tran  computer  program  CMPLXEIG.  This  program  was  developed  at  the  University  of 
Wisconsin  Computing  Center  [42], 


Program  CMPLXEIG  prints  out  the  eigenvalues  and  their  corresponding  eigenvectors. 
For  a  K-element  array,  there  will  be  K  eigenvalues  and  K  eigenvectors.  The  minimum 
value  0O2  that  the  eigenvalues  can  have  is  unity,  which  corresponds  to  receiver  noise 
power  level,  and  many  selected  distributions  o'i  sources  result  in  multiple  roots  (eigen¬ 
values)  equal  to  unity.  Such  unity  solutions  1  the  eigenvalues  are  called  nonunique* 
eigenvalues,  and  the  corresponding  eigenvectors  printed  out  are  generally  meaningless 
vectors.  The  useful  eigenvectors  printed  out  are  those  associated  with  the  unique  eigen¬ 
values  and,  for  CMPLXEIG,  these  are  defined  as  eigenvalues  greater  than  1.01,  which  cor¬ 
responds  to  an  eigenvector  beam  delivering  jammer  power  not  less  than  20  dB  below  the 
receiver  noise  level.  This  limit  for  qualifying  the  unique  eigenvalues  is  a  matter  of  judg¬ 
ment  and  may  be  selected  within  rather  wide  limits  of  perhaps  1.001  to  1.1  for  our 
purposes.  Based  on  this  limit  criterion,  the  unique  eigenvectors  corresponding  to  the 
unique  eigenvalues  are  culled  from  the  CMPLXEIG  output,  normalized  to  obtain  unit 
vectors  in  the  Hermitian  sense,  and  then  saved  in  a  data  file.  All  vector  data  output  cor¬ 
responding  to  the  nonunique  unity  eigenvalues  is  discarded  at  this  point;  Table  1  in 
Sec.  3.4  is  a  typical  illustration  of  the  unique  eigenvalue  and  unique  eigenvector  data 
saved  from  the  output  of  CMPLXEIG.  All  performance  characteristics  are  then  com¬ 
puted  from  the  saved  data. 


Although  the  nonunique  eigenvectors  are  unnecessary  or  calculating  system  per¬ 
formance,  situations  may  arise  in  which  it  becomes  desirable  to  operate  with  a  filled  Q 
matrix.  For  these  situations,  one  must  construct  the  missing  arbitrary,  nonunique  eigen¬ 
vectors.  A  convenient  approach  to  accomplishing  this  task  is  outlined  in  Appendix  B. 


Selected  Distribution  of  Interference  Sources— The  performance  characteristics  of 
several  selected  distributions  of  interference  sources  have  been  calculated  to  demonstrate 
the  behavior  of  our  eight-element  linear  adaptive  array.  The  distributions  selected  are 
listed  in  Table  2,  where  they  are  identified  by  case  symbols  A  through  H.  For  each  case, 
the  covariance  matrix  is  formed  as  described  above,  and  the  associated  unique  eigenvalues 
and  eigenvectors  are  computed.  From  these,  output  pattern  function  G(0,f)  is  evaluated 
by  using  Eq.  (4.119),  under  the  assumptions  discussed  above.  The  transient  performance 
is  usually  evaluated  from  the  increase  in  output  noise  power,  given  by  Eq.  (4.157)  or  a 
modificatic  thereof  to  include  bandwidth.  If  the  magni  j  of  the  adapted  weights 
changes  appreciably,  however,  transient  performance  is  evaluated  on  the  basis  of  the 
degradation  in  the  S/N  ratio,  given  by  Eq.  (4.159). 


Case  A  corresponds  to  a  single  narrowband  source  in  the  sidelcbe  region,  one  unique 
eigenvalue  (10,001),  and  one  unique  eigenvector.  This  case  was  discussed  in  considerable 
detail  in  Sec.  4.3  and  4.4,  with  the  adapted  pattern  shown  in  Fig.  31  and  the  transient 
performance  shown  in  Fig.  33.  It  should  be  emphasized  that  these  performance  plots  are 
very  sensitive  to  the  location  of  the  jammer  with  respect  to  the  quiescent  steered-beam 


♦It  is  possible,  though  rare,  to  get  nonunique  eigenvalues  (multiple  roots)  that  are  greater  than  unity. 
For  such  cases,  the  solutions  must  be  retained  as  if  they  were  unique. 
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Table  2 

Selected  Distributions  of  Interference  Sources 


Interference  Sources 

Number 

of 

Sources 

er 

os 

Location 

Angle 

(Deg) 

Bandwidth 

(%) 

Spectrum 

Lines 

Unique  Eigenvalues 


Number  Values 


1,250 

21 

0 

1,250. 

18 

0 

1,200 

22 

0' 

1,250 

18 

0 

125 

22 

0 

40 

18 

0 

125 

25 

0 

400 

33 

0 

1,250 

42 

0 

1,250 

42 

0 

1,250 

42 

2 

1,250 

42 

15 

1,250 

5 

0 

1,100 

36 

26 

1,100 

48 

21 

1,100 

66 

19 

1,100 

-66 

19 

1,100 

-48 

21 

1,100 

-36 

26 

1,100 

36 

26 

1,100 

48 

21 

1,100 

66 

19 
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pattern.  For  example,  if  the  single  jammer  happens  to  be  located  in  a  null  of  the  quiescent 
pattern,  then  even  though  the  one  degree  of  freedom  involved  has  been  consumed,  the 
single  eigenvector  beam  will  be  deleted  by  the  weighting  factor  Wql  =0,  the  adapted  pat¬ 
tern  will  suffer  no  change  from  quiescent,  and  there  will  be  no  transient  response.  The 
position  in  Fig.  31,  with  the  jammer  located  at  the  peak  of  the  urst  sidelobe,  represents 
a  typical  worst  case  for  sidelobe  jamming.  The  adapted  pattern  must  shift  its  sidelobes  to 
accommodate  the  jammer,  but  there  is  little  distortion  in  the  main  beam  because  only  one 
degree  of  freedom  of  the  available  seven  has  been  consumed. 


Case  B  has  two  narrowband  sources  in  the  sidelobe  region,  with  almost  equal  cower 
ratios  and  located  close  together;  two  unique  eigenvalues  (18,544  and  1,058);  and 
unique  eigenvectors  (see  Table  1).  This  case  was  discussed  in  detail  in  Sec.  4.4,  witi 
terns  shown  in  Figs.  34  and  35  and  transient  response  shown  in  Fig.  36.  Important  points 
demonstrated  by  this  case  include  the  widely  different  eigenvalues,  the  two  different 
eigenvector  beams,  and  the  relatively  slow  convergence  time  even  though  both  jammers 
are  strong. 


Case  C  is  the  same  as  Case  B  except  that  the  power  ratio  of  the  source  at  22°  is  re¬ 
duced  10  dB  below  that  of  the  source  at  18°;  there  are  two  unique  eigenvalues  (10,812 
and  190).  Note  again  that  the  ratio  of  the  eigenvalues  is  much  different  from  the  ratio  of 
the  two  jammer  powers.  Since  the  locations  of  the  two  sources  are  the  same  as  in  Case  B, 
the  steady-state  adapted  pattern  for  this  case  is  almost  exactly  as  in  Fig.  34.  The  unique 
eigenvectors,  however,  are  different  and  give  rise  to  different  retrodirective  eigenvector 
beams,  as  shown  in  Fig.  37.  Note  that  the  beam  shapes  remain  similar  to  those  in  Fig. 

35,  but  that  both  gj(0)  and  g'2(Q)  are  shifted  to  the  left  so  the  respective  peak  and  null 
fall  very  close  to  the  position  of  the  strongest  source.  This  provides  the  proper  power 
balance  for  achieving  decorrelation  between  their  outputs.  The  transient  response  shown 
in  Fig.  38  is  much  slower  than  that  of  Fig.  36  because  of  the  much  smaller  second 
eigenvalue. 


Case  D  has  four  narrowbai:*-  sources  in  the  sidelobe  region,  unequal  power  ratios,  and 
moderate  spacing;  there  are  fcmr  unique  eigenvalues  (11,616,  2,486,  406,  and  17).  The 
steady-state  adapted  pattern  is  shown  in  Figs.  39a  and  39b,  which  include  an  expanded 
scale  plot  in  the  vicinity  of  the  sources.  Here  we  note  that  the  four  sources  have  “cap¬ 
tured”  four  nulls,  or  four  degrees  of  freedom  out  of  seven  available,  and  that  this  large 
percentage  of  null  constraints  causes  an  appreciable  distortion  of  the  main  beam  as  well 
as  of  the  remaining  sidelobe  region.  This  is  a  good  illustration  of  the  fact  that  adaptive 
pattern  nulls  usually  do  not  aline  themselves  exactly  on  the  source  locations  unless  the 
sources  are  very  strong.  The  retrodirective  eigenvector  beams  are  illustrated  in  Figs.  40a 
and  40b  and  are  associated  with  the  eigenvalues  as  listed  below. 


Eigenvalue 

Eigenvector  Beam 

11,616 

2,486 

g'2$)  =  frff2*2(0) 

406 

g'Z(9)  =  fr,3*3(0) 

17 

=  #,4*40) 
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-10 


-15 


-20 


-25 


-30 


-35 


-40 
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(b)  Ream*  *g(0)  and  g'4{0) 

Fit.  40— Retrodireclive  eigenvector  beams  for  four  sources,  Case  D 
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Beam  ^(0)  receives  power  from  all  four  sources  but  functions  essentially  in  the 
manner  of  a  “sum”  beam  for  the  two  strongest  sources,  thus  accounting  for  the  largest 
eigenvalue.  Beam  g2  (#)  also  receives  power  from  all  four  sources,  but  it  functions  essen¬ 
tially  in  the  manner  of  a  “difference”  beam,  which  places  its  difference  null  so  as  to 
balance  off  the  three  weaker  sources  against  the  strong  source.  Beam.g3  (0)  essentially 
nulls  out  the  strongest  source  and  favors  the  summation  of  the  other  three.  It  is  a 
“single  lobe  in  notch”  type  of  pattern.  Beam  g4(0)  brackets  the  four  sources  with  a 
“double  lobe  in  notch”  pattern  in  which  the  two  strongest  sources  are  essentially  nulled 
out.  Its  low  gain  on  the  remaining  two  weak  sources  accounts  for  the  smallest  eigenvalue. 

Note  that  the  eigenvector  beams  become  progressively  more  complicated  and  less 
recognizable  as  to  their  function  as  we  proceed  toward  the  smaller  eigenvalues.  In  fact, 
it  seems  meaningless  to  refer  to  beam  g4  (6 )  as  a  retrodirective  beam  since  there  is  no 
main  lobe  as  such,  and  none  of  its  several  high  lobes  points  toward  the  sources.  How¬ 
ever,  the  dominant  consideration  is  still  its  characteristics  in  the  directions  of  the  sources 
so  as  .to  satisfy  Eqs.  (4.136)  and  (4.138),  and  in  that  sense  the  pattern  is  retrodirective. 

The  transient  response  for  this  case  is  shown  in  Fig.  41,  where  it  will  be  noted  that 
Ds„  is  the  ordinate  because  the  weight  magnitude  drops  about  26%  for  this  case.  Note 
also  that  the  time  scale  is  a  combined  linear  and  logarithmic  scale  with  the  transition 
point  at  t  -  6  ps.  The  log  scale  permits  better  assessment  of  the  long  convergence  time 
caused  by  the  two  smallest  eigenvalues.  Also  note  that  the  steady-state  degradation  in 
S/N  Is  about  1.5  dB. 


Fig.  41— Transient  response  for  four  sources,  Case  D;  eigenvalues  11,616,  2,486,  406,  and  17 
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Case  E  has  a  single  source  in  the  sidelobe  region,  42°  off  array  boresight;  various 
bandwidths  are  used  to  demonstrate  the  effect  of  this  parameter  on  our  frequency- 
sensitive  adaptive  array. 

El  (zero  bandwidth)  has  one  unique  eigenvalue  (10,001).  This  narrowband  case  is 
the  same  as  Case  A  except  for  source  location  angle,  and  the  two  cases  may  be  compared 
to  illustrate  the  dependence  on  local,  m.  The  steady -state  adapted  pattern  h.  shown  in 
Fig.  42a  and  serves  as  the  starting  point  in  this  bandwidth  series  of  patterns.  As  in  Case 
A,  the  source  “captures”  only  one  degree  of  freedom  (one  null),  corresponding  to  the 
single,  unique,  uniform-illumination  eigenvector  beam. 

Case  E2  (2-percent  bandwidth)  has  two  unique  eigenvalues  (9,986  and  16).  The 
steady-state  adapted  pattern  is  shown  in  Fig.  42b,  where  it  is  evident  that  a  significant 
change  has  occurred  in  the  sidelobe  region.  Even  with  this  small  amount  of  bandwidth, 
the  single  source  has  captured  two  degrees  of  freedom  (two  nulls),  and  the  corresponding 
two  unique  eigenvector  beams  would  be  similar  to  those  shown  in  Fig.  35  for  Case  B, 
i.e.,  a  “sum”  and  a  “difference”  beam.  A  helpful  property  of  this  case  is  that  the  two 
eigenvector  beams  function  effectively  up  to  a  source  bandwidth  of  approximately  10 
percent  under  the  conditions  assumed  here,  so  that  the  extra  captured  degree  of  freedom 
results  in  considerable  adjustment  accommodation  to  changes  in  bandwidth.  It  provides 
the  mechanism  for  obtaining  a  deep  sidelobe  notch  of  variable  width. 

Case  E3  (15-percent  bandwidth)  has  three  unique  eigenvalues  (9,529,  469,  and  5). 
The  steady-state  adapted  pattern  is  shown  in  Figs.  43a  and  43b,  which  include  an  ex¬ 
panded  scale  plot  in  the  vicinity  of  the  source.  The  equivalent  spatial  extent  of  the 
source  is  indicated  for  its  15%  bandwidth.  The  single  source  has  now  captured  three 
degrees  of  freedom  (three  nulls)  out  of  seven  available,  and  the  resulting  distortion  in  the 
adapted  pattern  is  very  evident.  The  corresponding  three  unique  eigenvector  beams  are 
shown  in  Fig.  44  and  are  associated  with  the  eigenvalues  as  follows. 


Eigenvalue 

9,529 

469 

5 


Eigenvector  Beam 

4(0)  =  *,1*1(0) 

4(0)  =  *,2*2(0) 

4(0)  =  *,3*3(0) 


The  “sum”  and  “difference”  beams,  g^(9)  and  4(0)  respectively,  are  representative  for 
any  of  the  bandwidths  tested  in  this  series,  from  2%  to  20%;  i.e.,  this  pair  of  beams  re¬ 
mains  practically  invariant.  Beam  g'^(0)  is  a  “deep-notch”  pattern  which  is  similar  to  its 
counterpart  in  Fig.  40b  except  that  the  single  lobe  within  the  notch  has  a  peak  of  -46 
dB  and  thus  does  not  show  here.  This  third  eigenvector  beam  comes  into  play  at  a 
bandwidth  of  roughly  7%  and,  except  for  the  notch  region  and  its  lobe,  the  pattern 
undergoes  very  little  change  for  the  bandwidths  tested,  up  to  20%. 

The  transient  responses  for  four  different  values  of  source  bandwidth  are  plotted  in 
Fig.  45.  In  general,  convergence  rate  slows  down  as  bandwidth  increases,  because  of  the 
small  eigenvalues  generated.  To  compute  these  curves  from  the  basic  output  noise  power 
equation  (Eq.  (4.159)),  we  modified  the  source  power  summation  to  incorporate  a 
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summation  of  spectrum-line  contributions  over  the  bandwidth  of  the  source,  in  the  same 
manner  as  discussed. in  the  previous  subsection.  Thus, 

!  R 

source  power  =  ^  PriGt2(6r,  t)  (4.169) 

|  r  =  l  /=  1 

} 

( 

where  Prl  is  the  spectrum-line  power  ratio  given  by  Eq.  (4.161)  and  G/(0r,f)  is  the  out¬ 
put  pattern  function  for  the  /th  spectrum  line  which  requires  the  modified  phase  factor 
|  urh  given  by  Eq.  (4.164),  for  its  evaluation. 

Case  F  has  a  single  narrowband  source  in  the  main-beam  region,  6°  off  array  bore- 
sight;  there  is  one  unique  eigenvalue  (10,001).  This  case  demonstrates  the  effects  of 
main-beam  jamming,  but  one  should  keep  in  mind  that  it  is  the  same  as  Case  A  and 
Case  El  except  for  source  location  angle.  The  source  captures  only  one  degree  of  free¬ 
dom  (one  null)  corresponding  to  the  unique,  uniform-illumination  eigenvector  beam.  The 
steady-state  adapted  pattern  is  shown  in  Fig.  46,  where  it  is  evident  that  main-beam  jam¬ 
ming,  produces  rather  severe  distortions  in  the  output  pattern.  An  inspection  of  Eq. 
(4.125),  which  specializes  the  output  pattern  function  to  the  single  narrowband  source 
case,  reveals  that  if  the  main-beam  steering  direction  90  happens  to  become  alined  with 
the  source  direction  0lt  then  the  main  beam  may  practically  disappear.  In  fact,  if  the 
quiescent  pattern  happens  to  be  of  uniform  illumination,  as  we  assumed  for  these  calcu¬ 
lations  (Eq.  (4.144)),  then  the  entire  output  pattern  function  would  indeed  disappear  (go 
to  zerc)  for  0O  =  6l. 


-90  -75  -60  -45  -30  -15  0  15  30  45  60  75  90 

SPATIAL  DEGREES 


Fig.  46— Steady-state  adapted  pattern  for  single  source,  main-beam  jammer.  Case  F 
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The  transient  response  is  shown  in  Fig.  47a.  There  are  two  characteristics  to  note, 
peculiar  to  main-beam  jamming.  The  first  is  that  the  transient  commences  from  a  high 
power  level  for  91  close  to  d0,  since  the  source  gets  the  gain  of  the  quiescent  main  bom 
at  t  =  0.  The  second  is  that  the  increase  in  output  noise  power  typically  drops  below 
quiescent  (0  dB)  in  the  steady  state,  since  the  adapted  weights  are  attenuated  for  0! 
close  to  0q.  The  reason  for  the  attenuation  of  the  weights  is  that  the  unique  eigenvector 
is  closely  alined  with  the  quiescent  weight  vector,  and  this  leads  to  small  or  zero  magni¬ 
tudes  for  the  Wk  in  Eq.  (4.65),  which  in  turn  results  in  small  magnitudes  for  the  real  Wk 
weights.  Such  attenuation  causes  the  increase  in  output  noise  power  to  be  unsatisfactory 
as  a  performance  index  for  main-beam  jamming  cases,  so  that  it  becomes  desirable  to 
consider  the  effect  of  the  attenuation  on  the  signal  as  well. 

Signal  effects  may  be  incorporated  by  computing  the  degradation  that  occurs  in  the 
S/N  ratio  Dsn,  formulated  in  Eq.  (4.159).  Figure  47b  plots  Dsn  for  this  case  for  ds  =  0°, 
and  we  see  that  there  is  a  net  steady-state  degradation  of  about  4.5  dB  in  output  S/N 
ratio.  Thus,  although  the  output  noise  has  decreased  below  quiescent,  the  output  signal 
power  drops  even  more  and  leads  to  a  net  degradation  in  S/N  ratio  for  main-beam  jamming. 

Case  G  has  three  sources  in  the  sidelobe  region;  the  locations  and  bandwidlhs  have 
been  chosen  to  result  in  complete  coverage  of  the  sidelobe  region  from  30°  to  90°.  This 
is  equivalent  to  spreading  out  33  narrowband  sources  rather  uniformly  in  sin  9  spacing 
over  this  sidelobe  region.  Six  unique  eigenvalues  are  associated  with  this  case.  Although 
this  case  may  not  represent  a  practical  interference  situation,  it  demonstrates  the  remark¬ 
able  effectiveness  with  which  the  adaptive  array  uses  its  degrees  of  freedom  to  cope  with 
such  widespread  interference.  Furthermore,  it  illustrates  clearly  that  when  interference 
sources  are  grouped  in  a  continuous  distribution,  the  eigenvector  beams  may  be  charac¬ 
terized  by  a  family  of  harmonic  pattern  “modes.”  This  case  involves  the  first  six  modes 
of  the  set. 

The  steady-state  adapted  pattern  is  shown  in  Figs.  48a  and  48b,  which  include  an 
expanded-scale  plot  demonstrating  the  remarkably  low  sidelobe  level  achieved  throughout 
the  entire  jamming  region  from  30°  to  90°.  With  only  seven  degrees  of  freedom  avail¬ 
able,  the  array  cannot  respond  to  the  33  interference  sources  on  an  ind'vidual  null  basis, 
but  it  can  respond  on  a  resolution  basis  because  of  the  close  spacings  of  the  sources. 

Thus,  it  depresses  that  entire  sidelobe  region  by  efficiently  using  six  degrees  of  freedom, 
or  six  eigenvector  beams. 

Previous  plots  of  the  eigenvector  beams  have  been  shown  in  relationship  to  the 
quiescent  pattern  and,  in  fact,  were  adjusted  in  power  level  to  the  quiescent  pattern  via 
weighting  factors  Wqi;  i.e  ,  the  patterns  plotted  have  been  g!(9 )  =  IVe^igi(9).  However,  the 
eigenvector  beams  g,(9)  are  determined  only  by  the  covariance  matrix  and  are  completely 
lependunt  of  the  quiescent  pattern.'  Furthermore,  these  beams  are  characterized  by 
c  ...  .•>.  shapes  or  modes  that  constitute  a  sort  of  harmonically  related  family.  To  stress 
i>  \bove  two  points  and  also  to  bring  out  the  role  of  the  Wqj  weighting,  the  eigenvector 
L'.  ams  for  this  case  are  plotted  as  gf(9)  directly  and  the  Wqi  magnitudes  are  listed  sepa¬ 
rately  in  the  following  tabulation,  as  associated  with  the  six  unique  eigenvalues. 
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(a)  Normal  scale  plot 


Fig.  48— Steady-state  adapted  pattern  for  cont  nusu:  Jistribution  sources,  Case  G 
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i  Unique 

Eigenvector 

IV 

•  Eigenvalue 

Bexun 

Weighting 

i  13,316 

8l(0) 

0.083 

:  9,692 

82$) 

0.416 

;  3,091 

83(0) 

0.07 

296 

84(0) 

1.03 

10.5 

g5(0) 

0.64 

1.16 

86(0) 

1.04 

The  six  eigenvector  beams  are  plotted  in  Figs.  49a -49c,  and  their  mode  characteris¬ 
tics  and  Wq{  weight  relationship  will  be  briefly  discussed.  Also,  their  role  in  the  transient 
response  shown  in  Fig.  50  will  be  pointed  out. 

Beam  gi(6)  is  the  sum  beam  mode  with  no  nulls  inside  the  source  region.  This 
beam  sums  all  33  sources  with  appreciable  gain,  thus  resulting  in  the  largest  eigenvalue, 
and  has  the  fastest  transient  response.  However,  note  that  the  $qi  weighting  is  only 
0.083,  which  relegates  this  high-power  beam  to  a  relative  ly  minor  role  in  determining 
transient  response  and  adapted  pattern. 

Beam  g2(d)  is  the  “difference”  beam  mode  with  one  null  inside  the  source  region. 
Except  for  the  null  region,  this  beam  also  sums  the  sources  with  appreciable  gain  and 
results  in  the  second  largest  eigenvalue  and  the  second  fastest  transient  response.  The 
Wq(  weighting  is  0.416,  which  means  that  this  high-power  beam  has  considerable  effect 
upon  both  the  transient  response  and  adapted  pattern.  It  accounts  for  most  ox  the  initial 
fast  drop  in  the  transient  response. 

Beam  g3(d)  is  a  “single  lobe  in  notch”  type  with  two  nulls  in  the  source  region.  This 
beam  is  obviously  of  lower  gain  than  the  first  two  and  results  in  the  third  largest  eigen¬ 
value  of  3,091.  Note  the  approximate  alinement  of  this  beam  with  the  first  two.  Its 
Wq 3  weighting  of  only  0.07  drops  it  into  a  minor  role,  even  though  it  carries  appreciable 
power.  It  is  interesting  at  this  point  to  look  at  the  adapted  pattern  at  time  t  =  8  ps, 
shown  in  Fig.  51a,  whic'.  incorporates  the  contributions  of  g\(0),  g2{0),  and  g3 (6).  Note 
that  these  first  three  eigenvector  beams  have  already  reduced  the  sidelobe  level  to  -30  dB 
for  most  of  the  source  region. 

Beam  g4(C)  is  a  “double  lobe  in  notch”  type  with  three  nulls  in  the  source  region. 
This  beam  is  of  low  gain  in  the  source  region  and  results  in  the  modest  eigenvalue  of  296. 
However,  note  that  the  Wq4  weighting  is  a  strong  1.03,  which  essentially  places  this  beam 
in  control  of  the  transient  response  after  the  initial  fast  drop  caused  by  g2(0).  The  slow 
decay  so  evident  in  Fig.  50  is  dominated  by  this  one  beam.  It  also  carried  the  adapted 
pattern  quite  far  along  toward  its  steady  state,  as  shown  in  Fig.  51b  for  time  t  =  120  Ms, 
which  incorporates  the  contributions  of  the  first  four  eigenvector  beams. 

Beam  £s(0)  is  a  “triple  lobe  in  notch”  type  with  four  nulls  in  the  source  region.  As 
shown  in  Fig.  49c,  this  beam  is  of  very  low  gain  in  the  source  region  and  results  in  the 
small  eigenvalue  of  10.5.  It  does  have  a  strong  Wqb  weighting  factor  of  0.64,  however, 
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Fig.  49— Retrodirective  eigenvector  beam#  for  continuous  distribution  sources,  Case  G 
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(c)  Beams  g5(0)  and  s6(0) 


Fig.  49  (Continued)— Retrodirect've  eigenvector  beams  for  continuous  distribution  sources,  Case  G 
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Fig.  50— Transient  response  for  continuous-distribution  sources,  Case  G 
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(b) t  =  120  ms 

Pig.  51~Adaptiv6  patterns  at  two  different  times  for  continuous-distribution  sources, 
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which  gives  it  appreciable  influence  in  the  low  sidelobe  region,  so  that  it  causes  the 
adapted  pattern  to  slowly  converge  very  close  to  steady  state.  It  requires  some  3000  (is 
for  its  transient  response. 


Beam  <g6(0)  is  a  “four  lobes  in  notch”  type  with  five  nulls  in  the  source  region.  As 
shown  in  Fig.  49c,  this  beam  is  of  extremely  low  gain  in  the  source  region  and  results  in 
an  eigenvalue  of  1.16,  barely  above  quiescent  noise  level.  This  small  eigenvalue  causes 
the  servo  gain  term  in  Eq.  (4.119)  to  have  a  value  of  only  0.075  and  results  in  an  almost 
negligible  contribution  from  the  beam,  even  though  its  Wq6  weighting  factor  is  a  strong 
1.04.  For  all  practical  purposes,  this  sixth  eigenvector  beam  could  be  dispensed  with, 
which  means  that  the  array  is  essentially  devoting  only  five  degrees  of  freedom  to  coping 
with  the  widespread  interference. 


Case  H  has  six  sources  in  the  sidelobe  region,  in  which  the  locations  and  bandwidths 
have  been  chosen  to  result  in  complete  coverage  of  the  sidelobe  region  from  30°  to  90° 
on  each  side.  It  is  essentially  an  application  of  the  previous  case’s  source  distribution  to 
both  sides,  to  demonstrate  how  the  array  can  take  advantage  of  symmetry  in  the  source 
distribution  to  better  use  its  limited  degrees  of  freedom  in  coping  with  a  jamming  situa¬ 
tion  which  occupies  60%  of  its  spatial  coverage.  All  seven  degrees  of  freedom  are 
consumed. 


The  steady-state  adapted  pattern  is  chown  in  Figs.  52a  and  52b,  which  include  an 
expanded-scale  plot,  again  demonstrating  the  remarkably  low  sidelobe  level  achieved 
throughout  the  entire  jam  ning  region  from  30°  to  90°  on  both  sides  (the  patterns  are 
symmetrical  for  this  case).  As  in  the  previous  case,  the  adaptive  array  is  responding  to 
the  widespread  jamming  by  depressing  the  entire  sidelobe  region.  The  transient  response 
for  this  case  is  chown  in  Fig.  53. 


A  complete  set  of  eight  unique  eigenvalues  is  associated  with  this  case.  They  are 
listed  below,  together  with  the  eigenvector  beam  identifications  and  the  Wqi  weighting 
magnitudes. 


Unique 

Eigenvector 

iwy 

Eigenvalues 

Beam 

Weighting 

13,532 

*10) 

7-10-14 

13,386 

*20) 

0.024 

12,619 

*30) 

8*10— 14 

9,682 

*4  0) 

0.60 

3,224 

*50) 

2,10— 14 

350 

*60) 

1.42 

14,7 

*70) 

3-10-12 

1.18 

*8(9) 

2.37 
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SPATIAL  DEGREES 
(a)  Normal  scale  plot 
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(b)  Expanded  sidelobe  plot  in  source  region 


Fig.  52— Steady-state  adapted  pattern  for  continuous-distribution  sources,  Case  H 
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Fig.  53— Transient  response  for  continuous-distribution  sources,  Case  H 


Note  the  severe  attenuation  in  four  of  the  weighting  factors,  which  completely 
eliminates  any  contribution  from  those  beams.  Nonetheless,  all  eight  eigenvector  beams 
are  shown  in  Figs.  54a  and  54b  and  will  be  briefly  discussed  because  of  their  departure 
from  the  usual  harmonic-mode  series  which  have  characterized  all  of  the  previous  cases. 
The  reason  for  the  departure  is  that  the  interference  sources  are  not  grouped  m  one  con¬ 
tinuous  distribution,  but  are  split  into  two  distributions. 

Beams  gj.  (0),  #2(0)>  Ss (0)  are  of  the  sum  mode  type,  each  having  appreciable  major- 
lobe  gain  in  the  source  regions,  thus  resulting  in  large  eigenvalues.  However,  note  that  the 
Wqi  weighting  factors  are  so  small  that  these  high-power  beams  are  eliminated  from  the 
transient  response. 

Beam  g4(0)  is  of  the  difference  mode  type,  with  appreciable  gain  in  the  source  re¬ 
gions  to  result  in  a  large  eigenvalue  of  9,682.  The  weighting  is  also  strong,  0.60, 
so  that  this  beam  alone  controls  the  initial  fast  drop  in  the  transient  response  of  Fig.  53. 
The  effect  of  this  single  beam  on  the  adapted  pattern  is  shown  in  Fig.  55,  which  is  plotted 
at  time  t  =  8  /is,  and  it  is  seen  that  sidelobes  are  already  reduced  to  a  level  of  about  -28 
dB  by  this  time. 

Beam  g5(0)  is  of  the  “double  lobe  in  notch”  type,  with  modest  gain  in  the  source 
regions  resulting  in  an  eigenvalue  of  3,223.  Like  the  first  three  beams,  its  weighting  factor 
is  so  small  as  to  eliminate  it  from  the  transient  response. 
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Fig.  55— Adaptive  pattern  at  /  =  8  H&  for  continuous-distribution  sources,  Case  H 

Beam  g6(0)  is  of  the  “double  lobe  in  notch”  type,  with  low  gain  in  the  source  re¬ 
gions  resulting  in  a  modest  eigenvalue  of  350.  However,  its  w  sighting  factor  is  a  strong 
1.42,  which  places  this  beam  in  complete  control  of  the  transient  response  after  the  initial 
fast  drop  due  to  g4(t?j.  The  slow  decay  noted  in  Fig.  53  is  caused  entirely  by  this  one 
beam,  and  after  its  decay  time  of  approximately  120  ns  the  adapted  pattern  will  be 
virtually  at  its  steady-state  condition. 

Beam  g7($)  is  of  the  “triple  lobe  in  notch”  type  with  very  low  gain  in  the  source 
regions,  resulting  in  the  small  eigenvalue  of  14.7.  Us  weighting  factor  is  very  small,  so  it 
plays  no  part  in  the  transient  response. 


,  Kt  - 


Beam  g8{0)  is  another  “triple  lobe  in  notch”  type  with  an  extremely  low  gain  in  the 
source  regions,  resulting  in  an  eigenvalue  of  1.18,  barely  above  quiescent  noise  level.  Be¬ 
cause  this  eigenvalue  is  so  close  to  unity,  the  servo  gain  term  in  Eq.  (4.119)  will  have  a 
value  of  only  0.083,  resulting  in  almost  negligible  contribution  to  the  transient,  even 
though  its  weighting  factor  is  a  strong  2.37. 


Although  this  beam  does  not  contribute  to  the  adapted  pattern  as  expressed  by  Eq. 
(4.119),  the  reader  should  consider  the  fact  that  for  this  case  the  interference  source  dis¬ 
tribution  has  consumed  all  seven  available  degrees  of  freedom,  and  jammer  power  is  being 
delivered  by  all  of  the  eigenvector  beams  except  g%(0).  If  one  observes  the  Wqi  weighting 
for  each  beam  and  applies  Eq.  (4.65)  for  the  Wk  orthonormal  weights,  g8(0)  will  be  the 
only  beam  capable  of  delivering  signal  power  to  the  output  of  the  orthonormal  system, 
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Fig.  56— Quiescent  sum  and  difference  tracking  patterns 


and  therefore  it  must  be  identical  to  the  steady-state  adapted  pattern  shown  in  Fig.  52; , 
i.e.,  for  this  case,  G(0,°°)  «  Wq8ga(d). 

Adaptive  Tracking  Patterns— It  is  of  interest  to  see  what  happens  to  typical  tracking 
patterns  when  both  the  sum  and  the  difference  pattern  are  instrumented  to  adapt  to  an ' 
interference  environment.  For  the  quiescent  patterns,  tapered  sum  and  difference  illumina¬ 
tions  were  chosen  as  follows: 
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These  illuminations  result  in  the  quiescent  patterns  shown  in  Fig.  56,  which  are  plotted  on 
a  linear  voltage  scale  to  emphasize  the  slope  and  crossover  associated  with  the  difference 
pattern. 
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From  Table  2,  Case  D  was  chosen  as  an  appropriate  sidelobe  interference  distribu¬ 
tion.  It  includes  four  narrowband  sources  of  unequal  power  ratios  and  moderate  spacing 
and  results  in  the  four  unique  eigenvalues  11,616,  2,486,  406,  and  17.  Its  associated 
eigenvector  beams  were  as  shown  in  Fig.  40. 

A  time  sequence  for  the  adapting  patterns  is  shown  in  Figs.  57a  -57d  for  times 
t  =  12,  120,  and  1,200  jus,  and  for  the  steady  state.  At  time  t  =  12  ms,  the  sidelobe 
structure  has  changed  considerably  on  the  right-hand  side,  but  the  main-beam  crossover 
region  has  not  yet  been  affected  much,  This  time  would  incorporate  the  effects  of  the 
first  two  eigenvector  beams  g1(0)  and  g2(®)- 

At  time  t  =  120  Ms,  which  would  incorporate  the  effect  of  the  third  eigenvector 
beam  g3(0),  note  that  the  negative  half  of  the  difference  pattern  is  collapsing.  This  col¬ 
lapse  changes  both  the  crossover  point  and  the  slope  at  crossover.  In  addition,  the  sum 
pattern  is  distorted  and  shifted  to  the  left. 

At  time  t  =  1,200  ms,  which  incorporates  the  effects  of  all  four  eigenvector  beams, 
the  collapse  of  the  difference  pattern  is  almost  complete,  the  crossover  point  is  4°  off, 
and  the  slope  has  changed  considerably. 

At  steady-state  adapted  conditions,  the  distortion  of  the  difference  pattern  is  so 
severe  as  to  render  it  useless. 

This  example  illustrates  that  tracking  patterns  may  undergo  severe  deterioration  when 
they  are  made  adaptive  to  sidelobe  interference.  One  method  for  coping  wit!  this  serious 
problem  is  to  incorporate  constraints  on  the  patterns  in  their  crossover  region  by  sacri¬ 
ficing  array  degrees  of  freedom.  For  example,  one  might  sacrifice  a  degree  of  freedom  in 
order  to  force  the  difference  pattern  to  always  maintain  a  fixed  crossover  position  regard¬ 
less  of  sidelobe  adaptation. 


4.6.  Hard-Limiter  Modification 


Alw  dynamic  range  and  transient  response  time  of  the  configuration  in  Fig.  24  can 
be  improved  by  incorporating  a  hard  limiter  or  fast  AGC  (automatic  gain  control)  in  the 
conjugate  signal  branches,  in  the  same  manner  as  described  in  Sec.  3.4  for  the  single 
adaptive  servo  loop.  The  modified  circuit  schematic  is  shown  in  Fig.  58.  Note  that  each 
servo  loop  is  arranged  in  exactly  the  same  manner  as  the  single  loop  of  Fig.  20.  Refer¬ 
ence  32  is  recommended  for  a  more  rigorous  mathematical  development  of  the  effects  of 
envelope  limiting  in  adaptive-array  control  loops. 

The  modification  changes  tb*  output  voltages  from  the  correlation  mixers,  in  that 
the  amplitude  variations  ir  the  conjugate  signals  are  removed  and  only  the  phase  varia¬ 
tions  are  retained;  i.e.,  Eq.  (4.15)  becomes 


=  kH\ 


(4.171) 
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Fig.  58— Hard-limiter  modification  of  linear  adaptive 
array  configuration 


where  h  is  an  amplitude  constant  associated  with  the  limiter  arrangement,  assumed  to  be 
equal  for  all  channels.  From  Eq.  (4.25)  the  |2?fc  |  are  given  by 


\Ek\  =iM>l2  +  £  Mil2  ; 


(4.172) 


i.e.,  the  rms  voltage  magnitudes  are  to  be  equal  for  all  channels.  Under  these  assumptions, 
one  finds  that  the  elements  of  the  new  covariance  matrix  differ  from  the  elements  of  M 
only  by  a  common  factor.  This  leads  to  the  following  equation,  equivalent  to  Eq.  (4.20): 


T°^  +  [,+  (ife)M]w‘?’ 


(4.173) 


where  y'  =  kzhG'.  Operating  upon  Eq.  (4.173)  with  the  Q- matrix  transformation  as  in 
Eqs.  (4.54)  through  (4.57)  results  in 
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dW 
r°  ~dt 


1  + 


y'fii 

\Ek\. 


8;, 


W  =  B* 


(4.174) 


Recall  the  equivalent  orthcnormal  adaptive  array  system  of  Fig.  25  and  note  that  Eq. 

(4.174)  would  be  tantamount  to  defining  a  new  input  signal  vector  E'  for  that  system, 
whereby 


(4.175) 


i.e.,  the  square  root  of  a  limiting  operation  is  applied  uniformly  to  all  element  signals 
before  entering  the  Q-matrix  transformation  network,  and  there  is  no  limiting  in  the 
orthonormal  servo  loops  themselves. 


Corresponding  to  Eq.  (4.45)  we  have  a  new  equivalent  servo  gain  factor  from  Eq. 
(4.174), 


(4.176) 


This  is  the  important  result  of  the  limiter  modification,  because  we  are  now  back  to  the 
form  of  Eq.  (4.47)  and  can  use  the  solutions  thereof  if  p,-  is  replaced  by  jtt).  From  Eq. 

(4.176)  define  a  quiescent  servo  gain  factor  n'Q  as  in  Eq.  (4.60), 


4>  = 


(4.177) 


where  we  anticipate  that  large  values  will  be  selected  for  Uq  (via  amplifier  gains)  for  the 
same  reasons  as  discussed  in  Sec.  3.4.  From  the  ratio  of  fxv.  to  n'0,  one  can  obtain  a 
more  convenient  form  for  Eq.  (4.176), 


(4.178) 


where  Pr  is  the  .usual  ratio  of  source  power  to  receiver  noise  power.  Comparing  this  ex¬ 
pression  against  the  one  in  Eq.  (4.142),  we  note  that  the  only  difference  is  the  square 
root  of  the  power  ratio  summation  in  the  denominator.  That  summation,  however,  pro¬ 
duces  a  number  of  interesting,  effects,  including  the  following: 


1.  The  strong  jammer  power  ratios  will  dominate  the  summation  and  result  in  lower 
values  for  any  nt  associated  with  the  weak  jammers. 
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2.  For  small  eigenvalues,  n\  will  generally  be  smaller  than  ju'o  and  can  approach  unity. 
This  effect  results  in  rather  large  magnitude  levels  for  the  adapted  weights  and  requires  a 
modification  of  the  eigenvector  beam  concept  equations. 

3.  For  maximum  power  summations,  the  transient  performance  with  limiters  will  be 
no  better  than  without  the  limiters  and  may  be  worse. 

4.  It  further  complicates  the  maximum  servo  gain  factor  condition  which  relates  to 
control-loop  noise. 

Equation  (4.146)  is  now  modified  to  the  expression 


u'  = 


Pm2\  _  (*BcTq\ 

ft)  V  10  / 


(3.179) 


A  little  thought  on  the  several  parameters  involved  in  this  expression  indicates  that 
it  will  be  difficult  to  make  exact  comparisons  with  the  transient  performance  characteris¬ 
tics  plotted  previously.  From  Table  2,  we  see  that  the  largest  eigenvalue  encountered  was 
18,544  in  Case  B,  with  a  total  power  ratio  of  2,450.  Using  these  values,  together  with  a 
reasonably  high  value  of  100  for  leads  to  a  preliminary  selection  for  (d'm  of  about 
40,000.  If  the  channel  bandwidth  Bc  is  kept  at  5  Me  as  before,  r0  must  be  increased  to 
25,500  ps.  Summarizing  these  selected  constants,  then,  we  have 

H'o  =  ioo 

H'm  *  40,000 

<• 

Tq  ~  25,500  (is 

Bc  =  5  Me. 

Before  calculating  transient  response,  it  is  pertinent  to  examine  the  steady-state  ortho¬ 
normal  weight  equation  (Eq.  (4.62)),  which  will  now  be  given  by 


JL i!% 


Note  that  since  the  value  of  /d-  can  be  much  smaller  chan  p'0,  it  is  evident  that  the 
ortbonormal  weights  can  be  much  larger  than  the  quiescent  weights,  and  this  can  carry 
over  into  the  final  real  weights.  Therefore,  the  usual  previous  performance  index  of  in¬ 
crease  in  output  noise  power  cannot  be  used  with  the  limiter  modification.  It  is  neces¬ 
sary  to  use  the  degradation  in  S/N  ratio,  defined  in  Eq.  (4.159),  instead.  Based  on  this 
latter  performance  index,  the  following  cases  from  Tabic  2  are  compared: 

Case  A— The  transient  performance  for  this  case  of  a  single,  narrowband  interference 
source  is  greatly  improved  by  the  limiter  modification,  as  shown  in  Fig.  59  for  jammer 
power  ratios  of  21  and  31  dB.  For  such  single-source  cases,  the  system  will  behave  in 
much  the  same  manner  as  the  single  adaptive  loop  discussed  in  Sec.  3.4,  because  only  one 
unique  eigenvalue  and  one  power  ratio  are  involved.  Note  particularly  the  improvement 
at  the  lower  power  ratio  of  21  dB. 
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Fig.  59— Transient  response  with  and  without  limiter  modification  for  single  source, 
Case  A,  for  power  ratios  of  21  dB  and  31  c.3 


i 

! 


To  illustrate  the  increase  in  the  adapting  weights  referred  to  above,  Fig.  60  is  a 
special  plot  of  the  increase  in  output  noise  power,  in  which  the  time  scale  becomes  loga¬ 
rithmic  after  t  -  7.2  ms.  After  the  initial  fast  transient  in  which  the  jammer  noise  power 
is  nulled,  the  weights  increase  in  magnitude  slowly  (while  keeping  the  jammer  nulled  out) 
until  they  reach  the  steady-state  level  given  in  Eq.  (4.1#°).  For  the  seven  nonunique  eigen¬ 
values,  n\  =  °-83;  the  ratio  in  Eq.  (4.180)  is  equal  to  26.4.  approximately  representing 
the  steady-scate  weight  magnitudes;  the  square  of  that  ratio  then  accounts  for  the  2S-dB 
increase  in  output  noise  power.  The  array  gain  also  increases  by  the  same  amount,  so 
that  the  effect  does  not  show  up  in  the  S/N  plots  of  Fig.  59. 


Case  B— The  transient  performance  for  this  two-source  case  is  shown  in  Fig.  61.  It 
is  the  same  with  limiters  as  without  (compare  with  Fig.  36).  The  reason  for  the  practically 
identical  behavior  is  that  this  case  was  the  basis  for  the  choice  of  nm ,  t0,  }i'm ,  and  t'q, 
and  results  in  almost  identical  transient  decay  factors.  For  example,  for  the  largest 
eigenvalue, 


1  +  _  18,545 

T0  ”  12,750 


1.455 


O'. 


m 


1  +  uf 

1 _ ZJS  _  37,458 

r0  ”  25,500 


1.469. 


The  only  major  difference  is  that  with  the  limiter  modification,  the  magnitudes  of  the 
weights  change  in  a  manner  similar  to  that  described  for  Case  A. 
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Case  C— ' This  two-source  case  is  the  same  as  Case  B,  except  that  the  power  ratio  of 
the  second  source  is  reduced  by  10  dB.  Transient  response  is  shown  in  Fig.  62  with  and 
without  limiters.  Note  that  there  is  little  difference,  a  reflection  of  the  fact  that  the 
transient  decay  factors  are  not  much  different. 


Case  D— This  four-source  case  has  a  transient  response  with  limiters  that  is  so  nearly 
identical  to  that  shown  in  Fig.  41  that  there  is  no  point  in  showing  it.  The  transient 
decay  factors  are  almost  equal  for  the  two  conditions. 


Case  E— The  single-source  case  with  various  bandwidths  shows  little  difference  with 
limiters.  Figure  63  illustrates  the  transient  responses  for  a  bandwidth  of  15%. 


Case  F~This  case  of  main-beam  jamming  is  little  different  from  that  shown  in  Fig. 
47b,  but  it  shows  distinctly  faster  transient  response  with  limiter  modification  if  the 
power  ratio  were  lowered  by  10  dB  or  more,  in  a  manner  similar  to  Case  A. 


Case  G— This  extensive  sidelobe-jamming  case  with  limiters  has  a  transient  response 
nearly  identical  to  that  given  in  Fig.  50. 


Case  H— With  symmetry  in  its  extensive  sideiobe  jamming  and  twice  the  total  power 
of  Case  G,  this  case  shows  some  difference  in  the  transient  response.  Note  in  Fig.  64  that 
the  response  is  worse  with  limiters  than  without.  The  considerably  poorer  performance  in 
the  vicinity  of  100  ps  is  partly  caused  by  increased  a,-  values  and  partly  by  the  effects  of 
the  slowly  increasing  we'ght  magnitudes. 


The  steady-state  adapted  patterns  for  the  above  cases  are  not  changed  significantly 
by  the  limiter  modification,  provided  that  fx'0  >  1.  However,  the  transient  pattern  be¬ 
havior  will  usually  be  different  because  of  different  transient  decay  factors  af  and  a 
modification  in  the  eigenvector  beam  summation.  To  bring  out  these  differences,  let  us 
rewrite  'he  steady-state  orthonormal  weights  of  Eq.  (4.180)  in  the  foliowing  form: 


(4.181) 


where 


A  +  Mb 


(4.182) 


ML  - 


(4.183) 


fl  +  >  P, 


Under  the  eigenvalue  condition  that  =  (302,  n'im  represents  the  minimum  value  of  /rj 
from  Eq.  (4.178),  and  it  will  generally  be  much  smaller  than  r'0.  For  example,  in  Case  A 
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Fig.  62— Transient  response  with  and  without  limiter  modification  for  two  sources,  Case  C 
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Fig.  63— Transient  response  with  and  without  limiter  modification  for  single  source  with 

16%  bandwidth,  Case  E3 
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Fig.  64— Transient  response  with  and  without  limiter  modification  for  continuous-distribution 

sources.  Case  H 


=  2.83,  and  this  results  in  a  value  of  a  =  26.4  for  the  constant  defined  in  Eq.  (4.182). 
The  constant  a  represents  the  increase  in  the  magnitude  level  of  the  nonunique  ortho¬ 
normal  weights. 

Substituting  this  new  expression  for  (°°)  into  Eq.  (4.48),  we  can  manipulate  it 
into  the  following  relationship  for  the  orthonormal  weights: 


-a,t  (Pi  ~  O 


#»;.  +  i 


K  -  ..(4.184) 


=  a\V/qi  -  (1-e  "'•){ 

Compare  Eq.  (4.184)  against  Eq.  (4.65),  and  note  that  the  first  two  terms  are  the  same  if 
we  assume  that  n'jm  plays  the  same  roll  as  p0-  However,  we  have  an  extra  transient  term 
on  the  right-hand  side  in  the  above  expression. 

Proceeding  with  the  same  arguments  used  in  developing  Eq.  (4.119)  results  in  the 
following  output  pattern  function,  from  Eq.  (4.184): 


G{0,t)  =  a  <G 


,(0)  -  £  a  -  «*“'')  (^77Y")^g,(0)  -  s'm!  (4.185) 
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»•'<*,  o  -  (i  -  if)  I!  ,*((«)• 


(4.186) 


Compare  Eq.  (4.185)  against  Eq.  (4.119),  and  note  that  except  for  g'(9,t)  the  expression 
is  similar  if  we  assume  equivalence  of  n'im  and  Hq,  and  that  under  steady-state  conditions 
the  results  should  be  virtually  identical.  This  explains  why  the  steady-state  adapted  pat¬ 
terns  are  not  changed  significantly  by  the  limiter  modification,  except  for  the  increase  in 
level  by  the  constant  a. 

However,  the  transient  behavior  is  modified  by  the  extra  pattern  function  g'(0,t), 
which  has  some  rather  interesting  characteristics  when  plotted  over  time.  For  example, 
note  that  at  time  t  =  0, 


*■«<*(*) 


(4.187) 


i.e.,  it  reduces  to  the  quiescent  pattern  and  results  in  G(0, 0)  =  Gq(9),  as  it  should.  If  we 
select  t  =  Tu ,  at  which  time  all  of  the  unique  transient  decays  have  just  finished,  then 
all  the  remaining  a;-  will  be  identical  and  equal  to  a'0,  such  that 


Xn 

*'<*.<>  -(i  E  v«0> 


(4.188) 


where  ar'0  =  (1  +  n'im  )/t'0  and  Xlt  X2,  ■  ■ . ,  Xn  denote  the  nonunique  members  of  the  set 
of  eigenvector  beams.  But  such  a  summation  will  usually  be  approximately  eoual  to  the 
steady-state  adapted  pattern;  i.e., 


rt 


(4.189) 


Thus,  for  time  t>Tu,  Eq.  (4.185)  reduces  to  the  approximation 


>,t)  «  ajl  -  (l  -  G(9t< 


a)  t>Tu 


(4.190) 


where  we  have  the  steady-state  pattern  essentially  established  by  time  Tu  and  increasing 
slowly  in  level  thereafter  to  its  steady -state  magnitude  of  aG(0,°°). 

For  precise  pattern  calculations  during  the  transient,  one  should  work  with  complete 
Eqs.  (4.185)  and  (4.186).  This  requires  computing  a  completely  filled  Q-transformation 
matrix  and  using  all  the  eigenvector  beams  that  have  finite  W„j  weighting.  For  example, 
in  Case  H,  eigenvector  beam  g8(0)  is  not  used  in  Eq.  (4.119)  oecause  of  the  (ju,-  -  p.0) 
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term  deletion,  but  in  Eq.  (4.186)  it  is  prominent  and  helps  tr  account  for  the  longer 
convergence  time  noted  in  Fig.  64. 

In  summary,  the  hard-limiter  modification  results  in  the  same  dynamic  range  improve¬ 
ments  extolled  in  Sec.  3.4.  It  gives  superior  circuit  operation  and  loop  stability,  but  little 
improvement  in  transient  response  except  for  the  special  case  of  a  single  narrowband 
interference  source.  The  advantages  of  dynamic-range  improvement  are  enough  to  recom¬ 
mend  incorporation  of  the  hard  limiter  (or  fast  AGC). 


5.  CONCLUSION 

An  adaptive  array  consisting  of  ?.  '..'-element  linear  array  with  K  adaptive  control 
ioops  of  the  Applebaum  analog  type  derives  feedback  control  error  signals  from  the  corre¬ 
lations  between  element  signals,  on  the  basis  of  the  covariance  matrix  of  the  set  of  sys¬ 
tem  inputs.  The  covariance  matrix  is  a  summation  or  repository  of  all  interference  dis¬ 
tribution  information  as  seen  by  the  array  in  its  operating  environment. 

The  set  of  linear  differential  equations  associated  with  the  control-loop  network  can 
be  solved  by  a  Q-matrix  transformation  into  orthonormal  eigenvector  space,  with  the  Q 
matrix  consisting  of  ;,he  eigenvectors  of  the  covariance  matrix.  In  this  report,  the  purely 
mathematical  Q-matrix  transformation  has  been  interpreted  in  terms  of  orthogonal  beam- 
forming  networks,  similar  in  principle  to  a  Butler  matrix  beam-forming  network.  This  is 
the  basis  for  an  equivalent  “orthonormal”  adaptive  control-loop  network  which  is  much 
easier  to  understand  than  the  real  network  itself. 

The  Q-matrix  beam-forming  network  produces  a  set  of  K  orthogonal,  normalized 
eigenvector  beams,  in  which  the  array  element  weights  associated  with  each  eigenvector 
beam  consists  of  the  components  of  the  eigenvector.  The  output  powers  from  these  beams 
are  proportional  to  the  associated  eigenvalues  and  are  decorrelated.  The  eigenvector  beams 
have  been  used  in  a  convenient  expression  for  the  output  pattern  function  for  the  array. 

As  developed  in  Eq.  (4.119),  the  expression  is 

G(0,t)  =  Gq{d)  -  £  (i  -e'a,r)(jr^)K^(6)- 

This  requires  only  the  quiescent  pattern  Gq(6)  plus  the  unique  eigenvalues  and  associated 
unique  eigenvectors.  The  term  (nt  -  Mo)  eliminates  all  of  the  eigenvector  beams  whose 
output  power  consists  of  receiver  noise  only.  The  expression  uses  the  concept  of  array 
degrees  of  freedom,  since  one  degree  of  freedom  must  be  consumed  for  each  unique 
eigenvalue  generated  by  the  covariance  matrix. 

Performance  characteristics  were  calculated  from  the  above  expression  for  eight  dif¬ 
ferent  distributions  of  interference  sources,  as  listed  in  Table  2.  These  distributions  were 
selected  to  demonstrate  the  effects  of  source  power  level,  source  location  with  respect 
to  the  quiescent  beam  pattern,. source  spacing  in  terms  of  array  resolution,  source  band¬ 
width,  and  continuous  source  distributions.  Some  of  the  major  effects  noted  were 
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1.  Sidelobe  jamming  distorts  the  output  pattern  in  proportion  to  the  number  of 
unique  eigenvalues,  i.e.,  the  array  degrees  of  freedom  consumed. 

2.  Main-lobe  jamming  produces  severe  pattern  distortion. 

3.  Slow  convergence  rates  are  associated  with  small  eigenvalues,  and  small  eigen¬ 
values  are  generated  by  closely  spaced  sources,  continuous  source  distributions,  sources 
with  bandwidth,  low-power  sources,  and  combinations  of  these. 

4.  Separated  narrowband  sources  usually  consume  one  array  degree  of  freedom  per 
source,  but  sources  with  bandwidth  may  consume  two  or  three  array  degrees  of  freedom 
per  source,  depending  on  the  source  locations  and  their  percentage  bandwidths. 
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5.  Narrowband  sources  are  usually  attenuated  by  developing  a  pattern  null  close  to 
their  locations.  These  nulls  typically  do  not  aline  exactly  with  source  position  unless  the 
source  is  very  strong,  with  a  power  level  of  30  dB  or  more  above  receiver  noise  power. 

6.  If  there  ate  many  sources  closely  spaced  so  as  to  form  a  continuous  distribution, 
the  system  attenuates  them  by  developing  a  “low  sidelobe  notch”  region  in  the  adapted 
pattern;  i.e.,  the  array  handles  them  on  a  resolution  basis.  This  permits  the  array  to  use 
its  degrees  of  freedom  so  efficiently  that  it  can  handle  a  large  number  of  sources,  far  in 
excess  of  its  number  of  available  degrees  of  freedom. 

In  computing  the  eight  source-distribution  cases,  the  unique,  retrodirective  eigen¬ 
vector  beams  were  plotted.  It  was  found  that  they  could  be  characterized  by  a  family  of 
harmonically  related  pattern  “modes.”  Case  G,  in  particular,  produced  a  classic,  har¬ 
monically  related  set,  up  through  the  sixth  harmonic.  Even  though  these  beams  are  not 
in  real  space,  they  give  one  an  excellent  perspective  for  the  synthesis  of  adaptive  pattern 
reactions  to  given  interference  distributions.  They  should  be  useful  in  developing  beam- 
transformation  algorithms  and  techniques  for  future  adaptive  systems. 

Applebaum  introduced  a  hard-limiter  modification  into  his  basic  control  loop  to  im¬ 
prove  the  circuit  characteristics,  particularly  the  system  dynamic  range.  This  modifica¬ 
tion  was  included  in  the  analysis  and  resulted  in  a  major  change  in  the  equiva'ent  servo 
gain  factors,  plus  the  addition  of  an  extra  transient  term  in  the  above-mentioned  output 
pattern  function  for  the  array. 

One  effect  of  the  modification  is  that  the  magnitude  levels  of  the  adapted  weights 
change  considerably  but  slowly,  so  that  there  is  no  significant  difference  in  the  steady- 
state  adapted  pattern  except  for  the  increase  in  level.  Also,  the  transient  behavior  is 
modified  somewhat,  but  not  by  any  large  factor  except  in  the  very  special  case  of  a 
single  narrowband  jammer  of  low  power  level,  for  which  a  much  faster  response  is 
obtained. 

The  present  state-of-the-art  in  adaptive  arrays  is  still  in  its  infancy,  and  many  prob¬ 
lems  remain  to  be  addressed,  in  theory  and  in  practice.  However,  progress  is  rapid  and 
the  future  for  these  systems  looks  very  promising.  For  example,  accelerated  convergence 
techniques  are  being  developed  to  overcome  sluggish  response  situations,  and  main-beam 
constraint  techniques  are  being  developed  to  overcome  the  objectionable  main-beam  dis¬ 
tortions  suffered  by  most  of  the  current  adaptive  techniques. 
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Appendix  A 

TRANSIENT  RESPONSE  OF  A  SIMPLE  RLC  CIRCUIT 


From  the  circuit  shown  in  Fig.  Al,  the  integro-differential  equation  may  be  written 


f  +  cf  ♦ 


L  T 


Fig.  Al  -Single-tuned  RLC  filler  circuit  for  W 


Taking  the  Laplace  transform  of  both  sides  and  assuming  that  initial  charges  and  currents 
are  zero,  we  get 


(k  * cs  +  s)  W(s>  - 


where  3  =  a  +  jio  is  the  usual  Laplace  complex  variable.  Then 


T(S)V(S), 


where 


Y(S)  = 


(S  +  a)2  +  02 


.  p2  =  Pi  -  «2 


natural  resonant  frequency 


-r  =  2RC  =  time  constant 


|  -'.3 
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1  1 

a  =  —  =  777:7,-  =  attenuation  constant . 


r  2  RC 


(A8) 


y(S)  is  the  Laplace  transform  expression  for  the  circuit  itself,  independent  of  any  excita¬ 
tion  function. 


and 


If  we  now  take  an  input  step  function  of  a  sinusoid  for  0,  so  that  v  =  0  for  t  <  0 


v  =  sin  (^ot  +  ^0  “  u)  fo r  t  >  0 , 

then  its  Laplace  transform  will  be  found  to  be 

r  r  mix  *n 

j  sin  (0 0  -  u)  IS  +’ 


(A9) 


m  • 


From  Eq.  (A3)  we  may  write 


ms)  » 


tan  (0o  ~  u)J 


S2  +  coo2 


(A10) 


CJq 

abG 

2aS  1  < 

sin  (0o  u)  S  +  ten  Wo  _  u) 

2 

(S  +  a)2  +  /32J 

S2  +  J02 

(All) 


This  may  be  rewritten  in  the  following  series  form,  based  upon  the  roots  of  S 


W(S)  =  abG 


2  [(S  +  of)  -  $  (S  +  a)  t-  Jl 3J 
Solving  for  the  constants  0,  B,  0,  and  A,  wo  have 


0  .=  0o 


-  u  +  tan  1  \-  £«coo  / 


B  = 


=  ~  +  tan-1 


■  (PV  *  ten-1  r. 

_a{/302  +  coq-)]  I' 


*  j 


I 
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Appendix  B 

CONSTRUCTION  OF  NONUNIQUE  EIGENVECTORS  FOR 
FILLING  OUT  THE  Q  MATRIX 


The  method  for  filling  out  the  Q  matrix  is  based  on  using  controlled  null  placement 
to  optimize  spatial  coverage,  in  a  manner  similar  to  the  operation  of  the  Davies  tree  net¬ 
work  discussed  in  Sec.  4.3.  The  vectors  so  derived  are  then  orthogonalizcd  by  the  Grara- 
Schmidt  procedure  and  normalized  to  produce  the  desired  eigenvectors. 


Using  Case  B  of  Table  Bl,  which  involves  two  sources  located  at  18  and  22  degrees 
respectively,  we  establish  these  tv,  o  source  locations  on  the  Z-plane  unit  circle  as  shown 
in  Fig.  Bla,  using  Eq.  (4.80):  Z  =  exp  (;2u).  These  shall  be  regarded  as  constrained  null 
locations  for  the  nonunique  beams.  It  is  necessary  to  have  as  many  constrained  nulls  as 
there  are  unique  eigenvalues;  we  shall  denote  this  number  by  the  symbol  Ku.  With  Ku 
null  positions  fixed,  there  remain  ( K  -  Ku)  null  positions  to  be  located  on  the  circle  end, 
although  the  choice  of  these  positions  is  theoretically  arbitrary,  it  seems  desirable  from  a 
beam-forming  point  of  view  to  seek  locations  which  optimize  the  spatial  coverage  by 
maximizing  the  separations  between  the  remaining  null  locations.  If  this  principle  is  ap¬ 
plied  in  Fig.  Bib,  the  remaining  null  locations  should  be  spaced  apart  equally  by  the  angle 


f  360-0  \  = 
\K  ~  Ku  +  i/ 


49.7° 


where  0  is  the  angular  separation  of  the  two  closely  spaced  constrained  nulls  on  the  circle. 


(a)  Coni  trained-null  (or  source 
position)  locations 


Fif.  Bl— Null  locations  on  Z-plane  unit 
circle  for  Case  B 
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(b)  Locations  chosen  for  six  remaining  nulls 

Fig.  Bl—  Null  locations  on  £-plane  unit 
circle  for  Case  B  (Continued) 


When  all  of  the  null  locations  L>i  through  Lg  have  been  selected,  the  next  step  is  to 
form  the  nonunique  beams  by  choosing  (K-  1)  of  the  K  null  positions  for  use  in  Eq. 
(3.83):  F(Z)  =  -  Z\)(Z  -  Z2){Z  -  Za)  ...  ( Z  -  Zjc-i).  The  constrained  null  posi¬ 

tions  L2  and  L3  must  be  used  in  each  beam.  If  the  ith  beam  is  denoted  by  G;(2),  then 

Gi(Z)  =  (Z  -  Zi)(z  -  z2yz  -  z3) ...  (Z -  ZK.X) .  (Bi) 

The  null  selections  for  the  six  beams  associated  with  Fig.  Bl  are  shown  in  Table  B2. 
Computer  programs  are  available  for  using  the  null  positions  of  Eq.  (Bl)  to  solve  for  the 
complex  coefficients  associated  with  the  standard  polynomial  form  of  G,(Z), 

G,(Z)  =  Aq  +  A\Z  +  A2Zz  +  Aa% 3  +  ...  +  Ak-\Zk~ 1 .  (B2) 

These  complex  coefficients  constitute  the  components  of  the  element  weighting  vectors 
associated  with  the  beams: 


K- 1 

G/(0)  =  {W /S)  =  £  AkSk+1 ,  (B3) 

k=0 

where 

W/  =  [A0,AltA2,Aa . Ajc-J  (B4) 

and  S  is  defined  in  Eq.  (4.8):  S'  =  [Si,S2, S;3,  •••. SK ] .  In  general,  the  W,-  vectors 
calculated  by  the  above  procedure  are  not  orthogonal  to  one  another,  and  it  is  necessary 
to  orthogonalize  by  the  Gram-Schmidt  procedure*. 


*F.B.  Hildebrand,  Methods  of  Applied  Mathematics,  2d.  ed.,  Prentice-Hall,  Inc.,  Englewood  Cliffs.  N.J. 
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Table  B1 

Selected  Distributions  of  Interference  Sources 


Interference  Sources 


Power 

Ratios 


1,250 

1,250 

1,200 

1,250 

125 

40 

125 

400 

1,250 

1,250 

1,250 

1,250 


Location 

Angle 

(Deg) 


1,250 

1,100 

1,100 

1,100 


1,100 

1,100 

1,100 

1,100 

1,100 

1.100 


Band¬ 

width 

(%) 


36 

26 

11 

48 

21 

11 

66 

19 

11 

■66 

19 

11 

■48 

21 

11 

-3o 

26 

11 

36 

26 

11 

48 

21 

11 

66 

19 

11 

Unique  Eigenvalues 

Number 

Values 

1 

10,001 

2 

18,544 

1,058 

2 

10,812 

190 

4 

11,616 

2,486 

406 

16.5 

1 

10,001 

2 

9,986 

16.4 

3 

9,529 

469 

4.7 

1 

10,001 

6 

13,316 

9,692 

3,091 

296 

10.5 

1.16 

8 

13,532 

13,386 

12,619 

9,682 

3,224 

350 

14.7 

1.18 

*  .TO,. 
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Table  B2 

Beam  Null  Locations 


Null 

Nonunique  Beams 

Positions 

Gi(6) 

02(A) 

O3(0) 

G4(0) 

GS(A) 

Oe(A) 

Z\ 

LZ 

' — - -  - 

Ll 

1 -  - 

Lx 

Lx 

Lx 

Z2 

Lz 

LZ 

L‘X 

Lo 

1 

l2 

L2 

Zz 

l4 

Lz 

Lz 

Lz 

Lz 

Lz 

Z4 

l6 

*5 

La 

La 

La 

La 

Z<3 

Lq 

l6 

Lz 

l5 

Lz 

Lz 

Zq 

L'i 

L'I 

l7 

L'i 

Lz 

Lz 

Z1 

Ls 

L» 

l8 

Ls 

Lz 

L'i 

In  our  example,  we  start  with  the  two  unique  eigenvectors  e4  and  e2 ,  the  compo¬ 
nents  of  which  are  printed  out  in  Table  B3.  These  two  are  already  orthogonal  and  normal¬ 
ized,  of  course.  We  next  take  the  first  vector  Wi ,  corresponding  to  the  nonunique  beam 
Gi(0),  and  form  the  new  vector, 

V3  =  Wi  -  qiei  -  q2e 2.  (B5; 

The  requirement  that  V3  be  orthogonal  to  ei  leads  to  the  relation 

i4%)  =  {**%)  -  91  =  0 


similarly, 


91  =  (e^Wt); 

92  =  (eJ'Wj) . 


Thus,  we  “subtract  off”  the  e4  and  e2  components  of  Wj,  obtaining  the  vector  Vj,  which 
is  now  orthogonal  to  both  ei  and  e2.  V3  is  then  normalized  by  dividing  by  its  Hermitian 
length  to  get  the  eigenvector  e3 : 

C3='ptr  (B8) 

Eigenvector  e3  is  therefore  our  first  nonunique  eigenvector. 

To  get  the  second  nonunique  eigenvector  e4,  we  take  the  second  vector  W2  from  our 
set  and  form  the  new  vector  V4,  orthogonal  to  each  eigenvector,  so  that 


imsuxa  owaM 
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Table  B3 

Eigenvalues  and  Eigenvectors  for  a  Two-Source  Case*! 


Eigenvector 

Coefficient 

Amplitude 

Phase  Angle 
(Deg) 

1,1 

0.340120 

69.7477 

1,2 

0.351501 

8.3626 

1,3 

0.359156 

-53.03 

1,4 

0.363004 

245.573 

1,5 

0.363004 

184.175 

1,6 

0.359156 

122.778 

1,7 

0.351501 

61.3852 

1,8 

0.340120 

0 

2,1 

0.537157 

250.308 

2,2 

0.387808 

188.701 

2,3 

0.234347 

126.989 

2,4 

0.078413 

64.5464 

2,5 

0.078413 

185.762 

2,6 

0.234347 

123.319 

2,7 

0.387808 

61.607 

2,8 

0.537157 

0 

*Jatnmer  power  ratios  1,250;  location  angles  18°  and  22°;  band¬ 
width  0%;  one  spectrum  line. 

^Unique  eigenvalues:  18,544.4  and  1,057.,58 


v4  =  w2  -  7iei  -  <22e2  -  93e3  ,  (B9) 

where  the  qi  constants  are  to  be  evaluated  as  before.  Dividing  V4  by  its  Hernrtian  length 
then  results  in  e4. 

The  above  process  is  simply  repeated  for  each  of  the  W;  in  turn,  until  all  of  the  non¬ 
unique  eigenvectors  have  been  computed,  whereupon  they  can  be  entered  into  the  Q 
matrix  to  fill  it  out. 


